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ABSTRACT 


Future  fleet  needs  for  real-time  intelligence  require  an  unmanned 
platform  capable  of  operations  from  small  surface  combatants 
without  the  need  for  extensive  support  personnel  or  equipment  and 
without  causing  disruption  to  the  operations  of  the  ship  from  which 
it  operates.  A  candidate  must  not  only  takeoff  and  land  vertically 
but  also  be  capable  of  high  forward  flight  speeds  and  efficient  on- 
station  performance.  The  design  and  initial  fabrication  of  a  Vertical 
Attitude  Takeoff  and  Landing  (VATOL)  Unmanned  Air  Vehicle  (UAV) 
airframe  was  completed  at  the  Naval  Postgraduate  School.  The 
vehicle,  called  Archytas,  was  a  combination  of  two  existing  UAV's-- 
the  AROD  and  Aquila—as  well  as  locally  manufactured  components, 
including  a  canard  support  structure  and  wing  spar.  The  objective  of 
creating  Archytas  was  to  provide  a  proof-of-concept  platform  for 
research  to  explore  performance  trade-offs  and  stability 
augmentation.  A  three-degree-of-freedom  simulation  was  used  as 
the  focus  of  the  design  efforts,  to  validate  design  decisions  made  in 
the  fields  of  propulsion,  aerodynamics,  structures  and  flight 
mechanics.  Engine  tests  were  conducted  to  determine  thrust  and 
control  power.  Structural  components  were  designed,  fabricated  and 
then  tested,  making  modifications  where  necessary  to  ensure 
sufficient  airframe  strength.  A  longitudinal  control  system  was 
designed,  validated  by  simulation,  and  tested  structurally. 
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1.  INTRODUCTION 


A.  MISSION  NEED 

Among  the  many  news  stories  reported  during  Desert  Storm  was 
that  of  a  group  of  Iraqi  soldiers  waving  a  white  flag  to  a  Unmanned 
Air  Vehicle  (UAV).  This  particular  UAV,  the  Pioneer,  had  been 
providing  targeting  information  to  a  battleship  firing  16  inch  shells 
from  over  the  horizon.  As  the  Pioneer  helped  refine  the  battleship's 
aim,  the  Iraqis  knew  that  the  small  craft  which  flew  over  their 
position  and  the  ship  raining  metal  upon  them  were  a  team.  Lacking 
any  direct  contact  with  those  operating  the  big  guns,  they 
surrendered  to  the  small  package  of  hardware  and  electronics 
orbiting  above  them.  The  utility  of  the  UAV,  long  known  to  military 
commanders,  was  graphically  demonstrated  on  world-wide  TV. 

UAV’s,  formerly  known  as  Remotely  Piloted  Vehicles  (RPV’s)  can 
serve  military  commanders  in  a  variety  of  missions.  C''rrying  optical 
and/or  infrared  (IR)  sensors,  they  provide  targeting,  intelligence  and 
reconnaissance  information.  They  can  also  be  equipped  with 
communications  intelligence  (COMINT)  and  electronics  intelligence 
(ELINT)  packages  for  relay  to  a  ground  or  sea-based  user.  Other 

systems  that  can  Le  carried  aloft  by  UAV's  include;  radar, 
communications  relay,  electronic  countermeasures  (ECM),  laser 
designators,  nuclear/biological/chemical  (NBC)  contamination 
sensors,  and  even  weapons,  all  without  sending  military  personnel  in 


B.  LIMITATIONS  OF  CURRENT  SYSTEMS 


The  utility,  then,  for  the  UAV  is  great  but  there  are  limitations. 
Many  systems,  such  as  the  Pioneer,  Sky  Owl,  and  the  Aquila  shown  in 
Figure  1.1,  require  special  equipment  or  procedures  for  launch  and 
recovery.  This  equipment  includes  pneumatic  rail-launchers  or 
rocket  assist  for  takeoff,  and  nets  for  capture  at  the  completion  of  the 
UAV's  mission.  The  landborne  systems  require  a  fairly  large  area  for 
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Figure  1.1 

Aquila  UAV  Sliown  Wiilt  Special  Launch  and  Recovery  Lquiptnent 
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takeoff  and  landing.  Seaborne  UAV's  can  impact  the  battlegroup’s 
operations  due  to  stringent  wind  and  sea  state  envelopes  for  launch 
and  recovery.  Because  of  these  disadvantages,  other  concepts  have 
been  explored  as  alternatives  to  the  classic  fixed-wing  UAV, 

Rotary  wing  UAV's,  such  as  unmanned  helicopters  and  ducted- 
fan  types,  have  the  advantage  that  they  can  operate  from  small 
forward-deployed  areas  or  ship  decks  without  special  equipment. 
However,  systems  with  exposed  rotors  present  a  formidable  safety 
hazard  to  nearby  personnel.  Additionally,  vertical-only  UAV's  are 
limited  by  their  relatively  slow  forward  flight  speeds  and  generally 
shorter  endurance  times  once  on-station.  They  simply  cannot  get  to 
the  target  area  fast  enough  or  stay  there  long  enough. 

C.  MOTIVATION 

Another  class  of  UAV's  that  has  recently  received  renewed 
attention  are  those  with  vertical  takeoff  and  landing  (VTOL) 
capability  and  the  ability  to  transition  to  wing-borne  forward  flight. 
These  vehicles  conceptually  have  the  strengths  of  both  fixed  and 
rotary  wing  vehicles:  able  to  takeoff  and  land  vertically  but  with 
higher  forward  flight  speeds.  This  class  of  aircraft  is  comprised  of 
several  types,  including:  lift  fan  (eg:  Ryan  XV-5),  tilt  rotor  (eg:  Bell 
XV-15  and  V'22),  tilt  wing  (eg:  LTV  XC-142A),  vectored  thrust  (eg: 
Harrier).  Lift  engines  plus  cruise  engines  (separate  engines  for  each 
function,  eg:  Dasault  Mirage  IIIV),  and  Vertical  Attitude  Takeoff  and 
Landing  (VATOL),  also  known  as  "Tailsitters"  [Ref.  1]. 
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The  goal,  then,  for  the  next  generation  of  UAV’s  will  be  to  retain 
the  capabilities  of  current  systems  but  to  provide  this  capability  in  a 
more  flexible  form. 
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II.  BACKGROUND 


A.  VERTICAL  ATTITUDE  TAKEOFF  AND  LANDING  (VATOL) 
1.  General 

As  discussed  in  the  introduction,  one  of  the  VTOL  concepts 
capable  of  wing-borne  flight  is  the  VATOL  type.  It  is  this  type  which 
is  the  subject  of  this  thesis.  In  1949  the  first  "convertible  aircraft 
congress"  was  held  and  from  this  meeting  came  several  designs,  as 
depicted  in  Figures  II.  1  and  II. 2  [Ref.  2].  Some  of  the  advantages  of 
the  VATOL  design  [Ref.  3]  are: 

1)  Simplicity  of  design.  The  same  engine  is  used  for  vertical 
lift  and  forward  flight,  normally  with  only  simple  modifications  to 
allow  hovering  control.  This  simple  design  generally  results  in  less 
weight  penalty  because  there  are  no  pivoting  wings,  nozzles  or 
engine  nacelles  with  the  accompanying  structural  weight. 

2)  High  thrust-to-weight  ratio  for  mission  use.  VATOL 
aircraft,  with  lower  structural  weight,  have  more  excess  thrust  which 
can  be  utilized  during  takeoff  and  landing  as  well  as  in  the  mission 
area. 

3)  Reduced  or  eliminated  degradation  of  vertical  lift 
capability  due  to  "suckdown"  or  exhaust  ingestion. 

4)  In  general,  better  control  at  high  angles  of  attack. 

In  a  1978  study  of  several  generic  designs,  the  David  Taylor 
Research  Center  [Ref.  4]  concluded  that  VATOI,  was  the  best  VTOL 
configuration,  with  a  12%  higher  payload/operating  weight  ratio  than 
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Figure  11.2;  Early  VATOL  concept 
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the  next  best  configuration,  as  shown  in  Figure  11.3. 

The  disadvantages  of  the  VATOL  concept  in  general  include: 

1)  In  the  case  where  special  launch-and/or-recovery 
platforms  or  equipment  are  required,  there  is  a  reduction  in  the 
operational  flexibility  compared  with  conventional  fixed  wing 
aircraft. 

2)  Most  VATOL  designs  have  no  STOL  capability,  limiting 
their  ability  to  handle  increased  takeoff  weight. 

2.  VATOL  historical  designs 

During  the  1950's  and  early  1960's  and  then  again  during 
the  Carter  administration,  there  was  much  research  on  VATOL 
designs,  including  several  successful  manned  and  unmanned 
programs.  The  prop-driven  Convair  XFY-1,  Figure  11.4,  was  a 
successful  design  flown  in  1953.  Over  400  flights  were  made,  all 
with  transitions  to  and  from  horizontal  flight.  Later  came  the  jet- 
powered  Ryan  XV-13,  Figure  II.5,  which  completed  136  flights,  of 
which  104  were  VATOL  [Re'.  1].  Both  programs  were  completed 
without  serious  incident.  It  is  important  to  note  the  configurations  of 
these  experimental  aircraft:  all,  except  for  the  least  successful 
Lockheed  XFV-1  about  which  less  has  been  written,  had  highly- 
swept  delta  wings  and  most  had  reaction  control  augmenting 
aerodynamic  control.  A  review  of  the  many  documents  concerning 
these  and  other  VATOL  studies  revc  led  several  potential  problem 
areas  for  VATOL  designs.  These  include: 

1)  Longitudinal  trim.  The  jet  designs  all  suffered  from  a 
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VTOL  COMPARISON 


VATOL  TILT  L+L/C  L/C 
TYPE  VTOL 


Figure  II. 3 

Payload  coniparision  of  several  VTOL  types. 
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Figure  11.4:  CONVAIR  XFY-1 -Propeller-Driven  VATOL  Research 

Airplane 


significant  pitch-trim  problem,  encountered  during  constant-altitude 
decelerating  transitions.  According  to  Girard  and  Everett  [Ref.  5],  the 
XV- 13  "could  not  be  flown  steady  state  at  attitude  angles  between  32 
and  70  deg  due  to  an  incompatibility  between  thrust  required  for 
constant  altitude,  and  that  required  for  longitudinal  control,  the 
latter  requirement  being  the  greater."  The  XV- 13  utilized  a  gimbaled 
nozzle,  much  like  that  on  a  rocket,  to  vector  thrust.  The  solution  to 
the  problem  was  to  simply  pitch  through  this  region  by  "following  a 
memorized  schedule  of  engine  rpm  versus  airplane  angle"  (to  avoid 
excessive  altitude  gain)  [Ref.  5]. 

2)  Longitudinal  flying  qualities.  Because  the  VATOL  design 
is  often  operated  at  relatively  high  angles  of  attack,  it  is  desirable 
that  the  vehicle  have  a  high  coefficient  of  lift  and  a  post-stall  region 
that  is  as  free  of  abrupt  changes  in  the  aircraft’s  stability  derivatives 
as  possible.  Longitudinally,  this  means  desiring  benign  stall 
characteristics.  The  typical  approach  has  been  to  utilize  fairly  low- 
aspect-ratio  delta  wings,  trading  off  a  lower  coefficient  of  lift  for 
good  high-alpha  flying  qualities.  Not  all  programs  experienced 
longitudinal  flying  qualities  problems,  however.  In  1961  NASA 
performed  tests  on  a  1 /5-scale  RPV  very  similar  in  configuration  to 
the  XV-13.  A  three-view  of  the  model  is  shown  in  Figure  II. 6.  The 
report  on  this  vehicle  stated:  "Transitions  from  hovering  to  normal 
forward  flight  and  back  to  hovering  flight  could  be  made  smoothly 
and  easily  in  the  full-scale  tunnel,  and  the  model  seemed  to  have 
stability  of  angle  of  attack  over  most  of  the  speed  range.  These 


flights.. .represented  slow,  constant-altitude  transitions  and  covered  a 
range  of  angles  of  attack  from  about  20  deg  to  90  deg."  {Ref  6]. 

3)  Lateral/directional  flying  qualities.  Another  problem 
encountered  in  the  high  angle  of  attack  regime  is  marked  changes  in 
aircraft  lateral -directional  flying  qualities.  In  many  airplanes  at  AOA's 
just  a  few  degrees  above  stall,  the  sign  of  directional  stability  can 
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Figure  II.6 

NASA  1/5  Scale  Model  of  a  Jet  Powered  VATOL  Research  Airplane 


change  and  aerodynamic  roll  control  power  can  be  drastically 
reduced.  The  XV- 13  experienced  this  problem.  Roll  control  required 
to  "handle  extreme  sideslip  angles  exceeds  the  control  available  in 
level  flight  by  a  considerable  margin  at  and  just  beyond  the  stall. 
Adequate  control  is  available  at  full  power  [the  XV- 13  utilized 
reaction  roll  control  to  augment  aerodynamic  control]  but  then  the 
airplane  would  begin  to  climb."  [Ref  5]. 

The  preceeding  review  of  flight  test  reports  for  both 
manned  and  unmanned  VATOL  research  projects  suggests  an 
additional  design  objective:  to  design  the  vehicle  to  be  flown  with  an 
angle  of  attack  (a)  in  regions  of  predictable  airplane  response.  This 
means  either  keeping  the  a  low— not  more  than  a  few  degrees 
beyond  stall— or  exploring  higher  a  behavior  in  wind  tunnel  or  flight 
test. 

B.  NAVAL  POSTGRADUATE  SCHOOL  UAV  PROGRAMS 
1.  General 

UAV  research  at  the  Naval  Postgraduate  School  (NPS)  is 
conducted  by  the  UAV  Flight  Research  Laboratory  (FRL).  Members 
of  this  lab  have  experience  in  several  UAV  projects  and  have 
obtained  assets  for  use  in  future  proof-of-concepi  designs. 
Additionally,  the  FRL  has  facilities  to  manufacture  and  test  UAV's. 
Resident  skills  include  metalworking,  composite  manufacturing  and 
model  building. 


2.  Aquila  UAV 

Developed  by  the  U.S.  Army  in  the  1970's  as  a  surveillance 
platform,  the  Lockheed  MQM-105  is  a  flying  wing  design  with  a 
shrouded  pusher  propeller,  as  shown  in  Figure  11.7.  The  Aquila  was 
made  from  Kevlar  composite  material  which  resulted  in  a  strong, 
lightweight  airframe.  The  Aquila  was  launched  by  pneumatic  ra  I 
and  captured  in  a  net,  as  shown  in  Figure  I.l.  Lateral  and 
longitudinal  flight  control  was  achieved  through  wing-mounted 
elevens.  The  Aquila  had  no  sustained  vertical  flight  capability. 
Several  Aquila  airframes  were  obtained  by  the  FRL  for  use  in  UAV 
research. 

3.  Airborne  Remotely  Operated  Device  (AROD) 

Designed  and  fabricated  by  Sandia  Research  Laboratory 

under  a  contract  for  the  Naval  Ocean  Systems  Center  (NOSC),  the 
AROD,  shown  in  Figure  II. 8,  was  a  ducted-fan  UAV  built  of  graphite 
and  Kevlar  composites.  It  was  powered  by  a  two-cylinder,  two-cycle 
engine  developing  approximately  26  BHP.  The  engine  is  more  fully 
described  in  Section  IV.  Capable  of  hovering  and  limited  forward 
flight,  the  AROD  was  controlled  by  four  vanes  mounted  in  the  wake 
of  the  engine  driving  a  three-bladed  composite  propeller. 
Stabilization  was  performed  by  an  onboard  computer  fed  by  rate 
sensors  and  gyros.  AROD  characteristics  are  shown  in  Table  ILL  All 
remaining  AROD  assets  were  obtained  in  1993  by  the  FRL  for  UAV 
flight  research. 


4.  Archytas  UAV  Notional  Concept 

The  work  of  theses  by  Blanchette,  Ellwood,  Brynstead,  Merz 
and  Kress  [Refs.  7-11]  have  led  to  t^’e  development  of  a  VATOL  UAV 
concept  combining  the  assets  of  Aquila  and  AROD  with  locally 
manufactured  components.  Figure  11.9  shows  the  notional  concept 
for  FRL's  VATOL  UAV,  named  Archytas.  It  is  designed  to  be  flown  in 
a  hover  in  the  same  manner  as  the  AROD,  by  using  the  vanes 
immersed  in  the  wake  of  the  propeller,  and  then  transitioned  into 
forward  flight  and  flown  as  a  conventional  fixed-wing  airplane,  much 
like  the  Aquila.  Following  return  to  the  landing  site  it  will  be 
designed  to  transition  back  to  a  vertical  attitude  and  land  from  a 
hover. 

In  his  thesis.  Lieutenant  Kress  designed  the  spar  to  connect 
the  AROD  to  the  Aquila  wing,  designed  and  built  an  engine  test  stand, 
and  performed  initial  engine  tests  [Ref.  11]. 
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Figure  II. 7:  Lockheed/US  Army  MQM-105  Aquila  UAV 
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Figure  II.8:  Sandia  Lab/NOSC  AROD  Ducted-Fan  UAV 


Table  Il.I  AROD  CHARACTERISTICS 


Figure  II. 9:  NFS  VATOL  UAV  Notional  Concept,  the  Archytas 


111.  ARCHYTAS  DESIGN-GENERAL 


A.  DESIGN  OBJECTIVES 

1.  General 

The  objectives  of  the  Archytas  program  are  to  design, 
produce  and  test  a  proof-of-concept  VATOL  UAV.  This  vehicle  will 
be  used  to  examine  the  technologies  and  techniques  applicable  to  a 
military  UAV.  A  conceptual  follow-on  UAV  to  the  Archytas  program 
would  be  capable  of  carrying  a  variety  of  military  payloads.  A  key 
feature  of  this  future  military  UAV  must  be  the  ability  to  conduct 
operations  with  minimal  support  equipment  and  personnel. 
Additionally,  these  operations  must  be  conducted  without  adversely 
impacting  the  operations  of  other  components  of  the  war-fighting 
team. 

2.  Specific 

As  in  any  design  effort,  there  are  contraints  which  the 
designer  strives  to  meet.  In  the  case  of  the  Archytas,  the  primary 
specific  design  objectives  were: 

•  Minimum  weight.  To  increase  the  maneuver  margins 
during  hovering  flight,  high  thrust  to  weight  (TAV)  is  desirable. 
Clearly,  a  TAV  less  than  one  will  prevent  flight  for  a  VATOL  airplane. 
Historical  information  [Refs  12,  13,  14]  suggest  that  the  minimum 
acceptable  TAV  is  in  the  range  of  1.05  to  1.08  with  a  desirable  TAV 
of  1.15.  The  thrust  of  the  Archytas  was  considered  to  be  fixed.  In 


other  words,  the  AROD's  propulsion  system--the  duct,  inlet,  propeller 
and  engine--would  be  incorporated  into  the  Archytas  design  with 
few  modifications.  This  left  weight  as  the  design  variable--low 
weight  was  a  significant  design  goal. 

•  Low  angle-of-attack  flight.  As  a  result  of  the  review  of 
the  historical  VATOL  programs  another  design  goal  emerged:  that  of 
operating  the  Archytas,  if  possible,  at  maximum  angles  of  attack  of  a 
few  degrees  beyond  stall  when  at  flight  speeds  beyond  slow  hover. 
This  limit  on  angle  of  attack  was  to  ensure  that  the  problems  and 
uncertainties  of  high  angle  of  attack  flight  were  avoided,  most 
notably  in  the  area  of  longitudinal  and  lateral/directional  flying 
qualities. 

•  Structural  support  for  a  100  lb  vehicle  at  8  g’s  ultimate 
normal  load  factor.  This  requirement  resulted  from  the  work  of 
Kress  [Ref.  11], 

•  Minimum  modification  of  AROD  duct,  to  prevent  loss  of 
structural  integrity  of  the  AROD. 

•  Favorable  static  longitudinal  stability  for  forward  flight. 

•  Favorable  longitudinal  control  power  for  forward  flight. 

•  Minimum  special  equipment  to  launch  or  recover  the 

vehicle. 

•  Minimum  modification  of  Aquila  wing. 

B.  DESIGN  TEAM 

The  subject  of  this  thesis,  design  and  fabrication  of  the  Archytas 
airframe,  is  one  part  of  a  larger  effort  at  NFS  to  produce  a 


controllable  VATOL  UAV  capable  of  autonomous  flight,  as  shown  in 
Figure  Ill.l.  Research  is  ongoing  by  members  of  the  NFS 
Aeronautical  Engineering,  Electrical  Engineering  and  Weapons 
Systems  curricula  to  achieve  this  goal.  The  topics  of  the  research 
work  include: 

•  Control,  including  advanced  six  degree  of  freedom  simulation. 

•  Navigation;  sensors  and  integration. 

•  Communication  and  data  acquisition.  Up-link  and  down-link. 

Design  of  the  airframe,  therefore,  required  interaction  with  these 

other  disciplines.  The  areas  of  most  interest  to  all  parties  included: 
weight  and  center  of  gravity  management,  configuration,  vehicle 
performance,  aircraft  control  power  and  aircraft  flying  qualities. 

C.  DESIGN  PROCESS  -ARCHYTAS  AIRFRAME 

Design  is  by  its  nature  an  iterative  process.  Design  of  the 
Archytas  was  no  exception  to  this  rule.  Several  disciplines  interacted 
closely  in  the  airframe  design  process.  As  depicted  in  Figure  III. 2, 
the  disciplines  included: 

•  Aerodynamics  research  included  wing  placement  and 
orientation.  This  involved  a  strong  interaction  with  Flight  Mechanics 
to  achieve  proper  longitudinal  characteristics.  Additionally,  strong 
interaction  with  the  Structures  discipline  was  necessary  to  supply 
predicted  flight  loads. 

•  Flight  Mechanics  theory  was  utilized  for  proper  sizing  and 
placement  of  the  canard  control  surface. 


Figure  Ill.l 

Arcliytas  Mission  Profile 
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•  Simulation  was  used  to  predict  results  of  all  displines'  effects 
on  design. 

•  Structural  analysis  was  performed  to  achieve  a  lightweight 
structure  capable  of  supporting  predicted  aero-loads.  This  portion  of 
the  project  included  test  of  key  structural  components. 

•  Propulsion  tests  were  conducted  to  characterize  the  installed 
propulsion  system. 


CANARD  DESIGN 

SIZING  AND  PLACEMENT 
STRUCTURAL  DESIGN  AND  TEST 


n  I 

Figure  III. 2 

Archytas  Airframe  Design  Process 

2  1 


D.  THESIS  ORGANIZATION 

This  thesis  is  organized  by  subject  area  as  shown  in  the 
preceeding  figure.  A  chapter  is  devoted  to  each  of  the  following 
subjects: 

•  Section  IV:  Engine  Tests 

•  Section  V:  Wing  spar  Structural  Design,  Fabrication  and  test 

•  Section  VI:  Canard  Design  and  test 

•  Section  VII:  Simulation 

Sections  VIII  and  IX  follow  with  conclusions  and 
recommendations. 
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IV.  ENGINE  EXPERIMENTS 


A.  PURPOSE  OF  TESTS 

The  purposes  of  performing  engine  tests  included: 

•  Characterization  of  thrust  versus  throttle  position  for  use  in 
simulation  program's  equations  of  motion.  This  led  directly  to 
estimation  of  desirable  thrust  to  weight  ratios  for  vertical  takeoff 
and  landing— which  in  turn  led  to  a  directed  weight  management 
effort. 

•  Measurement  of  noise  and  comparison  of  thrust  with  and 
without  factory-installed  muffler  system. 

•  Determination  of  roll  control  power  for  use  by  other  design 
team  members.  These  data  were  also  used  to  model  pitch  control 
power  in  the  simulation  program. 

•  Characterization  of  engine  acceleration  and  deceleration. 

•  Determination  of  fuel  flow  rates. 

•  Survey  of  flow  velocity  at  AROD  duct  exit  for  possible  use  in 
control  power  modeling. 

•  Measurement  of  the  effects  on  thrust  of  deflecting  the  vanes- 
-called  "Blockage  effect". 

•  Characterization  of  "blowback  effect"— changes  in  vane 
position  caused  by  propeller  wash. 

•  Gaining  of  experience  in  starting  and  running  AROD  in 
anticipation  of  flight  test. 
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B.  DESCRIPTION  OF  TEST  EQUIPMENT 

1 .  Engine  description  and  installation 

The  engine,  manufactured  by  Herbrandson  Enterprises,  was 
a  two-cylinder,  two-stroke,  horizontally-opposed,  air-cooled  engine 
with  a  displacement  of  290  cc  (17.4  cu  in)  developing  approximately 
26  BHP,  driving  a  graphite-composite  three-bladed,  fixed-pitch 
propeller.  The  uninstalled  engine  with  its  ignition  electronics  and 
without  mufflers  weighed  15.5  pounds.  [Ref  15].  Installation 
characteristics  were  previously  stated  in  Table  II. I.  The  engine  was 
fueled  with  a  mixture  of  two-cycle  engine  oil  and  unleaded  supreme 
gasoline  at  a  ratio  of  80:1. 

The  propeller  blade  angle  could  be  varied  by  the  loosening 
of  a  blade  retaining  assembly.  Propeller  blade  angle  was  left  at  the 
position  in  which  it  was  received— 14  degrees  at  the  tip.  Four  vanes 
were  immersed  in  the  wake  of  the  propeller  and  were  designed  to  be 
used  for  control  of  the  AROD  in  hovering  flight.  Each  of  the  four 
vanes  had  a  ten-inch  span  and  eight-inch  chord.  End-caps 
approximately  two  inches  high  were  placed  on  each  end  of  the  vanes, 
presumably  to  reduce  spanwise  loss  of  lift. 

Self-contained  electrical  power  was  designed  to  be  delivered 
from  an  engine  driven  generator.  For  engine  testing  purposes,  an 
external  power  source  was  also  available  [Ref  16].  The  engine  and 
engine  installation  in  the  test  AROD,  named  "Annie",  were  considered 
representative  of  the  eight  vehicles  obtained  by  NPS. 


2.  Support  and  control  equipment 

Control  of  vane  and  throttle  servos  was  affected  by  pulse 

width  modulated  signals  in  one  of  two  ways:  (1)  via  an  umbilical  line 

connected  to  an  IBM  PC  running  a  "C"  software  code,  or  (2)  via  radio 

(S> 

frequency  control  from  an  off-the-shelf  Futaba  radio  control  i,R/C) 
transmitter  sending  signals  to  a  receiver  mounted  inside  the  vehicle. 
Ignition  was  also  provided  in  two  different  manners;  (1)  external 
power  through  a  28  VDC  power  source,  or  (2)  on  board  engine  driven 
generator.  Several  problems  with  electronics  and  ignition  were 
encountered  during  engine  tests.  These  problems  and  the  control 
set-up  are  described  in  thesis  work  by  Lieutenant  Commander 
Moran  [Ref.  16]. 

3.  Test  stand 

The  test  stand  was  designed  to  measure  thrust  and  rolling 
moment  of  the  AROD.  Test  stand  design  and  manufacture  was 
accomplished  in  work  by  Lieutenant  Kress  [Ref.ll].  Modifications  to 
the  original  design  concerned  improvement  in  accuracy  of  force  and 
moment  by  replacement  of  the  original  "fish  scale"  apparatus  with 
sensitive  mechanical  force  gages.  Three  force  gages  were  obtained. 
Two  gages,  with  a  combined  capacity  of  175  lb,  were  used  in  parallel 
to  measure  force  and  one  50  lb  capacity  gage  was  mounted  to 
measure  rolling  moment  at  the  end  of  of  a  nine-inch  moment  arm. 
for  a  moment  capacity  of  37.5  ft-lbf.  Figure  IV.  1  shows  the  engine 
test  stand,  force  and  moment  measuring  equipment,  and  the  test 
article,  an  AROD.  Force  was  measured  by  adding  the  pair  of  force 
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gages  mounted  on  an  adjustable  plate  which  was  in-line  with  the 
thrust  shaft. 

C.  SCOPE  AND  METHOD  OF  TESTS 

The  tests  were  conducted  on  two  separate  occasions  for  a  total 
engine  run  time  of  approximately  four  hours.  Tests  were  performed 
at  the  engine  run  block-house  located  at  the  NPS  Laboratory.  Tests 
were  conducted  in  accordance  with  the  test  matrix  shown  in 
Appendix  A,  Table  A.I.  Data  included  temperature,  pressure,  thrust, 
engine  speed  in  rpm,  vane  position,  moment  reading,  and--when 
applicable— computer  inputs.  Vane  deflection  was  read  directly  off 
protractors  mounted  at  the  base  of  each  vane  and  calibrated  in  one- 
degree  increments.  Further  discussion  of  vane  numbering  and 
deflection  convention  follows.  Rpm  was  measured  by  an  oscilloscope 
connected  to  the  ignition  source  or  by  a  sensor  mounted  on  the  side 
of  the  R/C  transmitter  which  essentially  detected  light  oscillations 
caused  by  propeller  passage.  Velocity  of  airflow  in  the  area  of  the 
vanes  were  measured  with  a  hand-held  digital  wind-speed  meter. 
Double  hearing  and  eye  protection  were  worn  by  all  test  personnel. 
Data  were  recorded  on  pre-printed  data  sheets  and  by  a  video 
camera. 

D.  RESULTS  AND  DISCUSSION 

1.  General 

The  data  for  engine  tests  completed  are  shown  in  Appendix 
A,  Tables  A.II  and  A.III. 
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2.  Thrust  measurements 

The  first  task  was  to  determine  what  effect  the  mufflers  had 
on  thrust  and  noise  output.  The  engines  were  apparently  operated 
by  the  previous  user  with  mufflers  off  and  a  straight-tube  exhaust 
installed.  In  their  production  configuration,  the  Herbrandson 
mufflers  would  not  fit  inside  the  AROD  duct.  Several  sets  of  brand- 
new  mufflers  were  contained  in  the  spare  parts  obtained  by  NFS. 
FRL  personnel  built  an  exhaust  manifold  adapter  that  allowed 
installation  of  the  mufflers  in  the  constrained  AROD  duct.  Thrust 
reading  at  a  full  throttle  rpm  of  7595  with  the  mufflers  on  was  118 
Ibf.  The  engine  was  shut  down  and  the  mufflers  replaced  by  the 
original  straight  tube  exhaust  system.  Upon  restart  and  re¬ 
establishment  of  full  throttle,  again  at  7595  rpm,  the  thrust  had 
increased  over  four  percent  to  123  Ibf.  Temperature  and  pressure 
for  both  runs--separated  by  less  than  five  minutes--were  the  same. 
Noise  was  substantially  increased  by  removal  of  the  mufflers,  as 
discussed  in  a  later  paragraph,  but  since  initial  simulation  results  had 
showed  a  need  for  the  highest  available  thrust-to-weight  ratio 
(T/W),  the  decision  was  made  to  conduct  subsequent  tests  with  the 
straight-tube  exhausts.  This  test  was  repeated  twice  to  verify  a 
consistent  result.  The  surprisingly  large  difference  in  thrust  at  the 
same  temperature,  rpm  and  blade  angle  may  have  been  due  tc  drag 
on  the  mufflers,  which  were  in  the  wake  of  the  propeller.  The 
frontal  area  and  blockage  effect  of  the  mufflers  was  approximately 
32  square  inches  total. 


Thrust  measurements  were  considered  accurate  to  within 
1/2  Ibf.  Engine  speed  versus  thrust  is  shown  in  Figure  IV. 2.  A 
second  order  polynomial  curve  fit  is  overlaid  and  shows  good 
correlation.  Figure  IV. 3  shows  the  thrust  that  resulted  at  different 
throttle  settings  as  well  as  a  polynomial  curve  fit.  The  throttle 
position  on  this  plot  was  arrived  at  by  taking  the  computer’s  throttle 
position  and  mapping  it  linearly  to  a  zero  to  100  percent  scale.  This 
thrust  versus  throttle  position  curve  fit  was  used  in  the  simulation 
program,  as  discussed  in  Section  VII.  As  shown  in  Figure  IV. 3,  the 
thrust  output  at  idle  was  over  20  Ibf.  This  relatively  high  minimum 
thrust  setting  may  cause  the  vehicle  to  be  more  difficult  to  slow 
down  in  horizontal  flight.  Further  testing  may  be  required  with  the 
idle  adjustment  set  for  a  lower  rpm. 

3.  Vane  control  power  measurements 
a.  General 

The  rolling  moment  imparted  by  deflection  of  vanes 
was  measured  on  the  test  stand.  Moment  readings  tended  to 

"wander"  over  a  range  of  approximately  ±  3  Ibf  force  gage  reading 
which,  through  a  nine-inch  test  stand  moment  arm,  converted  to 
approximately  +  2  ft-lbf  accuracy  for  rolling  moment.  This  range  of 
readings  was  inconsistent  and  at  times  the  readings  wandered  5  ft- 
lbf  from  the  nominal  reading.  The  source  of  this  relatively  large 
variation  in  readings  could  not  be  determined.  Engine  speed 
fluctuations  were  ruled  out  in  that  aural  and  measured  rpm  we 
constant,  as  were  the  thrust  readings.  Wind  in  the  test  area  was  light 
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Figure  IV.2;  AROD  Engine  Test  Results-Thrust  Versus  Rpm 
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Figure  IV. 3:  AROD  Engine  Modelling— Thrust  Versus  Throttle  Position 
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and  variable.  The  control  vane  positions  were  visually  checked  as 
steady.  One  possible  explanation  for  the  variable  moment  readings 
would  be  unsteady  flow  within  the  duct— flow  over  the  engine  and 
its  support  structure  is  downstream  of  the  propeller  but  just 
upstream  of  the  vanes. 

In  keeping  with  MILSPEC  convention  for  roll  control 
deflection  [Ref  17]  positive  vane  deflections— clockwise  when  viewed 
from  the  vane  tip--were  those  which  produced  positive  rolling 
moments— clockwise  viewed  from  the  rear  of  the  AROD.  Zero  vane 
deflection  was  defined  as  the  deflection  at  which  a  vane  was  aligned 
with  the  longitudinal  axis  of  the  AROD.  The  "top"  of  the  AROD  when 
mounted  on  its  side  in  the  test  stand  as  shown  in  Figure  IV.  1  was 
considered  to  be  the  point  where  two  fuel  vents  were  next  to  each 
other.  Vanes  were  then  numbered  according  to  Figure  IV. 4  below. 


AROD  Duct 


Top  Of  AftOD  defined  by 
adjacent  fuel  vents 


VIEWED  FROM  REAR 


Figure  IV.4:AROD  Vane  Numbering  Convention  for  Engine  Tests 


b.  Dual-vane  roll-control  power 

Dual-vane  roll-control  power  was  obtained  by  deflection 
of  vanes  two  and  four,  while  one  and  three  were  set  to  zero 
deflection.  The  rolling  moment  response  to  this  vane  deflection, 
shown  in  Figure  IV.5,  was  fairly  linear  and  showed  no  propensity  to 
drop  off  at  high  deflection  angles.  The  combined  effects  of  engine 
and  propeller  torque  were  overcome,  for  the  two-vane  case,  by 
approximately  18  deg  of  deflection. 

c.  Four-vane  roll-control  power 

Four-vane  roll-control  power  was  obtained  by 
deflection  of  all  vanes  in  the  same  direction  so  as  to  produce  a  rolling 
moment.  The  response  to  this  vane  deflection  with  throttle  set  to  full 
open  (7520  rpm)  is  shown  in  Figure  IV. 6  and,  as  with  dual  vane 
power,  was  fairly  linear  and  showed  no  propensity  to  drop  off  at 
high  deflection  angles.  The  combined  effects  of  engine  and  propeller 
torque  were  overcome,  for  the  four-vane  case,  by  approximately  12 
deg  of  deflection. 

The  change  in  four  vane  roll  control  power  with  changes 
in  throttle  setting  was  investigated  for  two  additional  rpm  settings: 
7240  and  6580.  These  rpm  levels  corresponded  to  approximately  93 
and  75  percent  of  full  throttle  thrust  ,  respectively.  As  shown  in 
Figure  IV. 7,  there  was  little  effective  change  in  roll  control  power  as 
throttle  was  reduced,  even  though  the  thrust  was  decreased  by  25 
percent.  As  discussed  in  Figure  VII,  this  surprising  result  was 
researched  and  it  was  discovered  that  the  AROD  manufacturer. 
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Average  Vene  Deflection  for  Venes  2  &  4,  VPavg,  deg 


Figure  IV.5:  AROD  Engine  Test-Dual  Vane  Roll  Control  Power 
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Sandia  National  Lab,  had  similar  results  during  testing  [Ref  18]. 
d.  Pitch  control  power  modelling 

Data  from  the  engine  tests  were  used  by  other  team 
members  to  model  roll  control  power  [Ref  19].  Additionally,  the  data 
were  used  to  model  the  pitch  control  power  due  to  vane  deflection 
for  use  in  the  3-DOF  simulation,  discussed  in  Figure  VII.  Details  of 
this  transformation  from  roll  control  power  to  pitch  control  power 
follow. 

Pitch  control  power  from  the  vanes  in  the  Archytas 
airframe  will  be  implemented  by  deflection  of  vanes  two  and  four  as 
shown  in  Figure  IV .4.  They  will  both  be  deflected  so  as  to  cause 

pitching  motion.  In  other  words,  to  pitch  the  "nose"  of  the  Archytas 

down  vane  two  will  be  rotated  counter-clockwise  (negative  for  the 
roll  control  convention)  while  vane  four  will  be  rotated  clockwise 
(positive  under  the  roll  control  convention).  Since  this  will  involve 
deflection  of  two  vanes,  the  data  from  the  dual-vane  roll-control 
power  tests  were  used  to  model  pitch  control  power. 

Dual-vane  roll  control  was  assumed  to  occur  from  the 
lift  of  the  two  vanes  at  mid-span,  as  shown  in  Figure  IV. 8,  with  each 

vane  producing  an  equal  lifting  force,  F. 

The  rolling  moment  generated  by  these  pair  of  vanes, 

then,  is: 


L  =  2F 


(L  =  rolling  moment  in  ft  -  Ibf) 


(IV.A) 


34 


Equating  this  rolling  moment  model  to  the  curve  fit  of 
dual  vane  rolling  moment  test  data  from  Figure  IV. 5: 

L  =  \.5F  =  -7.86  +  (0.400)(VP) 

(Recall  that  VP  is  vane  deflection)  (IV.B) 

Solving  equation  IV.B  for  F: 

F  =  -4.24  +  (0.267)(VP)  (IV.C) 

Pitching  moment  will  then  be  a  result  of  twice  this  same 
lifting  force  applied  over  a  moment  arm  from  the  vane  quarter  chord 
to  the  longitudinal  center  of  gravity,  as  depicted  in  Figure  IV.9: 

Which  gives  a  pitching  moment  due  to  vane  deflection 
of: 


=  -2FixcgJ 

or,  substituting  for  F  from  equation  IV.C: 

=  -2{-4.24  +  (0.267)(VP)}(j:cgJ  (jV.D) 

The  number  -4.24  in  ft-lbf  from  equation  IV.C 
represents  the  rolling  moment  zero-vane-deflection  position,  as 
shown  on  Figure  IV.5.  Because  the  propwash  characteristics  will  be 
different  when  vanes  are  deflected  for  pitching  moment,  and  to 
simplify  user  input  to  the  simulation  program— i.e.  user  will  not  be 
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AROD  Engine  Tests—Change  in  Roll  Control  Power  With  Throttle 


iew  looking  down  AROD  longitudinal  axis 
with  vanes  1  and  3  removed  for  clarity. 
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required  to  know  that  a  certain  number  of  degrees  of  vane  deflection 

are  required  to  counter  the  swirl  of  propwash  in  the  duct--this 

quantity  was  eliminated  from  the  analysis.  Note  that  equation  IV.D 

accounts  for  the  sign  change  required  to  reflect  the  difference  in  sign 

convention  between  the  roll  axiS"Where  positive  deflections  produce 

positive  moments--and  the  pitch  axis—where  positive  deflections 

produce  negative  moments  [Ref  17].  Note  also  the  change  in  variable 
from  VP  to  6^,  used  to  define  vane  deflections  in  a  pitching  sense  for 

the  simulation  program.  Making  these  two  changes  to  equation  IV.D 
as  well  as  inserting  the  final  value  for  xcgy  of  1.5  ft,  which  was 

arrived  after  considering  issues  discussed  in  Figure  VI,  the  final 
model  for  pitching  moment  due  to  vane  deflection,  used  in  Figure 
VII,  was: 

M_=-(0.80)(<5_) 

where  5,^  h  vane  deflection  in  degrees  (IV  E) 

4.  Effect  on  thrust  from  blockage  in  AROD  duct 

The  changes  in  thrust  due  changes  in  vane  deflection, 
termed  "blockage  effect",  were  documented  by  noting  changes  in 
thrust  due  to  vane  deflection  at  three  throttle  settings.  This  analysis 
was  performed  primarily  to  predict  the  change  in  thrust  with  vane 
deflection  for  improved  computer  modeling  of  thrust,  particularly  in 


a  hover,  where  engine  thrust  levels  are  critical.  A  secondary  goal  of 
these  tests  was  to  attempt  to  characterize  the  swirl  of  propwash  in 
the  duct.  The  effect  on  thrust  due  to  vane  deflection  is  shown  in 
Figure  IV.  10.  Regardless  of  throttle  setting,  the  maximum  thrust 
was  achieved  at  about  -5  deg  of  vane  deflection.  The  differences  in 
indicated  thrust  were  likely  a  result  of  changes  in  the  local  angle  of 
attack,  and  hence  vane  drag,  due  to  swirl  in  the  duct.  Apparently  the 
local  velocity  at  the  vanes  was  approximately  in-line  with  the  vane’s 
symmetrical  airfoil  section  when  the  vanes  were  at  -5  deg.  At  vane 
deflections  off  this  peak-thrust-value  the  thrust  fell  off  steadily, 
decreasing  approximately  ten  percent  at  full  throw  in  each  direction. 
Incorporation  of  these  results  will  increase  the  fidelity  of  future 
simulations. 

5.  Impact  of  propwash  on  vane  position 

Early  in  engine  testing  it  was  observed  that,  when  using 
computer  control  of  vane  position,  the  vanes  would  be  at  a  different 
deflection  when  the  engine  was  running  than  they  had  been  during 
position  calibration  with  the  engine  shut  down.  Furthermore,  as 
throttle--and  presumably  local  velocity  at  the  vanes--was  increased 
the  vane  position  would  change  with  no  change  in  the  commanded 
position.  This  effect,  which  came  to  be  known  as  "blowback  effect", 
was  documented  during  testing  completed  on  April  16,  1993.  Vane 
position  was  calibrated  with  the  engine  off  through  a  program 
written  in  "C"  and  run  on  an  IBM  PC  386.  Reference  16  describes  the 
computer  control  set-up.  Figure  IV.  11  shows  the  change  in  vane 
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deflection  with  the  engine  running.  A  general  trend  in  the  errors  is 
apparent  in  Figure  IV.ll,  with  errors  approximately  linear  from 
minus  IS  to  plus  ten  degrees  and  then  dropping  off  at  higher 
commanded  vane  positions. 

Engine  tests  performed  on  April  27,  1993  were  conducted 
using  a  Futaba  R/C  set  to  command  vane  and  throttle  positions. 
When  using  radio  control,  this  blowback  error  did  not  occur. 

6.  Flow  speed  at  vane  locations 

The  speed  of  airflow  at  two  locations  near  the  vanes  was 
measured  with  a  hand-held  digital  wind  meter  for  future  use  in  duct 
aerodynamic  modeling  and/or  control  power  prediction.  Readings 
fluctuated  approximately  +  2  ft/sec.  These  fluctuations,  which  may 
have  been  due  to  actual  fluctuations  in  flow  speed  or  to  wind  meter 
inaccuracies,  were  considered  insignificant  to  the  goal  of  the  analysis. 
Flow  speed  was  gathered  over  a  range  of  speeds  at  two  different 
locations  as  shown  in  Figure  IV.  12.  Rotor  theory  predicts  an  induced 
velocity  at  the  disk-plane  of  a  rotor  according  to  the  following 
formula  [Ref  20]: 


where: 


(IV.F) 


s  Induced  velocity  in  ft/sec 
T  s  Thrust  in  Ibf 
p  =  density  in  slugs/ft''  3 
A  =  Rotor  area  in  ft^2 
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Figure  IV.  11 

AROD  Engine  Test—Blowback  Effect 

Survey  results  are  shown  in  Figure  IV. 13,  with  a  second 
order  polynomial  curve  fit  displayed  for  survey  position  #1.  Also 
shown  in  Figure  IV.  13  is  the  rotor  theory  prediction  for  a  12-inch 
radius  rotor  at  the  test  day  density  of  0.002382  slugs/ft'^3.  The  data 
followed  the  slope  predicted  by  theory  but  with  an  approximate 
average  increase  in  velocity  of  12  ft/sec.  This  velocity  above 
prediction  was  probably  due  to  a  combination  of  three  factors: 

•  The  duct  was  not  accounted  for  by  rotor  theory  but 
would  have  the  effect  of  containing  and  speeding  the  flow. 

•  The  value  for  thrust,  T,  used  in  equation  (IV.F  )  was  not 


actual  thrust  but  rather  engine  thrust  minus  drag  from  the  engine 
and  its  supports  well  as  duct  wall  and  vanes. 

•  The  hub  of  the  propeller,  at  approximately  15  percent  of 
the  prof)cller  diameter,  was  larger  than  that  on  most  rotor  systems. 


Measurement  position  #2: 

At  vane  mid-chord  and  mid-span, 
approximately  one  inch  from 
vane  upper  surface. 


Measurement  position  #1 : 

At  vane  mid-chord,  approximately 
equidistant  between  the  two  vanes 
and  between  center  of  duct  and 
duct  shdi. 


VIEWED  FROM  REAR 

Figure  IV.12:AROD  Engine  Tests-Locations  of  Flow  Velocity  Survey 


7.  Miscellaneous  engine  tests 
a.  Noise 

The  sound  level  was  measured  for  the  purpose  of 
determining  hearing  protection  required  for  test  personnel. 
Additionally,  the  noise  level  figured  into  the  decision  on  whether  to 
operate  the  AROD  with  or  without  its  mufflers,  as  discussed  in  a 
previous  paragraph.  Double  hearing  protection  must  be  worn  in 
areas  with  sound  levels  higher  than  90  dB.  With  mufflers  not 
installed,  and  when  within  about  20  ft  of  the  AROD  at  full  power,  the 
noise  level  exceeded  the  noise  meter  limit  of  120  dB.  Installation  of 


mufflers  reduced  the  sound  levels  by  about  6  dB.  For  the  mufflers- 
off  case,  a  separation  distance  of  greater  than  100  ft  was  required  to 
bring  sound  levels  below  the  double-hearing-protection  limit.  These 
measurements  were  taken  with  the  AROD  positioned  at  the  FRL  test 
cell,  which  is  surrounded  by  earth  berms.  Further  testing  will  be 
required  to  determine  sound  levels  for  flat  areas,  such  as  that 
expected  for  flight  test. 

b.  Fuel  flow 

Fue*  usage  was  recorded  for  three  separate  engine  runs 
at  sustained  high  power,  at  or  above  80  percent  of  maximum  rpm. 
The  average  fuel  consumption  for  these  high  power  conditions, 
approximately  representative  of  hovering  flight,  was  3.6  gallons  per 


44 


hour.  Qualitatively,  fuel  usage  at  reduced  power  was  significantly 
less.  No  attempt  was  made  to  lean  the  engine  to  reduce  fuel  flow. 

c.  Acceleration/deceleration  characteristics 

The  acceleration  and  deceleration  characteristics  were 
evaluated  by  attempting  to  record  rpm  versus  time  during  rapid 
throttle  movements.  The  response  of  the  engine  was  too  fast  to  allow 
observation/recording  of  the  rpm  on  an  oscilloscope.  Qualitatively, 
by  listening  to  engine  response,  the  spool  up  and  down 
characteristics  were  first  order,  with  a  time  constant  less  than  one 
third  of  a  second. 

d.  Reliability 

During  the  three  data  periods  and  several  other  practice 
and  engine  familiarization  sessions,  no  engine  related  failures 
occurred.  However,  several  failures  occurred  with  the  electronics 
which  caused  the  engine  to  shut  down.  These  failures  generally 
occurred  as  a  result  of  the  AROD’s  high  vibration  levels.  The  failures 
included  broken  solder  joints  and  wiring,  and  damaged  circuit  board 
components.  Electronics  failures  and  related  issues  are  discussed  in 
the  thesis  by  LCDR  Moran  [Ref.  16]. 


V.  WING  SPAR  DESIGN,  FABRICATION  AND  TEST 


A.  GENERAL 

The  purpose  of  designing  a  spar  was  to  connect  the  existing 
Aquila  wing,  and  the  wing  root  support  structure,  which  came  to  be 
known  as  the  "shoulder",  to  the  AROD  duct.  Previous  work  by  Kress 
had  specified  a  general  design  of  two  composite  spars  in  the  shape  of 
rings  that  would  surround  the  AROD  duct  [Ref.  11].  Initial  work  with 
cardboard  mock-ups  showed  this  design  to  be  impractical  for 
installation  on  the  AROD  duct,  due  to  the  AROD’s  longitudinal  struts 
and  landing  gear,  and  so  the  spar  was  redesigned.  The  new  design 
consisted  of  a  two  spar  configuration  with  a  shape  similar  to  the 
previous  design  but  spanning  only  90  deg  of  the  AROD  duct,  instead 
of  the  original  design's  360  deg. 

Several  structural  tests  were  conducted  to  determine  component 
strength.  A  final,  full  configuration  test  was  then  conducted  to 
determine  the  strength  of  a  production-like  installation. 

B.  WING  SPAR  DESIGN  AND  FABRICATION 

1 .  Design  load 

The  spar  was  designed  to  withstand  8  g’s  ultimate  load 
factor  and  5.33  g's  limit  load  factor  for  a  100  Ibf  vehicle.  This 
converted  to  400  Ibf  of  lift  load  per  wing  for  ultimate  strength  and 
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267  Ibf  for  design  limit.  For  purposes  of  spar  design  and  test,  this  lift 
load  was  assumed  to  occur  as  a  point  force  at  the  quarter  chord  of 
the  mean  aerodynamic  chord,  as  shown  in  Figure  V.l.  The  moment 
at  the  wing  root,  then,  was  this  lift  load  times  the  distance  from  the 
mean  aerodynamic  chord  (MAC)  to  the  wing  root,  30.5  inches.  This 
resulted  in  a  predicted  moment  of  12,200  and  8,128  in-lbf  for 
ultimate  and  design  limits,  respectively.  The  load  for  each  of  the  two 
spars  was  then  assumed  to  be  carried  as  a  function  of  its  longitudinal 
distance  from  the  point  lift  force.  Ultimate  and  design  limits  for  each 
spar  are  detailed  in  Table  V.l.  A  sample  calculation  for  the  aft  spar 
is  presented  below: 


AFT  SPAR 


design  moment 


_  (Wt)(Nz)(y) 


liESL 


x" 


(#  WINGS)  ^  , 

=  (100  lbO(5.33)(30.5  in)n6-5 

2  I  16 

=  5588  in  -  Ibf 


where : 


Wt  =  design  weight 
Nz  =  load  factor 

y  =  y  coordinate  distance  to  MAC 
X.p„  =  longitudinal  distance  between  spars 

X  =  longitudinal  distance  from  design  spar 
quaner  chord  of  MAC 


Table  V.I 

ARCHYTAS  WING  SPAR  STRUCTURAL  DESIGN  LIMITS 


SPAR 

ULTIMATE  MOMENT 

DESIGN  MOMENT 

(8  g's)  in-lbf 

(5.33  g's)  in-lbf 

Forward 

3812 

2540 

Aft 

8388 

5588 

TOTAL 

12200 

8128 

2.  Forward  spar  design  and  attachment 

Figure  V.I  depicts  the  two-spar  system  conceived  for  use  on 
the  Archytas  flight  article  [Ref.  21].  The  basic  structure  consisted  of 
a  combination  of  graphite  and  fiberglass  in  an  epoxy  matrix 
surrounding  a  foam  core.  The  two  spars  were  connected  to  the  AROD 
in  different  ways.  The  forward  spar  was  affixed  to  the  duct’s 
exterior  by  three  methods  (refer  to  V.I): 

•  The  spar  was  cut  to  fit  the  duct  surface  as  closely  as 
possible  and  was  attached  with  a  special  structural  adhesive  to  the 
outside  of  the  duct.  Adhesive  thickness  varied  depending  on  the 
quality  of  the  fit  between  the  spar  and  the  duct— approximately 
1/32  to  1/8  inch.  The  load  path  for  this  connection  method— termed 
"peel  strength” --was  from  the  surface  of  one  composite  component 
(the  AROD  and  the  forward  spar),  through  the  structural  adhesive 
and  into  the  other  component. 
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•  Grip  plates  were  manufactured  from  0.053 -inch  thick 
6061-00  aluminum  and  configured  with  0.25  in  nut  plates.  These 
grips  were  attached  to  the  forward  spar  surface  with  the  structural 
adhesive,  with  a  targeted  adhesive  thickness  of  5.0  mm.  As  specified 
in  industrial  references,  to  ensure  a  good  bond  between  the 
aluminum  grip  plates  and  the  composite/epoxy  spar,  the  aluminum 
was  bead-blasted  and  then  cleansed  by  a  commercially  available 
trichloroethane  product  ("Brake  Clean")  [Ref  22  and  23].  Research 
was  conducted  to  determine  the  feasibiliiy  of  first  etching  the 
aluminum  plates  with  chromic  acid  and  then  phosphoric  anodizing 
them,  in  an  attempt  to  optimize  the  bonding  surface  [Ref  22]. 
Discussions  with  personnel  in  the  NFS  Physics  department  as  well  as 
commercial  anodizing  contractors  revealed  that  this  process  was 
environmentally  risky  and  required  special  equipment  and  chemicals 
not  readily  available  [Ref  24  and  25].  A  decision  was  made  to  utilize 
a  mechanical  abrasion  method  (bead  blasting).  Additionally,  the  grip 
plates  had  number  25  holes  drilled  in  them  so  that  the  structural 
adhesive  could  flow  through  for  improved  bonding  [Ref  26].  The 
grips  were  attached  to  the  AROD  duct  by  0.25  in  bolts  passing 
through  their  own  load-bearing  plates,  which  were  adhered  to  an 
AROD  strut  with  the  structural  adhesive,  then  through  the  strut  and 
into  the  grip's  nut  plate,  as  shown  in  Appendix  B,  Figure  B.l 

•  The  last  connection  made  between  the  forward  spar  and 
the  duct  was  through  a  aluminum  grip  plate  of  the  same  material  as 
above  but  oriented  radially.  In  other  words,  the  base  of  the  grip 


plate  (with  an  integral  0.25  in  nut  plate)  rested  on  the  duct's  exterior 
surface  and  this  plate  was  attached  to  the  forward  spar  at  a  90-deg 
angle  to  the  previously  described  grips.  A  single  0.25  in  bolt  was 
passed  through  a  plate  mounted  on  the  inside  of  the  AROD  duct,  as 
shown  in  Appendix  B,  Figure  B.2,  and  into  the  nut  plate  on  the  grip. 
The  bolt  was  positioned  below  the  propeller  to  avoid  Foreign  Object 
Damage  (FOD)  to  engine  and  propeller  components.  The  purpose  of 
this  last  connection  was  to  improve  the  spar's  "peel  strength". 

3 .  Aft  spar  design  and  attachment 

The  aft  spar  was  similar  in  construction  to  the  forward  spar 
except  that,  owing  to  reduced  area  available  for  mounting  in  the 
wing  root,  it  was  2.0  in  thick,  1.0  in  thicker  than  the  forward  spar. 
The  aft  spar  had  aluminum  grips  mounted  in  the  same  manner  as  the 
forward  spar.  These  grips  were  bolted  to  the  AROD's  aluminum 
landing  gear  legs  and  constituted  the  entire  mounting  strength  for 
the  aft  spar.  Appendix  B,  Figure  B.3  shows  grips  installed  on  both 
spars. 

4.  Integration  of  NFS  spars  with  Aquila  wing  spars 
The  forward  and  aft  spars  manufactured  at  the  FRL  were 

mated  to  the  Aquila  wing's  root  by  a  combination  of  bolts,  plates  and 
structural  adhesive  as  shown  in  Figure  V.2.  The  Aquila  wing  was 
then  attached  to  its  wing  root  by  its  original  system  of  bolts  and  nut- 
plates. 
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Figure  V.l;  Archytas  Wing  Spar  Design 


adhesive  spar 

VIEW  LOOKING  FROM  WING-TIP  DOWN  WING  LINE 

Figure  V.2 

Archytas  FRL/Aquila  Wing  Spar  Integration 

C.  WING  SPAR  EXPERIMENTS 
1 .  Initial  component  testing 

Forward  and  aft  spars  were  first  fabricated  as  detailed  in 
previous  work  [Ref  11],  at  which  point  a  structural  test  was 
performed  on  these  test  specimens  to  compare  actual  strength  with 
those  predicted  in  the  finite  element  analysis  of  the  previous  thesis 
work.  The  structural  test  at  this  point  of  development  was  on  the 
individual  spar  component.  Figure  V.3  depicts  the  set-up  for  the 


structural  test  on  an  individual  spar.  The  upper  part  of  the  figure, 
part  a,  shows  the  spar  in  a  notional  installation.  The  lower  part  of 
Figure  V.3  depicts  the  experiment  to  test  spar  strength.  The  test 
spar  was  mounted  in  a  supporting  fixture  designed  to  test  the 
narrow  area  of  the  spar  and  the  bonding  of  the  structural  adhesive  to 
the  aluminum  plate,  as  depicted  in  Figure  V.2.  The  test  spar  was 
mounted  upside  down  so  that  the  weight  placed  on  the  beam  would 
simulate  a  lift  load.  Appendix  B  ,  Figure  B.4  is  a  photograph  of  the 
test  set-up. 

The  215  Ibf  weight,  a  stainless  steel  bar,  was  set  in  place  by 
a  hand  forklift  at  known  distances  marked  on  the  steel  beam.  The 
42  Ibf  steel  tube  was  then  placed  on  the  steel  beam  and 
progressively  moved  farther  out  the  beam.  Once  the  steel  tube  was 
at  the  beam’s  end  it  was  removed,  and  the  215  Ibf  weight  was  then 
moved  out  enough  to  produce  the  same  total  moment  that  had 
existed  when  the  bar  was  in  its  former  location  with  the  tube  at 

beam's  end.  The  tube  was  then  replaced  on  the  beam,  near  the  test 

spar,  and  then  again  progressively  moved  out  the  beam.  The  weight 
of  the  steel  beam  was  also  accounted  for  and  was  assumed  to  act  as  a 

point  force  at  half  its  length.  This  process  continued  until  either  a 

failure  occurred  or  a  desired  moment  value  was  reached. 


NPS 
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Aquita  Wing  Spar 


To  mean  aero 
chord,  7  =  30.5  in 


AROD  Struts 


!/ 


AROD  DUCT 


Figure  V.3 

Archytas  Wing  Spar  Installation  and  Test 
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With  the  combination  of  tube,  beam  and  bar  weights 
simulating  an  aerodynamic  lift  load,  the  moment  supported  by  the 
test  spar  was  calculated  as  follows  (refer  to  Figure  V.3): 

^  Moments  -  0 

215(4)  +  42(/,)  +  W'(/)-Af^=0 

or 

Mspar  =  215(4)  +  42(4 )  + 1^(0 

The  loading  was  placed  slowly,  i.e.  the  weights  were  not 

dropped  on  the  beam,  and  were  left  in  position  for  at  least  two 

minutes  without  any  failure  of  the  spar.  Failure  of  the  spar  was 

quantified  by  listening  for  cracking  noises  in  the  vicinity  of  the  spar. 
Typically,  when  the  load  was  applied  there  would  be  a  few  small 
cracking  noises  as,  presumably,  a  few  of  the  fibers  and/or  portions  of 
the  matrix  broke  and  "settled  in".  These  "micro-failures"  were  minor 
and  short  lived  and,  if  a  load  was  relaxed  and  re-applied,  generally 
did  not  reoccur.  Failure  was  generally  marked  by  a  slow,  continuous 
series  of  cracking  noises  which  did  not  cease,  although  the  time 

interval  between  cracks  was  often  on  the  order  of  a  minute,  finally 
culminating  in  an  actual  failure  of  the  whole  structure.  The  time 
period  from  application  of  the  load  to  the  actual  failure  took  as  long 
as  ten  minutes.  All  of  this  is  to  say  that  the  structural  testing  was 
considered  to  be  fairly  conservative,  in  that  the  moment  was  not 
rapidly  increased  in  an  effort  to  increase  the  apparent  strength  of 
the  component.  Table  V.II  shows  the  loading  sequence  for  the  initial 


test  of  the  forward  spar.  Recall  from  Table  V.I  that  the  goal  was  to 
reach  an  ultimate  moment  of  3812  in-lbf  for  this  component. 


Table  V.II 

ARCHYTAS  FORWARD  SPAR  TEST  (BEAM=16.5  Ibf,  L=74  in) 


Lb  (in)(l) 

Lg  (in)(  ^ ) 

Notes  in-lbf— Goal=3812) 

6 

Not  applied 

2  or  3  small  cracking  noises 

(M=2508^ 

6 

1  2 

No  cracking  (M=3012) 

6 

1  8 

10-15  small  cracking  noises, 

appeared  to  come  from  plate 

area  (M=3264) 

6 

24 

Small  cracking  noises  continuous, 

every  10-15  seconds.  Failure 

aficr  approximately  10  minutes 

(M=3516) 

Notes:  (1)  Refer  to  Figure  V.3  for  dimensions  and  experiment 

method 


The  moment  at  failure  was  92  percent  of  the  ultimate, 
corresponding  to  a  load  level  of  7.4  g’s.  Post-experiment  analysis  of 
the  test  spar  revealed  that  the  failure  may  have  started  at  a  lower 
loading  level  as  the  laminates  separated,  but  that  final  failure  was 
partially  masked  by  the  test  fixture  constraining  a  portion  of  the 


56 


spar.  The  failure  occurred  at  the  juncture  of  the  shear  web  and  spar 
cap  laminates  [Ref  11,  pp  31]  adjacent  to  the  plate  that  was  designed 
to  transmit  the  loads  from  the  Aquila's  wing  spar.  Figure  V.4  depicts 
the  location  of  this  failure. 


TOP  VIEW 


2024>T6  aluminum  plate 
attached  with  structural 
adhesive 


Tszs 


Shear 
web 
laminate 


Area  depicted 
below 


FAILURE-twistlng 
apart  of  structure 


Foam  core 


SIDE  yiEW 


Spar  cap  laminate 


Figure  V.4 

Archytas  First  Test  Spar  Failure  Description 


57 


2.  Redesign  of  basic  spar  structure 

As  a  result  of  the  failure  witnessed  in  the  first  structural 
test,  a  modification  was  made  to  the  basic  spar  design  in  an  effort  to 
increase  its  strength.  Specifically,  the  load  path  between  the  shear 
and  spar  cap  laminates,  which  heretofore  had  consisted  solely  of  the 
relatively  weak  matrix,  was  strengthened,  as  depicted  in  Figure  V.5. 
This  modification  was  made  to  the  tip  area  of  an  existing  aft  spar. 
These  modifications  changed  the  weight  of  this  particular  test  aft 
spar  from  0.77  to  1.51  Ibf,  an  increase  of  49  percent,  but  it  was  felt 
that  the  weight  of  a  spar  built  with  these  modifications  incorporated 
during  original  construction  would  not  gain  as  much  weight.  In  any 
case,  the  original  design  was  clearly  unsatisfactory  and  in  need  of 


modification. 


Modifications  made 
after  failure.  1  -Z  layers 
of  uni-directionai 
graphite  fibers  wrapped 
between  shear  and  spar 


Figuie  V.5:  Modifications  Made  to  Archytas  Spar 
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3 .  Testing  of  modified  spar 

The  modified  spar  was  tested  in  the  same  way  as  the 
original  test  spar,  except  that  the  beam  on  which  the  weights  were 
placed  was  replaced  with  a  shorter  version.  Table  V.III  describes 
the  loading  pattern.  The  goal  for  this  test  aft  spar  was  to  reach  an 
ultimate  moment  of  8388  in-lbf. 


Table  V.III 

ARCHYTAS  MODIFIED  AFT  SPAR  TEST  (BEAM=12  Ibf,  L=37  in) 


Lb  (in)(l) 

Lg  (in)(l) 

Notes  (Mgpaj.  in-lbf--Goal=8388) 

1  2 

36 

Nothing  noted  (M=4536) 

1  2 

48 

Nothing  noted 

12 

54 

Nothing  noted  (M=5292) 

1  8 

Nothing  noted  (M=4314) 

1  8 

30 

Nothing  noted  (M=5574) 

1  8 

36 

Nothing  noted 

1  8 

42 

Small  cracking  sounds  near  base 
of  test  fixture  (M=6078) 

1  8 

48 

tt 

1  8 

54 

Nothing  noted 

1  8 

60 

Nothing  noted  (M=6834) 

24 

Nothing  noted  (M=5604) 

24 

36 

Small  cracking  sounds  near  base 
of  test  fixture  (M=6078) 

24 

42 

As  above  (M=7368) 

24 

48 

As  above  (M=7620) 

24 

54 

Failure  of  spar  at  attach  point  to 
fixture  (M=7872--94%  of  goal) 

Notes;  (1)  Refer  to  Figure  V.3  for  dimensions  and  experiment 
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method 


During  the  testing  of  this  modified  spar,  cracking  sounds 
were  heard  and  appeared  to  be  coming  from  the  base  of  the  test 
specimen,  where  it  was  held  by  the  test  fixture.  As  shown  in  Figure 
V.6,  post-test  analysis  of  the  spar  showed  that  the  modified  areas 
were  without  permanent  deformation  and  uninjured.  Appendix  B, 
Figure  B.5  is  a  photograph  of  this  test  spar  showing  the  failure.  The 
failure  occurred  adjacent  to  the  area  that  was  modified.  This  area 
had  not  been  modified  because  it  had  been  assumed  that  the  test 
fixture  would  provide  the  necessary  strength  to  support  the  spar, 
obviating  the  need  to  modify  the  entire  test  spar.  As  mentioned,  the 
modified  area  was  unaffected  by  the  93  percent  of  ultimate  moment, 
correlating  to  7.5  g's.  While  the  results  of  this  second  test  appeared 
to  shadow  the  results  of  the  first  test  it  was  obvious  to  test  personnel 
that  the  modified  spar  was  significantly  improved  because  it  did  not 
bend  or  deform  as  had  the  previous  design.  Because  of  this 
confidence  in  the  modified  design  a  decision  was  made  to 
manufacture  a  series  of  new  spars  using  the  modified  technique  and 
to  proceed  to  the  next  stage  of  testing. 


Failure  point 
(adjacent  to 
modified  area) 


Modified  area  (shown 
larger  than  actual  size 
for  clarity)  No  permanent 
deformation  noted. 


Figure  V.6:  Archytas  Second  Spar  Test  Failure  Description 
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4.  Comprehensive  structural  test 

After  basic  design  problems  were  identified  and  corrected 
during  the  component  testing  phase,  a  test  was  conducted  on  an 
entire  duct,  wing  spar,  and  wing  assembly.  An  AROD  duct  which  had 
been  previously  partially  damaged  was  utilized  as  a  Structural  Test 
Vehicle  (STV).  The  damage  on  the  STV  was  limited  to  areas  not 
affecting  the  outcome  of  structural  testing.  As  depicted  in  Figure  V.7, 
the  STV  was  a  standard  AROD  duct  with  a  prototype  right-hand  wing 
assembly  mounted  to  it.  Components  included  the  NFS  manufactured 
wing  spars  and  their  attaching  hardware  as  described  earlier  in  this 
chapter,  the  canard  boom  with  support  hole  drilled  in  the  forward 
spar  for  it  and  boom  base  mounted  to  the  aft  spar  (see  Section  VI). 
and  the  Aquila  wing  root  bolted  to  the  NFS  spars.  The  STV  was 
assembled  using  the  techniques  and  materials  anticipated  for  use  in 
the  flight  article.  The  beam  used  for  load  application  was  mounted 
so  as  to  simulate  the  assumed  point  lift  force,  as  depicted  in  Figure 
V.l;  in  other  words,  it  was  placed  five  inches  in  front  of  the  aft  spar 
and  eleven  inches  aft  of  the  forward  spar.  Similar  to  the  structural 
tests  conducted  on  individual  components,  the  STV  was  loaded  by  a 
combination  of  large  (215  Ibf)  and  small  (42  Ibf)  weights.  Because 
the  wing  root  ("shoulder")  assembly  used  on  this  full  test  was 
intended  for  use  on  the  flight  article,  the  decision  was  made  to  limit 
the  load  applied  to  the  design  limit  of  5.33  g’s  or  8128  in-lbf.  The 
loading  sequence  and  comments  thereof  are  contained  in  Table  V.IV. 
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The  wing  shoulder  was  removed  and  inspected.  As  shown 
in  Figure  V.8,  the  wing  shoulder  had  begun  ripping  at  a  location 
adjacent  to  the  wing  spar.  No  failures  were  noted  in  the  FRL 
manufactured  parts.  A  decision  was  made  to  suspend  further  testing 
of  the  wing  shoulder  until  it  could  be  modified.  To  determine  the 
strength  of  the  FRL  manufactured  spars  and  their  attachment  with 
the  AROD  duct,  a  test  fixture  was  constructed  to  take  the  place  of  the 
wing  shoulder  and  allow  loading  of  just  these  components.  This 
experimental  set-up  is  shown  in  Figure  V.9.  Table  V.V  describes  the 
results  of  this  test.  The  decision  was  made  to  limit  the  load  to  design 
limit  (5.33  g)  in  order  to  preserve  the  capability  to  do  testing  with  a 
modified  shoulder. 


Table  V.V 

ARCHYTAS  STRUCTURAL  TEST  OF  FRL  SPAR  INSTALLATION 

(BEAM=12  Ibf,  L=43  in) 


Lb  (in) 

Lj  (in) 

Notes  (Mjpaj.  in-lbf— Goal=8128) 

1  2 

0 

Nothing  noted  (M=3096) 

1  2 

18-72 

Nothing  noted  (final  M=6120) 

24 

30  to  54 

Nothing  noted 

24 

60 

GOAL  REACHED,  test  suspended 

(M=81 96-101%  of  goal) 

Figure  V.9:  Structural  test  of  wing  spars  without  wing  shoulder 
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During  application  of  the  load,  there  was  approximately  five 
degrees  of  bending  in  what  would  be  the  wing.  When  the  load  was 
removed,  the  structure  returned  to  its  original  position.  The  spars 
and  attachments  were  inspected  after  the  test  and  no  cracking, 
permanent  deformation  or  loss  of  strength  were  noted. 
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VI.  CANARD  DESIGN,  FABRICATION  AND  TEST 


A.  GENERAL 

After  formulating  the  initial  concept  of  combining  the  wings 
from  an  Aquila  with  the  AROD  ducted-fan  UAV,  the  next  issue  was 
how  longitudinal  control  and  stability  would  be  provided.  Use  of  a 
conventional,  tail-aft  design  was  discarded  for  several  reasons, 
including  too  high  a  center  of  gravity  for  takeoff  and  landing  and  the 
structural  complexity  of  building  a  combination  tail  and  landing  gear. 
The  other  alternative  was  the  use  of  a  canard.  The  effort  of  the 
present  thesis,  then,  was  to  design  a  canard  and  design,  fabricate  and 
test  a  support  structure  for  the  canard,  that  would  have  the  following 
characteristics: 

•  Good  static  longitudinal  stability  with  a  15  percent  static 
margin,  and  the  aircraft  trimmable  throughout  expected  flight  range 

•  Good  longitudinal  control  power  for  forward  flight 

•  Location  of  the  canard  at  least  one  chord  length  in  front  of  the 
nearest  AROD  structure  to  minimir'  aerodynamic  interference 

•  Lightweight  construction  capable  of  withstanding  expected 
flight  loads 

The  process  that  resulted  in  a  final  canard  size  and  location  was 
iterative  in  nature.  For  example,  static  margin  was  easily  increased 
by  moving  the  canard  forward,  but  because  stronger  (heavier)  booms 
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would  be  required,  this  added  weight  to  a  weight-critical  design.  The 
iterative  nature  of  the  canard  and  boom  design  process  was  as 
follows: 

(1)  The  wing  was  modeled  using  vortex  lattice  techniques  to 
determine  what  moments  the  canard  would  be  required  to  counter. 

(2)  Knowing  the  wing  characteristics,  the  canard  could  then  be 
positioned  and  sized  to  provide  a  certain  static  margin.  This  was 
accomplished  largely  through  the  creation  of  a  MATLAB  program 
"xplot",  which  displayed  the  aircraft's  neutral  point  and  center  of 
gravity  for  a  range  of  canard  sizes  and  locations  [Ref.  27].  Note  that 
at  this  point  the  canard  and  boom  weight  were  assigned  a  nominal 
value,  per  linear  foot  for  the  boom  and  per  square  foot  for  the 
canard. 

(3)  Once  the  canard  location  and  size  were  selected,  further 
analysis  of  the  static  longitudinal  trimmability  was  carried  out 
through  creation  of  another  MATLAB  program,  named  "incidence", 
whose  output  of  pitching  moment  versus  angle  of  attack  graphically 
depicted  the  incidence  angles  required  to  fly  over  a  range  of  speeds, 
or  in  other  words,  the  stability  and  trimmability  of  the  chosen 
design. 

(4)  With  the  canard  configuration  chosen,  the  boom  was 
optimized  by  selecting  the  lightest  tubing  size  capable  of 
withstanding  the  expected  loads. 

(5)  Components  were  then  manufactured  and  tested  to  assure 
that  they  could  indeed  withstand  the  expected  loads. 
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B.  WING  MODELING 

The  airfoil  of  the  Aquila  wing  was  a  modified  NACA  23015 
section.  Modification  to  the  airfoil  consisted  of  a  slight  reflex  at  the 
trailing  edge,  affecting  the  last  ten  percent  of  the  upper  surface  and 
30  percent  of  the  lower  surface  [Ref.  28].  Table  VI.l  shows  the 
Aquila  airfoil  coordinates.  Figure  VI.l  plots  these  coordinates  as  well 
as  the  mean  camber  line,  and  shows  polynomial  curve  fits  for  two 
segments  of  the  mean  camber  line. 

The  Aquila  wing  also  incorporated  a  linear  geometric  twist  of 
three  deg,  with  the  tip  at  a  lower  angle  of  incidence  than  the  root 
[Ref  28].  A  planform  view  of  the  basic  wing  geometry  used  for 
analysis  is  shown  in  Figure  VI, 2.  The  Aquila  wing  sweep  angle  of  28 
deg  was  retained  for  the  Archytas  design  for  several  reasons.  The 
primary  reason  was  the  relative  structural  simplicity  and  associated 
weight  savings  of  mating  the  Aquila  wing  spar  in-line  with  the  spar 
manufactured  by  FRL  personnel.  Additionally,  lower  sweep  angles 
had  the  adverse  effect  of  reducing  the  static  margin  for  a  given 
center  of  gravity,  causing  longitudinal  stability  and  control  problems. 
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Table  Vl.I 

AQUILA  WING  AIRFOIL  COORDINATES  (PERCENT  OF  CHORD) 


Station 

Upper  surface 

Lower  surface 

(x/c--percent) 

(z/c--percent) 

(z/c--percent) 

0 

- 

0 

1.25 

3.34 

1.54 

2.5 

4.44 

2.25 

5 

5.89 

3.04 

7.5 

6.90 

3.61 

10 

7.64 

4.09 

15 

8.52 

4.84 

20 

8.92 

5.41 

25 

9.08 

5.78 

30 

9.05 

5.96 

40 

8.59 

5.92 

50 

7.74 

5.50 

60 

6.61 

4.81 

70 

5.25 

3.79  (modified) 

80 

3.73 

2.43  (modified) 

90 

2.26  (modified) 

1.00  (modified) 

95 

1.69  (modified) 

0.239  (modified) 

100 

1 .00  (modified) 

0.016  (modified) 
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Figure  VI.l 

Aquila  airfoil,  modified  NACA  23015 
With  Mean  Camber  Line  Curve  Fits 


Leading  edge  sweep  angle  = 
28.7”  V 


Figure  VI, 2 

Aquila  Wing  Geometry 
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An  analysis  was  performed  on  the  basic  wing  planform  using  a 
NASA-developed  vortex  lattice  FORTRAN  program,  named  "SUB", 
resident  on  the  NFS  VAX  computer  system  [Ref  29].  The  wing  was 
divided  into  100  panels  in  a  ten-by-ten  equally  spaced  grid.  A 
FORTRAN  program,  named  "camber.for"  and  shown  in  Appendix  C, 
converted  the  derivative  of  the  equation  for  the  mean  aerodynamic 
chord,  given  in  Figure  VI.  1,  to  a  local  incidence  angle  for  each  of  the 
100  panels.  The  output  of  the  program  "camber.for"  did  not  include 
the  three  deg  geometric  twist.  This  twist  was  calculated  manually 
for  each  panel.  The  final  input  file  to  program  "SUB"  is  also  given  in 
Appendix  C.  The  results  of  the  vortex-lattice  analysis  are  listed  in 
Table  VI.II. 

Table  VI.II 

RESULTS  OF  VORTEX  LATTICE  ANALYSIS  OF  AQUILA  WING 


CHARACTERISTIC 

VALUE 

Lift  curve  slope,  Cj^^ 

3.771  per  radian 

0.06581  per  deg 

Slope  of  pitching  moment/Cj^ 

curve, 

dCm/dCj^  =  CmcL 

-0.7578 

Intercept  of  Cm/Cj^  curve  for 

^macwb 

0.00761 
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C.  CANARD  SIZING  AND  LOCATION 


1 .  Static  Margin 

a.  General 

With  the  wing  characteristics  known,  a  suitable 
balancing  force,  in  the  form  of  a  canard,  could  now  be  designed.  A 
static  margin  of  15  percent  was  targeted  during  the  longitudinal 
stability  analysis,  to  provide  a  good  balance  between  stability  and 
pilch  response  [Ref  30].  Figure  VI.3  depicts  the  coordinate  system 
and  variables  used  during  the  static  margin  analysis.  Note  that  the 
variables  shown  on  Figure  VI.3  are  written  in  the  form  used  in  the 
static  longitudinal  stability  program  "xplot". 

b.  Computer  modelling 

The  static  margin,  essentially  the  distance  between  the 
whole  aircraft's  aerodynamic  neutral  point  and  its  center  of  gravity, 
was  determined  by  creating  a  MATLAB  program  named  "xplot".  The 
code  for  this  program  is  shown  in  Appendix  D.  The  algorithm  for  this 
program  follows: 

(1)  The  aircraft  neutral  point  was  assumed  to  be  due  to 
a  combination  of  two  items:  the  Aquila  wing  and  the  to-be-designed 
canard.  As  with  the  simulation  effort  described  in  Section  VII,  no 
aerodynamic  allowance  was  made  for  the  AROD  duct.  The  following 
equation  was  used  to  define  the  neutral  point  [Ref.  31]: 


Xac  alrcr^,  f  +  Vh  •  ^2=2. 


1- 


V 


d£ 


(VI.A) 
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where: 


^acainr^z  *  Longitudlnal  position  of  aircraft  aero  center  (neutral  point),  normalized  by  c 
^aewb  *  Aerodynamic  center  of  wing/body  combination 
c  s  Mean  aerodynamic  chord  of  wing  (32.2  in) 

V//  s  Tail  volume  coefficient 

=  2^ 
cS 
where: 

s  Longitudinal  position  of  canard' s  aero  center  (user  specified  in  "xplot") 
s  Area  of  canard  (a  vector  of  values  in  "xplot") 

S  s  Wing  reference  area  (29.2  ft^) 

3-fl  =  3  ■  dimensional  lift  curve  slope  of  canard 
s  Lift  curve  slope  of  wing/body  combination 

downwash  (assumed  zero) 
da 


(2)  With  reference  to  Figure  VI. 3,  allowances  were 
made  for  the  weight  and  location  of  the  following  components;  AROD 
duct  including  landing  gear,  Aquila  wings  (including  spars  and 
hardware  discussed  in  Section  V),  fuel,  electronics  pod,  canard  with 
servo  and  supporting  booms,  and  an  optional  ballast  box  at  the  tip  of 
the  canard's  supporting  boom,  A  full  description  of  these 
components  and  their  weights  in  contained  at  the  end  of  program 
"xplot".  Appendix  D. 

2.  Longitudinal  trimmability 

The  purpose  of  program  "xplot",  discussed  in  the  preceding 
section,  was  to  specify  the  canard  size,  given  a  desired  static  margin 


and  canard  location.  The  next  step  in  the  overall  longitudinal 
analysis  was  to  determine  if  this  canard  size  was  satisfactory  to  trim 
the  Archytas  over  its  expected  flight  range.  To  accomplish  this 


Figure  VI. 3 


Archytas  longitudinal  stability  dimensions 


portion  of  the  analysis  a  program  named  "incidence"  was  written  in 
MATLAB  to  depict  the  incidence  angles  required  to  trim  the  Archytas 
over  a  range  of  angles  of  attack.  A  printout  of  the  code  for  program 
"incidence"  is  contained  in  Appendix  D.  The  germane  equations  were: 

^moaircr(tft  ~ 3-D ' k  (VI.B) 

where: 


~  AifCTaft  pitching  moment  at  zero  deg  absolute  angle  of  attack 
=  Wing  pitching  moment  coefficient  at  zero 

deg  angle  of  attack  (from  vortex  -  lattice  analysis) 

4  =  canard  incidence  (measure  positive  with  leading  edge  down) 


and 


(h-hac^h)-VH 


where: 


(Vl.C) 


- =  Slope  of  aircraft  pitching  moment  versus  absolute  angle  of  attack  curve 

dec 

Cu,  =  Lift  curve  slope  of  wing  (from  vortex  lattice) 
h  =  distance  from  leading  edge  of  mean  aerodynamic  chord,  c, 
to  center  of  gravity,  normalized  by  c 

hac^  =  distance  from  leading  edge  of  c  to  wing/body  aerodynamic  center, 
normalized  by  c 

^Laca^d  =  Lift  curve  slope  of  canard 


The  program  "incidence"  simply  calculated  equations,  VLB  and 
Vl.C  over  a  range  of  canard  incidence  angles,  i^.  ,  and  displayed  the 


results  on  a  plot  of  pitching  moment  coefficient  versus  angle  of 
attack.  Typical  flight  angles  of  attack  were  calculated  for  a  cruise 
condition,  called  "CR"  (sea  level,  standard  day,  speed  of  135  ft/sec,  95 
Ibf  weight),  and  a  slow  flight  condition,  called  "PA"  for  "powered 
approach"  (speed  of  59  ft/sec).  These  angles  of  attack,  2.3  deg  and 
10.6  deg  respectively,  were  also  plotted  by  "incidence"  to  give  the 
user  a  graphical  display  of  the  incidence  angles  to  trim  the  Archytas 
over  a  range  of  angles  of  attack,  and  therefore  speeds. 

3.  Results 

In  earlier  thesis  work  the  electronics  pod,  containing 
components  for  servo  control  and  vehicle  stability  augmentation,  had 
been  located  below  the  engine.  Figure  VI. 4  shows  this  aft  pod 
location  as  "configuration  1".  The  original  AROD  had  the  electronics 
pod  located  as  shown  on  Figure  VI.4  and  labelled  as  "configuration 
2".  Initial  longitudinal  stability  estimations  for  the  Archytas 
configuration  were  made  with  the  electronics  pod  in  the 
"configuration  1"  position.  The  results  showed  that  an  Archytas  with 
ten  Ibf  of  electronics  that  far  aft  required  canard  support  booms  on 
the  order  of  8  ft  to  achieve  the  15  percent  static  margin  goal.  Even 
with  this  long  a  moment  arm  for  the  canard,  its  chord  was  less  than  3 
in  (remember  that,  due  to  installation  reasons,  the  canard  span  was 
considered  constrained  between  the  two  booms,  a  distance  of  31.5 
in).  Use  of  the  program  "incidence"  showed  that  the  aircraft 
configured  with  this  long  boom  and  small  canard  required  incidence 
angles  beyond  +.  20  deg  to  trim  the  vehicle  over  the  operational 
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airspeed  range. 

It  appeared  at  this  point  that  an  aft  electronics  configuration 
was  not  feasible.  However,  before  abandoning  the  aft  pod  placement, 
a  tradeoff  study  was  conducted  to  estimate  the  amount  of  ballast 
required  to  bring  the  center  of  gravity  for  configuration  1  far  enough 
forward  to  satisfy  stability  considerations.  The  study  compared  an 
aft  pod  (configuration  1),  with  ballast  placed  on  the  tip  of  the  boom, 
versus  a  pod  placed  forward.  Multiple  executions  of  programs 
"xplot"  and  "incidence"  were  conducted  to  achieve  similar 
longitudinal  stability  and  trimmability  characteristics.  Figure  VI. 5 
shows  pitching  moment  coefficient  plots  for  each  of  the  two  different 
configurations.  The  two  plots  on  Figure  VI.5  are  similar  in  that  the 
chord  of  the  canard  for  each  configuration  is  about  the  same,  12  in 
and  the  aircraft  can  be  trimmed  over  the  same  range  of  airspeeds 
with  similar  incidence  angles.  But,  the  similarity  ends  there  in  that 
the  aft  pod  configuration  required  7  Ibf  of  ballast!  This  increased  the 
weight  of  the  vehicle  a  corresponding  amount,  thereby  reducing  the 
thrust-to-weight  ratio  by  over  seven  percent,  with  no  added 
capability,  just  ballast.  As  discussed  in  the  simulation  section  of  this 
thesis,  that  kind  of  reduction  in  thrust  to  weight  would  likely  prove 
to  be  unacceptable.  Therefore,  the  aft  pod  configuration  was 
abandoned. 

The  forward  pod  configuration  was  then  considered  in 
detail.  To  minimize  aerodynamic  interference  between  the  pod  and 
the  canard,  a  minimum  separation  distance  of  one  canard  chord 
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between  the  two  was  chosen  as  a  design  constraint.  Figure  VI. 6 
shows  the  location  of  the  forward  pod  configuration  and  the 
separation  constraint.  Program  "xplot"  and  "incidence”  were  then  run 
for  a  range  of  canard  locations. 


Figure  VI.4 

Archyias  electronics-pod-location  options 
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Pitching  moment  coefficient  about  c.g.,  Cmcg  Pitching  moment  coefficient  about  c.g.,  Cmcg 


Absolute  angle  of  attack,  alpha  a,  deg 


Absolute  angle  of  attack,  alpha  a,  deg 


Figure  VI.5 

Static  Longitudinal  Stability  Comparison  Between  Pod  Configurations 


80 


Canard 


Scale  Drawing 


Figure  VI. 6 

Archytas  Forward  Electronics  Pod  Configuration--Deiail 


To  review,  the  constraints  and  goals  at  this  point  were; 

•  15  percent  static  margin  (analysis  tool:  "xplot") 

•  Minimum  weight  (analysis  tool:  "xplot") 

•  One  chord  length  separation  between  canard  and  pod 

•  Trimmable  over  flight  range  (analysis  tool:  "incidence") 

•  Canard  span  limited  to  31.5  in  by  installation  restrictions 

The  result  of  runs  for  various  canard  locations  (variable 

"Canardlocation"  in  program  "xplot")  on  the  two  programs  is 
summarized  in  Figure  VI. 7.  As  noted  on  the  figure,  with  a  canard 
located  52  in  in  front  of  the  AROD  duct  entrance,  all  constraints  were 
satisfied  and  the  weight  was  minimized.  This  canard  location  and 
size  was  therefore  selected.  Figures  VI.8  and  VI. 9  show  the  results 
of  running  programs  "xplot"  and  "incidence"  for  the  selected  canard 
configuration. 


Location  of  canard,  Canardlocation,  in  forward  of  AROD  duct 


Figure  VI. 7:  Summary  of  Longitudinal  Stability  Analysis 
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Figure  VI. 8:  Final  results  from  program  "xplot" 
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AOA  vs.  Cmcg  for  varying  incidence  angles 


33UI3  ‘  S  o  inoqc  uidpijjsoo  iiiomoiu  Siiiip)^ 


Figure  VI. 9:  Final  results  from  program  "incidence" 
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D.  CANARD  BOOM  STRUCTURAL  DESIGN 
1.  General 

Once  the  canard  had  been  sized,  and  placed,  by  the 
longitudinal  stability  and  control  analysis  described  in  previous 
paragraphs,  a  final  design  for  the  boom  to  support  the  canard  could 
be  performed.  An  initial  effort  was  made  to  manufacture  a  graphite- 
epoxy  tubular  boom  in  an  attempt  to  minimize  weight.  Finite 
element  analysis  for  composite  structures  was  unavailable  and  so  a 
trial-and-error  method  ensued.  A  test  boom  was  manufactured  by 
wrapping  graphite  in  an  epoxy  matrix  around  a  waxed  plastic 
mandrel.  The  resulting  structure  was  too  flexible  and  easily 
crushable  and  was  considered  too  weak  to  support  expected  loads. 
To  prevent  delay  of  other  aspects  of  the  project,  manufacturing  of  a 
composite  boom  was  abandoned  in  favor  of  aluminum.  A  tubular 
shape  was  chosen  as  a  good  aerodynamic  shape  and  to  allow  routing 
of  wiring  and  tubing  to  the  canard  servo  and  planned  sensors  such  as 
angle  of  attack  and  airspeed. 

The  installation  for  the  canard  boom  is  shown  in  Figure  V.l 
and  consisted  of  a  support  of  the  boom  by  a  hole  in  the  forward  spar 
and  a  base  for  the  boom  attached  to  the  aft  spar.  Because  the 
majority  of  the  boom  would  be  forward  of  the  forward  wing  spar,  the 
boom  was  modeled  as  a  cantilever  tubular  beam  as  show.i  in  Figure 
VI. 10. 
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Figure  VI.  10 

Archytas  Canard  Boom  Modeled  as  Cantilevered  Beam 

2.  Shear  stress  derivation 

Figure  VI.  11  shows  an  elemental  portion  of  this  tubular 
beam.  Summing  forces  in  the  x-direction: 

(T  +  dt)  t dx-rt dx  -  \OxdA-\-  \{ax-^dax)dA  =  0 

A  A 

X  in  the  first  two  terms  cancels  and  (J  from  the  last  two 
terms  cancels. 
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Figure  VI.ll:  Archytas  Canard  Boom  E  emental  Portion 
With  Forces  and  Coordinate  System  shown 

The  following  substitutions  are  then  made: 
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which  results  in: 


^'itax  -\-dL]  ^dA  =  0 


rfTt+J-^^=0 
A  ' 


Converting  to  polar  coordinates  (i.e.  y=rsin0  and  dA=:trd0)  results 


drt+Y  jrsin9trd6  =  0 

^  e 


Eliminating  t  from  the  above  equation  and  removing  all  non-0 
terms  from  the  integral  and  then  evaluating  the  integral  results  in: 


0  +  dd 

dr+^  jsinft/0  =  O 

dT+^[-cose]^-^^^  =0 


(VLB) 
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Consider  for  a  moment  the  bracketed  term  to  be  evaluated: 


[-cos0]^‘^‘^^  =-COS(0  +  i/0)+COS0 

=  -cos6cosd6+sinOsind0+cos9 
(making  a  small  angle  assumption :)  (Vi.C) 

=  -  cos  0 + sin  6lti0  -  cos  0 
=  sin0d6 


inserting  equation  VI.C  into  VLB  and  then  solving  for  dx  results  in: 


dr  = 


Fr^sinOdO 

/ 


integrating  both  sides  results  in: 

f  .  r  FF  sin  ddd 

jdr  =  -j - - - 

_  Fr^  COS0 


3 

the  moment  of  inertia  for  a  tube  (I=7tr  t)  is  then  inserted,  giving 
the  final  expression  for  shear  stress: 


FcosO 

nrt  (VI.D) 
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4.  Failure  criteria 

The  Von  Mises  stress  criterion  was  chosen  because  it  is  less 
conservative  than  the  Tresca  criterion  [Ref  32].  A  less  conservative 
approach  was  considered  appropriate  because  of  the  unmanned  and 
weight  critical  nature  of  the  project.  A  factor  of  safety  was  also  built 
in  because  material  analysis  was  based  on  the  yield  stress,  not  the 
ultimate  stress.  Following  is  the  development  of  the  Von  Mises 
stress  for  the  loading  shown  in  Figure  VI. 10  and  developed  in 
previous  paragraphs. 

aka:  X 

^  Von  Mises  P  P  ^  ®  ®  0 

=  -  (fy?  +  (fy  +(CT,  +6(^^ 
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Substituting  equations  VI.D  and  VI.E  gives  the  final  general 
expression  for  Von  Mises  stress: 


^  _  IF^x^  sin^  6  j_  6F^  cos^  6 

^vm  -,2.4.2  '  „2_2.2 


n  r  r 


n  rr 


(Vl.F) 


To  find  the  position  of  the  highest  Von  Mises  stress,  the 
derivative  of  equation  Vl.F  is  taken  with  respect  to  0  and  set  equal 
to  zero  and  constant  terms  eliminated.  The  result: 


dO 


(  2 
7  ^ 

v'* 


lsin0cos0)=O 


(VI.G) 


Which  is  satisfied  at  0  =  0,  90,  180  and  270  degrees—and  is 
termed  the  first  derivative  criterion.  To  determine  if  each  point 
satisfying  the  first  derivative  criterion  was  a  maximum  or  minimum 
point,  the  second  derivative  of  equation  Vl.F  with  respect  to  0  was 
taken.  The  result: 


Always  >  0 

/ 


-3 


COS  20 


>  or  <  or  =  0?? 


(VI.H) 
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The  maxima  will  occur  when  equation  VI. H  is  less  than  zero. 
The  first  factor  in  equation  VI. H  is  always  greater  than  zero  and  so  is 
excluded  from  the  analysis. 

The  second  term  is  less  than  zero  only  when 

(x^/r^  -  3)  <  0  (VI. I) 

The  largest  allowable  radius,  r,  due  to  physical  limitations  to 
the  installation  of  the  wing  spar  on  the  AROD,  was  about  three  inches. 
A  3.0  in  radius  substituted  satisfies  the  inequality  of  equation  VI. I 
only  when  x  is  less  than  5.2  in,  which  is  clearly  an  impractically  short 
boom  length  and  therefore  the  second  term  of  equation  VI. H  is  also 
excluded  from  the  analysis. 

The  range  of  0  that  will  make  equation  VI.H  less  than  zero 
are  given  as  follows: 

cos20  <  0  when: 

45  degrees  <0  <135  degrees 
or 

225  degrees  <0  <315  degrees 

This  range  of  angles  is  termed  the  second  derivative 
criterion.  There  are  only  two  angles  which  satisfy  both  the  second 
and  first  derivative  criterion:  90  and  270  degrees.  Substituting  90  or 
270  degrees  into  equation  VI.F  gives  the  final  expression  for  the 
maximum  Von  Mises  stress: 
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vtflmax.imum 


7t  r  t 


(VI.J) 


5 .  Canard  boom  loading  model 

Figure  VI.  12  depicts  the  loading  on  the  canard's  support 
boom  and  a  static  diagram  of  the  loads  and  moment  arms  involved. 
The  canard  lift  was  assumed  to  be  the  maximum  lift  that  an  airfoil  of 
the  canard's  size  could  generate.  This  value  for  maximum  lift  was 
calculated  as  follows: 


l^amrd  ^  4  =  O.SpV^S.Cj,  =  170  lb, 
where: 

p=0.00237688  slugs/ft^  (std  sea  level) 

V=100  knots=169  ft/sec 

Sc=canard  area=(31.5)(14)  in^ 

Cj  =canard  coeff.  of  lift 

=C’L««stall  =(2^X15  degrees) 

The  value  for  velocity  came  from  the  simulation  work.  The 
canard  area  was  the  result  of  the  longitudinal  stability  analysis.  The 
canard's  coefficient  of  lift  was  considered  as  a  worst-case  scenario. 
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6.  Material  selection 

Equation  VI.J  gives  the  Von  Mises  stress  as  a  function  of  the 
applied  force  (one  half  of  L  from  equation  VI.K--because  there  are 

two  booms  to  support  the  canard's  lift),  the  distance  from  force 
application  (a  maximum  value  of  59  inches  as  shown  in  Figure 
VI.  12),  and  the  radius  and  thickness  of  the  tube.  This  equation  was 
implemented  on  a  spreadsheet  for  31  commercially  available 
aluminum  tubing  configurations  [Ref.  33].  The  spreadsheet  is  shown 
in  Appendix  A,  Table  A. IV  and  the  Von  Mises  strength  of  each  of  the 
configurations  is  shown  on  Figure  VI.  13,  along  with  the  yield  stress 
values  for  the  two  different  types  of  aluminum  [Ref.  32]. 
Configurations  which  fall  below  the  appropriate  yield  line  are  those 
that  "passed"  the  Von  Mises  criterion.  The  configurations  that  passed 
this  criteria  were  then  compared  by  weight.  The  weight  per  linear 
foot  was  calculated  simply  by  knowing  a  configuration’s  volume 
times  its  density  [Ref  32],  as  shown  in  Appendix  A,  Table  A. IV.  The 
configuration  eventually  selected  was  the  lightest  one:  6061 -T6 
aluminum  with  a  2.5  inch  outer  diameter  and  a  wall  thickness  of 
0.035  inches,  as  shown  in  Figure  VI.  14. 

Prior  to  purchasing  the  selected  tubing  a  buckling  analysis 
was  conducted.  With  a  thickness  ratio,  r/t,  of  35,  the  chosen 
configuration  was  not  considered  a  "thin  walled"  structure  (the 
literature  focused  on  tubes  with  thickness  ratios  beginning  at  500) 
and  therefore  comparison  with  empirical  information  showed  that 
buckling  was  extremely  unlikely  [Refs  34  and  35]. 
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Configurations  below  yield  are  those 
predicted  to  not  fail 


(|sd)  8sej)s  PI»|A  s»ss|yy  uoa 


Figure  VI. 13 

Canard  support  boom  strength  comparison 
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Figure  VI.  14 

Canard  Boom  Configurations  Compared  By  Weight 

7 .  Structural  testing  of  boom  support  structure 

The  canard  booms  were  designed  to  be  supported  as  shown 
in  Figures  VI.  10  and  VI.  1.  They  were  manufactured  as  specified  in 
reference  21.  The  strength  of  these  supports,  the  carry-through  hole 
in  the  forward  spar  and  the  boom  base  on  the  aft  spar,  were  each 
tested.  The  tests  were  conducted  in  a  similar  manner  as  that 
described  in  Section  V  for  the  wing  spar.  The  load  was  applied 
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gradually  and  incrementally  using  a  two  weight  system.  Figure 
VI.  15  shows  the  test  spar  in  its  fixture.  The  two  separate 
components  tested  with  this  one  test  spar  are  shown:  the  canard- 
boom  carry-through  hole  and  the  canard-boom  base. 

With  reference  to  Figure  VI.  12,  a  simple  statics  problem 

showed  that  at  a  design  load  (canard  lift)  of  85  lb,  the  force  at  the 
boom  base  Cubase  Figure  VI. 12)  would  be  228  lb  and  the 

support's  reaction  ("^support  Figure  VI.12)  would  be  313  lb.  The 

boom  base  was  tested,  as  shown  on  the  upper  portion  of  Figure 
VI.  16,  to  243  percent  of  the  design  load  without  any  failure  or 
permanent  deformation.  Appendix  B,  Figure  B.6  is  a  photograph  of 
the  test  of  the  boom  base.  The  boom’s  forward-spar  carry-through 
hole  was  tested,  as  shown  on  the  lower  portion  of  Figure  VI.  16,  to 
207  percent  of  the  design  load  without  any  failure  or  permanent 
deformation.  Appendix  B,  Figure  B.7  is  a  photograph  of  the  test  of 
the  carry-through  hole.  The  boom  base  was  then  retested  by 
removing  it  from  the  vise  shown  in  Figure  VI.  15  and  adhering  it  to  a 
test  spar  with  standard  epoxy.  This  test  spar  was  constrained  and 
the  test  of  the  boom  base  repeated  with  the  same  results. 
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Test  boom  base 


Test  spar  a.  Boom  base  feft  soar)  structural  test 

with  hole 

through 


b.  Boom  support  hole  (forward  soar)  structural  test 


Figure  VI.  16:  Archytas  canard-boom  test  set-up 
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VII.  SIMULATION 


A.  GENERAL 

The  equations  of  motion  were  developed  for  a  three-degree-of- 
freedom  (3-DOF)  simulation  of  Archytas  motion— two  displacements 
(downrange  and  vertical)  and  one  rotation  (pitch).  The  simulation 
did  not  address  either  cross-range  displacement  or 
lateral/directional  rotations.  The  simulation  program,  named  AWAl, 
will  be  discussed  following  derivation  of  the  equations  of  motion  and 
explanation  of  the  modeling  of  germane  forces  and  moments. 
Several  assumptions  made  during  formulation  of  the  simulation  are 
contained  in  Appendix  E. 


B.  EQUATIONS  OF  MOTION-LINEAR  MOMENTUM 
1 .  Derivation  of  equations 

Figure  VII.  1  depicts  the  general  coordinate  system  and  the 
nomenclature  used  during  formulation  of  the  conservation  of  linear 
momentum  equations  of  motion  for  the  Archytas. 

With  reference  to  Figure  VII.  1,  the  first  two  equations  of 
motion  are  simple  vector  resolutions  of  the  velocity  vector  into  the 
inertial  reference  frame; 


x=  Vcosy 
h  =  Vsiny 


(VILA) 

(VII.B) 
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Figure  VII:  Coordinate  System  and  Nomenclature  for 
Archytas  Conservation  of  Linear  Momentum  Equations  of  Motion 


The  next  equation  of  motion  is  obtained  by  applying 
Newton's  second  law  for  conservation  of  linear  momentum  along  the 
Velocity  vector: 


^Falong  velocity 


vector 


(ni)(aCCClcrationajong  velocity  vector) 


Summing  forces  and  making  the  following  substitution: 
^  =  V,  results  in: 

Tcosa  -  D  -  mgsiny  =  mV  v. 
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gives: 


V/  Tcosa  D 
V  =  -  -  —  -  gsiny 

m  m  ^  ^  (VII.C) 

The  next  equation  of  motion  is  obtained  also  with  Newton's 
second  law,  this  time  applied  along  the  vertical  inertial  axis: 


=  rnh 


(VIl.D) 


Consider  for  a  moment  the  expression  h; 


h  -  d(h) 
dt 


and  with  h  =  Vsiny  (equation  VII.B) : 


h  = 


_  dCVsiny) 
dt 


=  vd(sinr)  siny  W 

dt  ^  dt 

=  Vcos/i''+  sinyV 


(VII.E) 


Substituting  equation  VII.E  into  equation  VIl.D  and 
summing  forces: 

Lcosy  -  Dsin  +  Tsin0  -  mg  =  m(Vcos7X  +  sinyV) 
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Solving  for  y  : 


•  _  Lcosy  Dsiny  ^  Tsin^  g  _  Vsiny 
^ "  mVcosy  mvcosy  mVcosy  Vcosy  Vcosy 

_  L  Dtany  ^  Tsin^ _ g  Vtany 

mV  mv  mVcosy  Vcosy  V 


Substitution  of  the  following; 

sin0  =  sin(y+  a)  =  sinycosa  +  sinacosy 


in  equation  VII.F  gives: 


.  _  L  Dtany  T(tanycosa-f  sing)  g  Vtany 
^  mV  mV  mV  Vcosy  V 

(VII.G) 


Equations  VII.A,B,C  and  G,  then  form  four  equations  of 
motion  for  displacement  and  velocity  in  two  axes.  The  variables 
solved  for  by  finite  differencing  these  four  equations  are  as  follows: 

X  =  downrange  distance 
h  =  height 

V  =  Speed  along  velocity  vector 
y  =  flight  path  angle 
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2.  Modeling  of  linear  momentum  equation  elements 
a.  Thrust  modeling 

The  thrust  is  modeled  based  on  engine  tests  (see  Section 
IV),  Data  from  these  tests  were  used  to  transform  a  user  input  of 
throttle  position  (on  a  scale  from  0%  for  idle  to  100%  for  full  power) 
to  engine  thrust.  Details  of  the  analysis  to  produce  this 
transformation  are  contained  in  Section  IV.  When  plotted,  the  model 
appears  as  shown  in  Figure  VII.2. 


Figure  VII.2 

Archytas  Throttle  Position  vs.  Thrust  model 
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b .  Lift  modeling 

The  lift  of  the  wing  was  modeled  as  a  function  of  angle 
of  attack.  The  pre-stall  relationship  between  the  trim  coefficient  of 
lift  and  the  trim  angle  of  attack  was  obtained  from  Aquila  wind 
tunnel  tests  [Ref.  28].  However,  no  information  was  available  on  the 
lift  curve  past  the  stall  angle  of  attack.  In  the  absence  of  high-angle- 
of-attack  lift-curve  data  for  the  Aquila  wing,  the  post-stall  lift  versus 
angle-of-attack  relationship  was  assumed.  The  relationship  between 
coefficient  of  lift  and  angle  of  attack  is  graphically  depicted  on  Figure 
VII. 3.  This  figure  has  two  parts:  part  (a)  shows  the  lift  curve  as  a 
function  of  angle  of  attack  in  radians;  because  no  one  curve  fit  would 
work  over  the  entire  range  of  angles  of  attack,  the  lift  curve  was 
broken  into  several  sections  and  polynomial  curve  fits  were 
developed  for  each  segment.  Part  (b)  is  the  same  plot  but  with  curve 
fits  removed  and  angle  of  attack  in  degrees. 

The  lift  was  then  calculated  from  the  formula: 

L  =  \pV'^SCl 

where: 

p  =  density 
V  =  speed 

S  =  reference  wing  area 
Cl  =  Coefficient  of  lift 
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Trim  Coefficient  of  Lift,  Cl 


SOURCE:  Aquila  wind  tunnel  data  for  unstalted  region;  assumed  profile  for  stalled  regions. 
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Trim  Angle  of  Attack,  ALPHA,  rad 

a.  Lift  model  with  curve  fits  shown 


Trim  Angle  of  Attack,  ALPHA,  deg 


b.  Lift  model  with  curve  fits  removed 


Figure  VII. 3  (two  graphs):  Archytas  Coefficient  of  Lift  Modeling 
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c.  Drag  modeling 

The  coefficient  of  drag,  Cj^  ,  is  traditionally  expressed  as 
a  function  of  the  coefficient  of  lift: 


Co=Cd„+KCl^ 


(VII. H) 


But  for  the  case  of  a  vehicle  potentially  operating  at  high  angles  of 

attack,  such  as  the  Archytas,  this  equation  predicted  a  decrease  in 

drag  at  angles  of  attack  beyond  stall.  Unfortunately  this  is  not  the 

reality  and  so  an  alternate  model  for  the  coefficieiu  of  drag  was 
assumed.  The  new  model  assumed  that  C|^  would  continue  to 

increase  as  a  second-order  function  of  angle  of  attack  (a).  A  further 
assumption  was  that  the  relationship  between  Cj^  and  a  would  be  the 

same  for  pre-stall  and  post-stall  regions,  allowing  use  of  equation 
VII.H  as  follows. 

From  Figure  VII.3a,  for  angles  of  attack  less  than  stall, 
the  following  relationship  exists  between  the  coefficient  of  lift  and 
angle  of  attack: 


Q  =  0.0952+ 3.38a 


(VII.I) 


The  following  constants  were  used: 


1 

;zi4Re 


0.0936 


(VII.J) 
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where: 


AR  s  aspect  ratio  =  4.25 
e  =  span  efficiency  factor  =  0.8 


C Do  =0.06 


(VII.K) 


Substituting  equations  VII.I  ,  VII.J  and  VII.K  into  VII. H 

results  in: 


Cd  =  0.0608  +  0.0602a + 1 .01 


(VII.L) 


Equation  VII.L  is  graphically  depicted  in  Figure  V1I.4: 
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Figure  VII.4:  Archytas  Coefficient  of  Drag  Model 
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The  drag  was  then  modeled  traditionally  as: 


D-^pV^SCo  (VII.M) 

C.  EQUATION  OF  MOTION-ANGULAR  MOMENTUM 
1 .  Derivation  of  equation 

Figure  VII. 5  depicts  the  general  coordinate  system  and  the 
nomenclature  used  during  formulation  of  the  conservation  of  angular 
momentum  equation  of  motion  for  the  Archytas. 


Figure  VII.5:  Coordinate  System  and  Nomenclature  for  Archytas 
Conservation  of  Angular  Momentum  Equations  of  Motion 
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With  reference  to  Figure  VII. 5,  the  equation  of  motion  for 
the  pitch  axis  is  developed.  Applying  Newton's  second  law  for 
conservation  of  angular  momentum: 


Pitch  Moments  about  c.g.  =  (Moment  of  Inertia)  *  (pitch  angular  acceleration) 


or. 


^ctWlNG  ^ctCANAKD  ^ cfVANE  ^ ctPITCH  RATE  ~  ^ yy^ 


(VII. N) 


solving  equation  VII.N  for  &  and  substituting  results 


in: 


•  _  ^ctW/NG  ^ctCANARD  ^ctVANE  ^ctPrTCH  RATE 

fyy 


(VII.O) 


With  reference  to  Figure  VII. 5, 
expressed  as  follows: 

Mowing  =  iLcosa){xcgyy) 
^CfCANARD  ~  ^^^ftcANARD  COS(X){XCg(-) 

^cgVASE~  ~ 

^cgPnCH  RATE  ~  ~ 


the  individual  moments  are 


(VII. Pa) 

(VII. Pb) 

(VII. Pc) 

(VII. Pd) 


Note  that,  unlike  Section  Vi's  formulation  for  static 
longitudinal  stability  where  dimensions  forward  of  the  leading  edge 
of  the  mean  aerodynamic  chord  were  negative,  distances  in 
equations  Vll.Pa  through  VII. Pd  are  absolute  distances  and  the  sign 
of  each  term  is  accounted  for  by  the  moment  convention  shown  on 
Figure  VII.5. 

Substituting  equations  Vll.Pa  through  VII. Pd  into  equation 
VII.O  gives  an  equation  of  motion  for  pitch  acceleration: 


Q  = 


_  Mqc  -  iLcosa)ixcg^)  +  cosa){xcgc )  - 


•yy 


(VII.Q) 

Combining  equation  VII.Q  with  the  following  equations, 
defines  pitch  angle,  rate  and  angular  acceleration. 


^  =  ^  (VII.R) 

0  =  a  +  y  (VII. S) 

2.  Modeling  of  Angular  Momentum  Equation  Elements 
a.  Pitching  Moment  due  to  wing 

The  lift  of  the  wing  was  modeled  as  previously 
described  (Section  VI).  With  reference  to  equation  Vll.Pa,  the 
pitching  moment  of  the  wing  about  the  center  of  gravity  was 
modeled  as  a  function  of  the  moment  about  the  aerodynamic  center, 
taken  to  be  the  quarter  chord  of  the  mean  aerodynamic  chord,  and 
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the  lift  times  the  distance  from  the  quarter  chord  to  the  center  of 
gravity.  The  moment  arm  for  the  lift  force,  xcgw.  was  a  function  of 
where  the  wing  was  installed  with  respect  to  the  center  of  gravity. 
The  wing's  longitudinal  position  was  limited  due  to  installation 
constraints. 

b.  Pitching  Moment  due  to  canard  deflection 

For  the  simulation,  the  canard  airfoil  was  assumed  to  be 
the  same  as  the  main  wing,  to  provide  lift  as  well  as  longitudinal 
control.  The  lifting  force  of  the  canard  was  modeled  based  on  the 
coefficient  of  lift  versus  local  canard  angle-of-attack  relationship 
shown  in  Figure  VII. 6.  The  canard  angle  of  attack,  as  shown  on 
Figure  VII,6,  was  a  combination  of  the  vehicle's  angle  of  attack  and 
the  canard  deflection.  Changes  in  this  angle  due  to  pitch  rate  were 
considered  but  discounted  because,  at  the  highest  pitch  rate 
observed  during  simulation  runs,  the  change  in  canard  angle  of 
attack  due  to  the  pitch  rate  was  less  than  ten  percent.  At  nominal 
pitch  rates  this  change  in  local  angle  of  attack  was  approximately  one 
percent.  Additionally,  the  assumption  was  made  that  upwash  from 
the  wing  had  no  effect  on  the  canard.  This  was  assumed  because  the 
canard  was  over  two  chords  in  front  of  the  wing  and  the  ducted-fan 
was  assumed  to  negate  some  of  the  interaction  between  the  wing  and 
the  canard. 
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Local  canard  angle  of  attack,(ALPHA-Delc),  deg 
Figure  VII.6 

Canard  coefficient  of  lift  versus  local  angle-of-attack  model 


c.  Pitching  Moment  due  to  vane  deflection 

The  pitching  moment  due  to  vane  deflection  was 

modeled  based  on  engine  test  data.  These  tests  produced  rolling 

moment  data,  which  were  analyzed  to  predict  pitching  moment 

effectiveness.  Details  of  this  analysis  are  contained  in  Section  IV.  As 
discussed  in  Section  IV,  the  results  of  the  engine  tests  and  analysis 
showed  essentially  no  change  in  control  power  for  thrust  levels  from 
approximately  90  to  120  lb  (approximate  thrust-to-weight  ratios  of 
0.9  to  1.2),  which  was  in  agreement  with  tests  performed  by  the 
AROD’s  manufacturer  [Ref  18].  Because  the  thrust  levels  required 

during  hovering  and  initial  transition  were  anticipated  to  be  in  this 
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constant-control-power  range  of  90  to  120  lb  and  with  the 
assumption  that  airflow  over  the  vanes  would  remain  approximately 
constant  following  transition  to  horizontal  flight,  with  increases  in  the 
local  velocity  at  the  vane  due  to  forward  flight  speed  balanced  by 
throttle  decrease,  control  power  due  to  vane  deflection  was  modeled 
as  constant.  The  relationship  between  vane  deflection  and  pitching 
moment,  for  all  speeds  and  throttle  settings,  was  therefore  modeled 
as  shown  in  Figure  VII. 7.  The  vanes  are  assumed  to  produce  zero 
moment  at  zero  deflection.  As  discussed  more  fully  in  Section  IV, 
engine  tests  showed  a  swirl  of  the  airflow  in  the  vicinity  of  the 
vanes,  which  resulted  in  some  moment  being  produced  by  the  vanes 
at  zero  deflection.  This  swirl  effect,  however,  was  difficult  to  predict 
and  no  tests  were  performed  to  analyze  the  effect  of  two  vanes, 
deflected  in  a  pitching  sense,  on  the  airflow  patterns  near  the  vanes. 
The  vane  control  power,  therefore,  was  referenced  to  an  arbitrary 
neutral  position. 
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Figure  VII.7 

Control  Vane  pitching  moment  model 


d.  Pitching  moment  due  to  pitch  rate 

The  moment  due  to  pitch  rate,  also  known  as  pitch 

damping,  was  modeled  conventionally  as  per  equation  Vll.Pd,  using 
the  dimensional  derivative  M^,  defined  as  follows  [Ref.  36]: 

’  2V  (VII.T) 

where: 

^mq  s  pitching  moment  coefficient  due  to  pitch  rate 

q  =  dynamic  pressure  (do  not  confuse  with  q,  pitch  rate) 

S  s  Wing  area 
c  =  mean  aerodynamic  chord 
V  s  speed 
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With  no  better  data  available,  the  value  for  used 

was  that  given  for  the  Aquila  airframe  of  -1.47  per  radian.  [Ref.  28]. 
D.  COMPUTER  CODE 

1.  Program  features 

A  main  program,  called  "AWAl",  and  several  function 
programs  (similar  to  subroutines)  were  written  on  MATLAB®  to 

perform  simulation  of  the  three-degree-of-freedom  equations  of 
motion.  A  MACINTOSH®  IIX  computer  was  used  to  ease 
programming,  data  analysis  and  processing.  The  program  included 
the  following  features: 

•  Plot  of  vehicle's  downrange  and  height  position  as  well  as 

pitch  angle  were  graphically  displayed  in  an  easily  interpretable 

manner.  Data  to  enable  replay  of  flight  were  stored  and  replay  was 
available. 

•  Initial  conditions  could  be  specified  by  user,  or  default 
values  were  easily  input. 

•  Program  was  suspended  at  user's  discretion  so  that 

changes  in  control  inputs  could  be  made  if  desired. 

•  Plots  of  all  stored  variables  were  available  for  later 

analysis. 

•  Vehicle's  flight  path  and  all  aerodynamic  data  were 
stored  for  later  analysis  or  re-use. 

2.  Program  description 

By  way  of  introduction,  the  programs  and  their  basic 
functions  are  listed  below: 
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•  "AWAl  "--performed  variable  initialization,  file 
management,  numerical  differencing,  input/output  control  and 
overall  simulation  management. 

•  "archytas3"--stored  aircraft  constants,  e.g.,  wing  area, 
canard  location  relative  to  center  of  gravi,._y  pitch  damping 
coefficient. 

•  "yprimes4"— contained  equations  of  motion  VILA,  VII. B, 
VII.C,  and  VII.G;  calculated  the  derivatives  of  downrange  position, 
height,  speed,  and  flight  path  angle. 

•  "control3"--contained  the  model  for  the  conservation  of 
angular  momentum,  described  in  equations  VII. Q  and  VII.R.  Output 
was  updated  pitch  angle,  rate  and  acceleration. 

•  "THROTTLE"--converted  user  input  from  throttle 
position,  in  percent,  to  thrust  in  Ibf,  based  on  engine  test  data. 

•  "labeller""provided  labels  for  plots  requested  by  the 

user. 

•  "plotter" -provided  a  series  of  eight  plots  of  control  input 
and  vehicle  response.  Used  following  program  execution. 

•  "Movie"— replayed  the  vehicle's  flight  path  and  pitch 
angle  from  most  recent  program  execution  (or  a  saved  simulation 
run). 

A  flowchart  of  the  program  and  several  of  its  subroutines  is 
given  in  Appendix  F.  Appendix  G  provides  a  listing  of  the  code  for  all 
simulation  programs.  The  basic  program  scheme,  as  depicted  on  the 
flowchart,  was  as  follows: 
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1)  Aircraft  constants  were  imported  from  "archytas3". 

2)  The  user  was  prompted  for  how  many  iterations  were  to 
be  performed,  what  time  increment  to  use,  vehicle  weight,  control 
mode  (i.e.,  use  vanes  or  specify  a  fixed  pitch  rate)  and  whether  or  not 
the  user  wanted  to  suspend  the  program  to  make  control  inputs  (if 
answer  to  last  question  was  "yes"  then  user  was  asked  "how  many 
seconds  between  suspensions?").  The  user  also  had  the  option  of 
typing  "1"  which  input  default  values. 

3 )  The  user  was  prompted  for  vehicle's  initial  height,  speed, 
pitch  attitude  and  flight  path  angle. 

4)  Initial  derivatives  of  the  information  passed  in  step  3) 
were  calculated.  Initial  drag  was  calculated  based  on  drag  model. 

5)  The  output  file,  called  "OUTFILE"  which  stored  virtually 
all  important  values,  was  initialized.  This  file  was  updated  each 
time  step  to  allow  later  plotting  of  all  variables. 

6)  The  following  vector  of  primary  variables  was  defined: 

downrange  position 
h  _  height 

=  speed 

y\  flight  path  angle 

The  equations  of  motion  for  downrange  distance,  height, 
speed  and  flight  path  angle  were  solved  using  predictor-corrector 
finite  differencing.  Euler  first  forward,  half  step  and  Richardson 
Extrapolation  were  combined  as  follows  [Ref  37]: 
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Yil\  =  Y.+Att 

Y<‘\  =n+^  y„ 

n+j  2 

L‘ 

y  =9y(2)_ya) 

7)  The  user  was  presented  with  a  picture  of  the  vehicle's 
motion,  as  shown  in  the  sample  output.  Figure  VIL8.  This  display 
provided  the  user  with  a  display  of  the  aircraft's  position,  displayed 
as  a  circle,  and  also  pitch  attitude,  displayed  as  a  line  drawn  through 
the  circle.  The  user  also  received  speed  cues  by  observing  the  space 
between  successive  positions  in  this  constant-time-increment 
process. 


Figure  VII. 8:  Sample  Plot  From  Program  AWAl— Continuously 
Updated  Position  and  Pitch  Attitude 
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8  At  user-specified  intervals,  the  user  was  then  allowed  to 
view  a  plot  of  any  variables  desired  and  to  make  control  input 
changes.  Plots  were  presented  as  requested.  The  user  had  control 
over  canard  and  vane  deflections  and  throttle  position. 

9)  The  program  returned  to  step  5)  and  repeated  this 
process  until  the  user-specified  program  run  time  had  elapsed  or  the 
vehicle  crashed  (height  less  than  zero). 

10)  Following  program  execution  the  user  could  type 
"plotter"  which  provided  eight  plots  on  two  pages.  The  plots 
provided  (all  versus  time  except  last  plot): 

•  angle  of  attack 

•  speed 

•  pitch  attitude 

•  pitch  rate 

•  canard  deflection 

•  vane  deflection 

•  throttle  position 

•  height  versus  downrange  position 

11)  The  user  also  could  type  "Movie"  and  be  presented  with 
a  replay  of  the  vehicle's  flight  path. 

12)  Ancillary  programs  included: 

•  "Mplot"--to  produce  a  plot  of  all  pitching  moment 
equation  values  versus  time.  Used  for  troubleshooting. 

•  "filesaver"— to  save  program  runs  for  future  analysis 

•  "resurrect"— to  re-load  files  saved  by  "filesaver" 


121 


E.  RESULTS  AND  DISCUSSION 
1.  Time  step  determination 

The  optimum  time  step  for  use  during  runs  of  program 
"AWAl”  was  determined  by  varying  the  time  increment,  called 
"DELTAT"  in  the  program,  and  observing  simulation  behavior.  To 
speed  the  simulation  a  large  time  step  was  desirable  but  too  large  a 
time  increment  destabilized  the  numerical  differencing,  producing 
erroneous  results.  Time  increments  of  0.2,  0.1,  0.05  and  0.025  were 
compared  while  performing  identical  maneuvers.  In  each  case,  the 
maneuver  was  initialized  with  the  vehicle  in  a  vertical  attitude  at  a 
speed  of  1  ft/sec.  After  allowing  the  vehicle  to  accelerate  for  a  few 
seconds,  a  gentle  pushover  transition  was  initiated.  Figure  VII. 9 
depicts  a  comparison  of  the  vehicle’s  angle  of  attack  response  with 
varying  time  increments.  As  shown  in  the  figure,  at  time  increments 
of  0.2  and  0.1  the  angle-of-attack  response  showed  a  tendency  to 
"spike".  This  was  true  after  any  maneuver  was  initiated.  In  fact,  at  a 
0.2-sec  time  increment  the  simulated  vehicle  would  always  depart 
controlled  flight,  as  defined  by  angles  of  attack  exceeding  90  deg  in 
one  to  two  seconds.  As  shown  in  Figure  VII. 9,  the  response  was 
smooth  at  time  increments  less  than  or  equal  to  0.05  sec.  This 
smoothness  was  evident  during  all  maneuvers.  However,  the 
simulation  ran  noticeably  slower  with  the  0.025  increment,  with  no 
significant  differences  in  response.  As  a  result,  0.05  sec  was  chosen 
as  the  optimum  time  increment  to  balance  stability  of  the  numerical 
scheme  with  program  run-time. 
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Figure  Vn.9 

Archytas  Simulation:  Comparison  of  Time  Increments 
For  Identical  Maneuvers 

2 .  Takeoff,  transition  to  forward  flight 
a.  Thrust-to-weight  envelope 

In  order  to  characterize  the  effects  of  increasing  weight 
on  the  vehicle,  simulations  were  conducted  at  a  range  of  thrust-to- 
weight  ratios  (T/W).  In  theory,  any  thrust-to-weight  ratio  greater 
than  1.0  is  acceptable,  given  an  extremely  robust  control  system  and 
the  time  to  wait  for  a  given  flight  condition  to  be  reached.  However, 
the  reality  pilots  understand  is  that  more  thrust  is  better;  the 
question  is,  "how  much  is  enough?”.  Several  assumptions  were  made 


123 


when  considering  a  minimum  T/W. 

The  thrust  was  assumed  to  be  a  constant,  as  delineated 
in  Section  IV,  and  the  weight  was  varied.  Equivalent  takeoff 
maneuvers  were  simulated  at  full  thrust  with  weights  of  115,  110, 
105  and  100  lb,  corresponding  to  T/W  values  of  1.04,  1.09,  1.14,  and 
1.2  respectively.  The  maneuver  consisted  of  starting  in  a  vertical 
climb  at  1  ft/sec,  accelerating  straight  up  to  5  ft/sec  and  then 
beginning  a  very  gentle  "push-over",  using  3  deg  of  vane  deflection. 
The  maneuver  was  considered  complete  when  the  vehicle  had 
achieved  stable  forward  flight  with  less  than  10  deg  angle  of  attack. 

The  figures  of  merit  for  this  T/W  comparison  were 
flight  path  and  angle  of  attack.  Flight  paths  were  compared  based  on 
the  altitude  change  from  the  point  the  "pushover"  was  begun.  Angle 
of  attack  profiles  for  each  T/W  level  were  also  compared.  It  was 
assumed  that  high-angle-of-attack  (high-a)  regions  were  to  be 
avoided  due  to  the  difficult-to-control  and  unpredictable  nature  of 
flying  qualities  most  airplanes  exhibit  at  high  a.  Angles  of  attack 
above  20  deg  were  considered  unacceptable  for  a  vehicle  controlled 
remotely  by  a  person  without  the  proprioceptive  and  visual  cues  that 
the  pilot  of  a  manned  airplane  has. 

Of  the  two  figures  of  merit  for  this  analysis,  the  angle  of 
attack  became  the  primary  parameter  for  comparison.  The  flight 
paths  at  different  T/W  were  all  generally  acceptable.  In  other 
words,  at  all  T/W,  owing  to  the  gentle  nature  of  the  pushover,  the 
vehicle  did  not  lose  altitude  during  the  transition.  The  vehicle  was 
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only  about  75  feet  higher  at  the  highest  T/W  than  it  was  for  a 
comparable  maneuver  at  the  lowest  T/W.  A  comparison  of  the  angle 
of  attack  profiles  for  the  various  T/W  is  shown  in  Figure  VII.  10.  At 
T/W  of  less  than  1.09  the  vehicle  exceeded  the  established  angle-of- 
attack  limit. 

Because  the  vehicle  weight  predicted  by  the  analysis  of 
Section  VI  was  approximately  100  lb,  a  T/W  value  of  1.2  was  utilized 
for  all  subsequent  simulations. 


Figure  VII. 10 

Archytas  Simulation:  Thrust-To-Weight  Comparison 


b.  Takeoff  transition  tradeoffs 

Over  500  runs  of  program  AWAl  were  performed  for 
the  transition  from  a  hover  to  horizontal  flight  and  several  basic 
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profiles  emerged.  These  profiles  and  references  for  the  resulting 
plots  are  given  in  Table  VII.  For  each  maneuver  the  corresponding 
plots  of  downrange  distance  versus  height  ("X-H  PLOT")  and  time 
history  ("PARAMETER  PLOT")  are  referenced. 

Table  VILI 

ARCHYTAS  SIMULATION:  TAKEOFF  TRANSITION  PROFILES, 
DESCRIPTION  AND  LOCATION  OF  RESULTS 


No. 

MANEUVER 

description(I) 

X-H 

PLOT 

PARAMETER 

PLOT 

T" 

Low  vane 
deflection 
pushover 
from  hover 

From  hover,  3®  Delv  until 
10°  pitch  attitude,  then 
pull  up. 

Figure 
VII. 11 

Appendix 
H,  Figure 
H.l 

2 

High  vane 
deflection 
pushover 
from  hover 

From  hover,  10°  Delv 
until  0°  pitch  attitude, 
then 

pull-up  to  level  flight. 

Figure 
VII. 12 

Appendix 
H,  Figure 
H.2 

3 

Low  speed, 
low  altitude 
climbing 
gentle 
pushover 

From  hover  accelerate 
straight  up  until  V=5 
ft/sec,  then  3° 

Delv  until  45°  pitch 
attitude  then  pull-up  to 
level  flight. 

Figure 
VII. 13 

Appendix 
H,  Figure 
H.3 

1 

Vertical 
acceleration 
to  gentle 
pushover 

From  hover,  accelerate 
straight  up  until  V=20 
ft/sec  then  3°  Delv  until 

45°  pitch  attitude,  then 
accelerate  and  climb. 

Figure 
VII. 14 

Appendix 
H,  Figure 
H.4 

Note:  (1)  "Delv"  is  vane  deflection 
(2)  "V"  is  speed 
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Figures  VII.  11  through  VII.  14  depict  the  results  of  a 
particular  profile.  The  pull-up  maneuver  used  consistently  was  a 
medium-effort  pull-up  defined  by  use  of  10  deg  of  canard  deflection 
and  0  deg  of  vane  deflection.  Appendix  H,  Figures  H.l  through  H.4, 
give  detailed  time  histories  of  several  parameters  as  well  as  control 
inputs. 

The  first  profile,  in  which  the  vehicle  was  pushed  over 
using  a  small  vane  deflection,  was  characterized  by  over  100  ft  of 
altitude  loss  and  an  angle  of  attack  that  peaked  at  15  deg.  The 
second  profile,  similar  to  the  first  except  that  the  vane  was  deflected 
about  three  times  as  much,  was  marked  by  more  altitude  loss  (more 
than  200  ft)  and  about  the  same  angle-of-attack  response.  The  third 
profile  involved  initiating  a  gentle  pushover  while  climbing  but  at  a 
relatively  slow  speed.  The  third  profile’s  angle  of  attack  increased 
both  during  the  initial  part  of  the  transition  and  also  during  the  pull- 
up.  There  was  a  slight  loss  in  altitude  from  the  apex  of  the 
transition.  The  fourth  profile,  similar  to  the  third  except  the  vehicle 
was  allowed  to  accelerate  to  a  higher  speed  prior  to  pushing  over, 
was  characterized  by  very  low  angle  of  attack  and  a  continuous  climb 
requiring  less  than  75  ft  of  vertical  distance  to  complete. 

The  fourth  profile  was  the  most  repeatable  and  had  the 
following  advantages; 

•  The  angle  of  attack  was  low  during  the  entire  profile 
and  was  easily  controlled  with  slight  changes  in  vane  deflection 
during  the  initial  part  of  the  transition  and,  as  speed  built,  by  canard 
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deflection  during  the  latter  portion.  Controlling  angle  of  attack 
during  the  other  profiles  was  not  as  carefree,  requiring  constant 
attention  to  ensure  that  a  high  angle  of  attack  region  was  not 
entered. 

•  Altitude  required  to  complete  the  last  profile  was  well 
less  than  100  ft  and  it  was  all  altitude  gain  .  By  contrast,  the  other 
profiles  involved  altitude  losses  up  to  200  ft.  In  particular, 
performing  high-vane-deflection  pushovers  similar  to  the  first 
profile  generally  involved  having  to  climb  the  vehicle  in  a  hover  to 
an  extremely  high  altitude  (up  to  500  ft)  to  ensure  enough  distance 
below  to  complete  the  maneuver  without  striking  the  ground.  The 
altitude  required  for  these  high-vane-deflection  pushovers  was  very 
unpredictable. 

Profiles  1,  2,  and  3  had  no  apparent  advantages  and 
therefore  the  fourth  profile  was  recommended. 
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Weight  =  100 

deflection  pushover  from 

hover" 

Description: 

from  hover,  3®  Delv 

until  1 0®  pitch  attitude,  then 

"pull-up"  to  level  flight. 
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igure  VII.ll;  Archytas  Simulation:  Takeoff  Transition  Profile  #1 
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Delta  T  =  0.05 
Weight  =  100 


Maneuver:  "High  vane 
deflection  pushover  from 
hover" 

Description: 
from  hover,  1 0°  Delv 
until  0°  pitch  attitude,  then 
"pull-up"  to  level  flight. 
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Note:  scales  not  equal 
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Figure  VII.  12: Archytas  Simulation:  Takeoff  Transition  Profile  #2 
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Accelerating  straight  up,  3° 
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then  "pull-up"  to  level  flight. 
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Figure  VII.  13:  Archytas  Simulation:  Takeoff  Transition  Profile  #3 
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Maneuver:  "Vertical  acceleration  to 
gentle  pushover". 

Description: 

from  hover,  accelerate  straight  up  until 
V=20  ft/sec  then  3"  Delv  until  45“  pitch 
attitude,  then  accelerate  and  climb. 


Note:  scales  not  equal 


Delta  T  =  0.05 
Weight  =  100 


50  100  150  200 

Downrange  distance,  x,  ft 


250 


300 


Figure  VII.  14:  Archytas  Simulation:  Takeoff  Transition  Profile  #4 
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3.  Horizontal  flight 
a.  Trim  conditions 

The  horizontal  flight  characteristics  of  the  Archytas 
were  simulated  at  two  separate  trim  conditions  listed  in  Table  VII.II. 
These  were  considered  the  standard  conditions  for  all  horizontal 
flight  analysis.  Condition  1  was  considered  an  average  cruise 
condition  for  employment  of  the  vehicle  and  Condition  2  was  tested 
to  compare  with  Aquila  data  at  similar  flight  conditions. 


Table  VII.II 

ARCHYTAS  SIMULATION  HORIZONTAL  FLIGHT  TRIM  CONDITIONS 


PARAMETER 

CONDITION  I 

CONDITION  2 

Altitude 

150  ft 

150  ft 

Airspeed 

153  ft/sec 

103  ft/sec 

Angle  of  attack 

— 

3.2  deg 

Pitch  attitude 

IBHRHI 

3.2  deg 

Canard  deflection 

-7.2  deg 

-9.2  deg 

Vane  deflection 

0 

0 

Throttle  position 

9.5  percent 

Thrust 

49.6  Ibf 

25.5  Ibf 
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b.  Cruise  flight  short  period  characteristics 

The  short  period  characteristics  were  explored  in 
horizontal  flight  in  Condition  1,  as  specified  by  Table  VII. II,  with  a 
sinusoidal  pitch  doublet  input.  Figure  VII.  15  depicts  the  input,  the 
response  and  an  expanded  view  of  the  response  along  with  the 
analysis  of  the  dynamics. 

c.  Predicted  short  period  characteristics 

In  order  to  compare  the  simulation's  behavior  to  that  of 
the  Aquila,  the  short  period  characteristics  of  the  Archytas  were 
explored  in  horizontal  flight,  at  a  slower  speed:  Condition  2  as  listed 
in  Table  VII.II.  The  simulation  was  excited  by  a  sinusoidal  doublet 
near  the  natural  frequency.  Figure  VII.  16  depicts  the  input,  the 
response  and  an  expanded  scale  view  of  the  response  along  with 
analysis  of  the  dynamics.  These  are  compared  in  Table  VII.III  to  the 
Aquila's  short  period  characteristics,  as  derived  from  wind  tunnel 
data  and  analytical  methods  [Ref  38]. 
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Angle  ofattack,  ALPHA,  deg  Angleof  a.«ck,ALPHA,deg  Canard  Dell, Dele.deg 


Time,t,sec 

Figure  VIL15:  Predicted  Short  Period  Response  in  Cruise  Condition  1 
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Table  VILIII 

SHORT  PERIOD  CHARACTERISTICS: 
COMPARISON  BETWEEN  AQUILA  AND  ARCHYTAS 


CHARACTERISTIC 

aouila(^) 

ARCHYTAS^^) 

Speed  (ft/sec) 

101 

103 

Gross  weight  (lb) 

130 

100 

■SllllllillHi 

21 

-6.5 

Period  (sec) 

1.06 

0.798 

0.54 

0.40 

Natural  frequency 

(rad/sec) 

5.9 

8.57 

Notes:  (1)  Wind  tunnel  and  analytical  analysis  [Ref  28] 


(2)  Simulation  program  AWAl 


Angle  of  attack,  ALPHA,  deg  Angleof  attack,ALPHA,dcg  Canard  Defl,Delc,deg 


Time,t,sec 


Figure  VII,  16:  Archytas  Short  Period  Response  at  Condition  2 
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In  view  of  the  use  of  the  Aquila's  pitch  damping 

coefficient  in  the  Archytas  simulation,  the  similarity  of  damping 

ratios  between  the  two  is  a  favorable  indication  of  simulation 

fidelity.  The  difference  between  the  two  aircraft’s  period/frequency, 
a  strong  function  of  the  derivative  can  be  attributed  to  the  15 

percent  static  margin  of  the  Archytas  compared  to  that  of  the  Aquila, 
which  has  a  static  margin  on  the  order  of  five  percent  [Ref .28]. 
d.  Long  period 

The  long  period  or  phugoid  characteristics  of  the 
Archytas  were  simulated  at  both  trim  conditions  listed  in  Table 
VII. III.  The  phugoid  was  excited  by  slowing  the  Archytas  from  a 
level  flight  trim  condition  by  approximately  15  ft/sec  and  then 
returning  the  controls  to  the  trim  point  setting  [Ref.  39].  Figure 
VII. 17  depicts  the  control  input,  the  airspeed  response  and  an 
expanded  plot  of  the  airspeed  response  with  phugoid  characteristics. 
Table  VII. IV  shows  a  comparison  of  phugoid  characteristics  between 
Archytas  simulation  and  Aquila. 

The  simulation's  phugoid  period  was  within  nine 
percent  of  the  following  approximation  for  the  period  [Ref.  36]: 


'  phugoid 


=  V27r 

=  V2;r 


V 

8 

104 

32.2 


=  14.34 


The  simulation's  long  period  characteristics  also 
compared  well  to  the  information  available  for  the  similarly  shaped 
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Aquila.  The  only  significant  disparity  between  predictions  for  the 
two  vehicles  involved  damping  ratio.  The  difference  for  the  two 
aircraft's  damping  ratios,  a  strong  inverse  function  of  the  lift-to-drag 
ratio  (L/D),  may  be  explained  by  either  a  real  or  modelled  difference 
in  lift  or  drag.  In  as  much  as  the  lift  model  was  based  on  a  vortex 
lattice  analysis  that  agreed  well  with  data  available  for  the  Aquila,  a 
likely  reason  for  the  difference  in  damping  ratio  was  the  drag  model. 
However,  both  sets  of  long  period  characteristics  are  representative 
of  aircraft  flying  at  the  simulation  speed  and  will  pose,  at  worst,  only 
a  nuisance  to  a  control  system  or  pilot  attempting  to  counter  the 
phugoid. 

Table  VII.IV 

LONG  PERIOD:  COMPARISON  BETWEEN  AQUILA  AND  ARCHYTAS 


CHARACTERISTIC 

AQUILA^ 

ARCHYTAS^^^ 

Speed  (ft/sec) 

101 

104 

Gross  weight  (lb) 

130 

100 

CG  (percent  MAC) 

21 

-6.5 

Period  (sec) 

15.8 

16.9 

0.044 

0.16 

Natural  frequency 

(rad/sec) 

0.37 

0.40 

Notes:  (1)  Wind  tunnel  and  analytical  analysis  [Ref  28] 


(2)  Simulation  program  AWAl 
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4.  Transition  from  horizontal  flight  to  vertical  landing 
a.  General 

The  transition  from  horizontal  flight  to  a  hover,  for  the 
purpose  of  a  vertical  landing,  was  evaluated  during  approximately 
150  runs  of  the  simulation  program  AWAl.  Two  types  of  transitions 
were  attempted,  described  as  follows: 

•  A  level  deceleration  to  hover  parameters.  From 
stabilized  horizontal  flight,  the  profile  involved  attempting  to  slow 
the  vehicle  without  gaining  altitude  while  progressively  increasing 
the  pitch  attitude  to  90  degrees.  This  profile  necessitated  the  ability 
to  generate  high  angles  of  attack  during  the  deceleration  to  avoid 
gaining  altitude  as  the  pitch  attitude  was  increased. 

•  A  "zooming"  transition  to  hover  parameters.  From 
stabilized  horizontal  flight  the  vehicle's  pitch  attitude  was  increased 
while  maintaining  unstalled  angles  of  attack.  This  profile  involved  a 
climbing  or  "zooming"  flight  path.  As  the  vehicle  approached  high 
pitch  attitudes,  the  power  was  reduced  from  the  level  flight  setting 
to  allow  the  vehicle  to  decelerate.  Pitch  inputs  were  then  made  to 
capture  approximately  90  degrees  for  both  pitch  attitude  and  flight 
path  angle. 

For  standardization,  the  vehicle  was  considered  to  have 
entered  a  hover  when  the  airspeed  was  less  than  0.5  ft/sec,  the  pitch 
angle  and  flight  path  angle  were  within  5  degrees  of  90  degrees  and 
the  pitch  rate  was  less  than  2  deg/sec. 
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b.  Level  deceleration  profile 

As  is  common  with  canard  designs,  wherein  the  canard's 
aerodynamic  loading  is  higher  than  the  wing's,  the  Archytas 
simulation  could  not  be  slowed  below  approximately  65  ft/sec  (38 
KTAS)  in  level  flight.  Attempts  to  further  slow  the  airspeed  resulted 
in  a  pitch  attitude  decrease  with  an  accompanying  increase  in 
airspeed  and  loss  of  altitude—a  classic  canard  stall.  This  occurred 
despite  full-effort  control  inputs  for  both  the  canard  and  the  vanes. 
Figure  VII.  18  shows  a  situational  plot  of  an  attempt  to  slow  the 
vehicle  with  full  effective  canard  and  vane  deflections.  Appendix  H, 
Figure  H.5  displays  the  aircraft  response  and  control  inputs  for  this 
attempt.  The  canard  was  deflected  to  produce  the  maximum 
coefficient  of  lift  for  the  canard,  which  occurred  at  approximately  15 
degrees  relative  angle  of  attack.  As  shown  on  Figure  VII.5,  the  local 
angle  of  attack  at  the  canard  was  a  combination  of  the  aircraft's  angle 
of  attack  and  the  canard  deflection.  The  vane  was  deflected  to  30 
degrees,  the  maximum  effective  deflection.  As  shown  in  both  Figures 
VII.  18  and  Appendix  H,  Figure  H.5,  the  vehicle's  pitch  attitude 
simply  could  not  be  slowed  below  the  65  ft/sec  range  because  the 
canard  would  stall  prior  to  stalling  of  the  wing.  Figure  VII.  19  shows 
the  balance  of  moments  that  prevented  the  pitch  attitude  from  being 
increased.  This  figure  shows  that  the  canard  and  vane,  both  at 
maximum  effort  to  pitch  the  nose  up,  were  countered  by  the  wing's 
nose  down  pitch  moment.  The  high  nose  down  pitching  moment  was 
due  to  the  aircraft's  high  static  margin  coupled  with  the  canard 
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Attitude,  h,  ft 


design.  Another  characteristic  of  the  level  deceleration  profile  was 
the  large  downrange  distance  used  to  slow  the  vehicle. 


Downrange  distance,  x,  ft 


Figure  VII.18:Level  Deceleration  Attempt:  Situational  Plot 
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t  Figure  VII.  19 

Level  Deceleration  Attempt:  Balance  of  Moments 

c.  ''Zooming"  profile 

Attempts  made  using  the  second  type  of  profile  met 

j 

with  success.  Figure  VII.20  shows  a  transition  first  to  horizontal 
flight  and  then  back  to  a  hover.  The  transition  back  to  a  hover  began 
with  the  vehicle  at  approximately  70  ft/sec,  at  which  time  power 
was  added  and  the  nose  raised  with  canard  and  vane  as  the  vehicle 
climbed  at  angles  of  attack  less  than  stall.  As  the  pitch  attitude 

I 

I 

I 
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neared  90  degrees  power  was  reduced  and  an  opposite  vane/canard 
input  was  made  to  stop  the  pitch  rate.  Appendix  H,  Figure  H.6 
depicts  the  aircraft  response  and  control  inputs  made  during  this 
particular  run. 

The  altitude  required  to  perform  this  type  of  maneuver 
varied  from  350  to  900  feet,  depending  on  the  speed  at  which  the 
zoom  was  initiated  and  the  throttle  scheduling  used. 


Figure  VII. 20 

Transition  to  Horizontal  Flight 
Followed  by  "Zooming"  Transition  Back  to  Hover 
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VIII.  CONCLUSIONS 


A.  GENERAL 

The  design  and  initial  fabrication  of  a  Vertical  Attitude  Takeoff 
and  Landing  (VATOL)  Unmanned  Air  Vehicle  (UAV)  airframe  was 
completed  at  the  Naval  Postgraduate  School  (NPS).  The  vehicle, 
called  Archytas,  was  a  combination  of  two  existing  UAV’s—the  AROD 
and  Aquila— as  well  as  locally  manufactured  components,  including 
wing  spars  and  support  structure  for  a  canard.  The  objective  of 
creating  Archytas  was  to  provide  a  proof-of-concept  platform  for 
research  to  explore  performance  trade-offs  and  stability 
augmentation.  A  three-degree-of-freedom  simulation  was  used  as 
the  focus  of  the  design  efforts,  to  validate  design  decisions  made  in 
the  fields  of  propulsion,  aerodynamics,  structures  and  flight 
mechanics.  Engi..e  tests  were  conducted  to  determine  thrust  and 
control  power.  Structural  components  were  designed,  fabricated  and 
then  tested,  making  modifications  where  necessary  to  ensure 
sufficient  airframe  strength.  A  longitudinal  control  system  was 
designed,  validated  by  simulation  and  tested  structurally. 

Specific  conclusions,  by  section,  follow. 
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B.  SPECIFIC 

1 .  Engine  tests 

An  existing  engine  test  stand  was  upgraded  with  improved 
force  and  moment  measuring  devices.  The  existing  engine  from  the 
AROD  was  tested  and  shown  to  provide  sufficient  thrust  and 
longitudinal  control  power  to  power  the  Archytas  design. 
Specifically,  thrust  and  control  power  were  quantified  as  a  function 
of  rpm  and  throttle  position.  Additionally,  several  miscellaneous 
characteristics  were  explored  including  duct  blockage,  vane 
movement  in  the  presence  of  propwash,  the  speed  of  the  propwash, 
noise  levels,  fuel  usage  and  engine  dynamics. 

2.  Wing  spar  design,  fabrication  and  test 

The  structural  elements  to  connect  the  Aquila  wings  to  the 
AROD  duct,  including  wing  spars  and  associated  hardware,  were 
designed,  fabricated  and  tested  on  an  individual  basis.  These 
individual  tests  revealed  design  deficiencies  which  were  corrected 
and  the  components  then  retested  with  satisfactory  results.  A  full 
configuration  test  was  then  conducted  to  simulate  the  flight  loads  on 
an  actual  vehicle.  This  full  test  demonstrated  that  the  basic  spar 
design  was  sound  and  identified  deficiencies  with  the  Aquila  wing 
root  structure. 

3.  Canard  design;  support  structure  fabrication  and 
test 

A  static  longitudinal  stability  analysis  was  conducted  to 
choose  the  size  and  location  of  a  canard.  Detailed  weight  and  balance 
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analysis  as  well  as  aerodynamic  modeling  dictated  a  shift  in  the 
location  of  an  electronics  pod  from  an  aft  center  of  gravity  location  to 
a  forward  one.  Canard  location  and  sizing  analysis  considered 
trimmability  and  controllability  criteria  over  the  Archytas'  expected 
airspeed  range. 

Following  canard  design,  a  structural  analysis  of  a 
supporting  boom  structure  was  detailed.  Empirical  analysis  of 
graphite  booms  was  conducted  and  eventually  abandoned  in  favor  of 
an  aluminum  tube  boom  structure  with  graphite  composite  support. 
Aluminum  tubing  was  chosen  based  on  a  combined  strength  and 
weight  comparison.  Boom  support  structural  members,  boom 
support  hole  in  forward  spar  and  boom  base,  were  manufactured  and 
tested  with  favorable  results. 

4.  Simulation 

Beginning  with  conservation  of  linear  and  angular 
momentum  equations,  a  non-linear  three-degree-of-freedom 
simulation  program  was  created.  The  series  of  programs  that 
constituted  the  simulation  allowed  for  analysis  of  the  Archytas' 
longitudinal  flight  path  and  pitch  axis.  Simulation  was  used  for 
validation  and  analysis  of  a  wide  variety  of  issues  including:  wing 
location  and  sizing;  canard  location  and  sizing;  weight  and  balance 
and  static  margin  optimization;  vehicle  thrust-to-weight  ratio 
analysis;  and  transition  profile  prediction  and  analysis.  Additionally, 
dynamic  characteristics  for  horizontal  flight  were  predicted  and 
compared  with  theory  and  Aquila  characteristics. 
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C.  ARCHYTAS  CONFIGURATION 

The  airframe  that  results  from  the  design  work  of  this  thesis  is 
the  Archytas.  With  the  addition  of  an  avionics  suite  for  command 
and  control,  the  Archytas  will  be  a  Vertical  Attitude  Takeoff  and 
Landing  (VATOL)  Unmanned  Air  Vehicle  (UAV)  capable  of  hovering 
flight,  horizontal  flight  and  transition  to  and  from  horizontal  flight. 
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IX.  RECOMMENDATIONS 


A.  GENERAL 

The  Archytas  program,  comprised  of  students  and  professors 
from  several  departments  at  NPS,  has  so  far  been  marked  by  a 
concerted  effort  to  approach  testing  in  a  careful  manner,  often  called 
"build-up"  in  the  test  community.  This  is  evident  in  the  testing  of 
individual  components  and  use  of  low  risk  platforms,  including  a 
remotely  controlled  car  and  conventional  airplane,  for  validation  of 
systems.  The  writer  strongly  encourages  continuation  of  this 
laudable  approach,  particularly  as  the  program  enters  flight  test, 
with  all  its  potential  hazards. 

B.  SPECIFIC 

1.  Engine  tests 

(a)  The  combination  of  engine  tests  and  simulation 
indicated  a  satisfactory  thrust-to-weight  ratio  for  the  Archytas. 
However,  if  the  vehicle's  weight  rises  beyond  about  110  Ibf,  further 
engine  tests  may  be  required  to  determine  optimum  thrust 
conditions.  Specifically,  a  matrix  of  different  blade  angles,  rpms,  and 
carburetor  mixtures  could  be  performed  to  optimize  thrust. 

(b)  More  accurate  fuel  flow  measurements,  perhaps  by 
using  a  commercially  available  fuel  flow  sensor,  could  be  taken 
during  future  engine  tests.  Such  fuel  flow  data  would  be  useful  in 


predicting  vehicle  range  and  endurance. 

(c)  As  discussed  in  thesis  work  by  Lieutenant  Commander 
Moran,  several  electronics  failures  were  experienced  during  engine 
testing  [Ref.  16].  Future  electronics  manufactured  must  be  built  in  a 
hardy  fashion.  Redundancy  in  design  is  also  encouraged  where 
possible. 

2.  Wing  spar  design 

(a)  The  spar  design  manufactured  by  the  FRL  appears  to  be 
satisfactory.  Modifications  to  the  Aquila  wing  root  will  be  required 
before  a  flight  article  can  be  assembled.  These  modifications  might 
include:  adding  graphite  and  fiberglass  material  to  the  wing  root  spar 
to  increase  its  crush  strength,  thereby  allowing  the  mounting  bolts  to 
be  tightened  more  than  was  feasible  during  testing;  use  of  skin  as 
structural  members;  and  use  of  structural  adhesive  to  supplement 
the  present  system  of  bolts  and  plates. 

(b)  All  design  modifications  should  be  tested  to  determine 
actual  strength.  Such  tests  should  eventually  include  a  full 
configuration  test. 

(c)  The  Department  of  Aeronautics  and  Astronautics  at  NFS 
should  have  a  finite  element  analysis  program  capable  of  performing 
analysis  on  composite  material  intended  for  manufacture.  The 
Department  presently  has  a  program  called  PAL2  but  use  of  this 
copy  is  for  educational  purposes  and  is  specifically  restricted  from 
being  used  to  assist  in  the  design  of  any  components  to  be 
manufactured  [Ref.  40]. 
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3.  Canard  design 

A  canard  must  be  designed  and  manufactured  to  sizes 
specified  in  this  thesis.  In  view  of  the  fact  that  canards  provide  lift 
for  the  vehicle,  a  slightly  cambered  airfoil,  such  as  a  NACA  2412,  is 
recommended.  The  recommended  method  of  fabrication  is  hot¬ 
wiring  of  blue  foam,  then  covering  the  foam  with  a  layer  of  graphite 
weave  and  then  a  finish  layer  of  fiberglass  weave.  Structural  tests  of 
a  sample  should  be  conducted. 

4.  Simulation 

The  fidelity  of  the  simulation  program  may  be  improved  as 
further  information  about  the  Archytas  is  determined.  Appendix  E 
discusses  the  assumptions  made  in  formulation  of  the  simulation 
program  and  is  therefore  a  good  source  of  potential  improvements  to 
the  program.  High-angle-of-attack  aerodynamics  is  perhaps  the  area 
that  deserves  the  most  attention. 

Most  simulation  runs  were  made  using  a  thrust-to-weight 
ratio  of  1.2.  If  the  flight  article's  thrust-to-weight  is  less,  the 
simulation  should  be  repeated  to  determine  acceptable  profiles  at  a 
higher  weight. 

Simulation  results  indicate  several  recommended  Archytas 
flying  techniques: 

•  Takeoff  transition  should  be  performed  by  allowing  the 
vehicle  to  accelerate  nearly  straight  up,  allowing  speed  to  build  to 
approximately  20  ft/sec  and  then  gently  pushing  over  with  vane 
control.  As  the  pitch  attitude  reaches  45  deg  the  vanes  may  be 


relaxed  to  neutral,  and  the  transition  to  level  flight  completed  with 
canard  control. 

•  Short  period  characteristics  in  horizontal  flight  are 
predicted  to  have  moderate  damping  and  frequency.  Care  should  be 
taken  to  avoid  excitation  of  the  short  period,  particularly  with  inputs 
near  the  natural  frequency. 

•  Phugoid  characteristics  in  horizontal  flight  are  predicted 
to  be  benign  and  should  present  no  adverse  flying  qualities  to  the 
Archytas  pilot. 

•  Due  to  the  vehicle's  large  static  margin  and  its  canard 
design,  the  transition  from  horizontal  flight  to  hover  will  necessarily 
have  to  be  made  using  a  "zooming"  profile.  The  simulation  showed 
that  a  deceleration  in  level  flight  to  a  hover  will  not  be  possible.  The 
preferred  method  will  be  to  "zoom"  the  vehicle  from  horizontal  flight, 
as  shown  in  section  VII. 

5.  Weight  control 

Thrust-to-weight  ratio  in  any  VTOL  design  is  critical.  The 
weight  of  candidate  components,  electronic  or  otherwise,  must  be 
high  on  the  designer's  list  of  priorities  and  continually  stressed  to 
team  members.  A  minimum  thrust-to-weight  ratio  of  approximately 
1.09  was  identified  during  simulation  and,  as  a  result,  an  absolute 
maximum  weight  of  110  lb  is  recommended.  Accurate  records 
should  be  kept  of  the  weight  of  all  components  added  to  the 
Archytas  to  assist  in  a  weight-reduction  analysis,  in  the  event  the 
vehicle  becomes  overweight. 


6.  Lateral  and  directional  flight  mechanics 

This  thesis  essentially  ignored  the  lateral  and  directional 
axes  of  the  Archytas.  Great  care  is  recommended  when  considering 
the  areas  of  horizontal  flight  directional  stability  and  high-angle-of- 
attack  lateral/directional  stability  and  control.  Wind  tunnel  tests 
may  be  useful  in  determining  design  r-'.quirements  in  these  areas. 
Specifically,  a  study  should  be  conducted  to  determine  directional 
stability  in  the  Archytas'  current  configuration  and  then  assess  the 
need  for  addition  of  a  lightweight  vertical  surface.  Active  control  of 
sideslip  through  use  of  the  AROD's  vanes,  essentially  augmented 
directional  stability,  is  possible  but  probably  would  not  work  well  in 
the  event  of  engine  failure. 

7.  Flight  test 

The  author  has  some  experience  in  flight  test  and  so  a  few 
words  on  the  subject  are  appropriate.  Prior  to  approving  any  flight 
test,  including  tethered  hover,  a  thorough  test  plan  is  recommended. 
This  test  plan  should  include  a  strong  description  of  exactly  what 
tests  are  to  be  performed  and  exactly  what  data  are  to  be  taken. 
The  safety  of  test  personnel  and  bystanders  must  be  thoroughly 
considered  and  "what-if"  questions  asked  and  answered.  A  matrix  of 
events,  showing  the  order,  description  and  purpose  of  each  event  as 
well  as  data  required,  will  be  a  useful  tool.  Build-up  of  events  should 
be  considered  and  specifically  addressed,  in  other  words,  low-risk 
events  should  precede  high-risk  ones.  Checklists  for  operation  of  the 
vehicle  and  its  systems  are  encouraged  as  a  means  of  ensuring 
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standardization  and  continuity  in  the  project.  If  possible,  the  test 
plan  should  be  reviewed  by  someone  with  experience  in  UAV  testing 
and/or  flight  test.  The  test  plan  should  be  approved  by  competent 
authority  prior  to  the  commencement  of  any  test  plan  events  and 
any  deviations  from  the  test  plan  must  first  receive  approval  from 
the  same  authority— the  time  for  creative  thinking  is  as  you  write 
the  test  plan,  not  as  you  execute  it. 
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APPENDIX  A:  TABLES 
TABLE  A.I 

AROD  THRUST  STAND  TEST  MATRIX 


EVENT 

THRUST  (1) 

DATA  (2) 

REMARKS 

A 

TSG.  MSG.  etc 

CALIBRATE  THRUST 
STAND  WITH  KNOWN 
LOADS. 

B 

Test  Team.  Date, 
AROD  number 
and 

Configuration 
(ignition 
source,  etc, 
etc???) 

Prestart.  See  note  (A) 

1 

0-1 00%  by  10% 
increments 

REMOVED 

TS,  rpm,  T,  TSG, 
M,  MSG,  Pa,Ta, 
event  number 

Warm-up,  basic  cal. 

First  run  up  to  full 
power  to  define  100%, 
note  rpm.  and  TS, 
divide  into  increments. 
Post:  Plot  TS  and  rpm 
vs.  T  and  decide  if  any 
non-linearities  exist. 

If  linear,  use  similar 
TS’s  for  all  events.  (1) 

2 

0-100%  by  10% 

O/O/O/O 

as  above  +  VS, 

VP 

Vanes  on  warm-up  and 
cal. 

3 

100% 

Vary  all  in  same 
direction  and 
magnitude 

Same  as  event  2 

Attempt  to  find  max 
thrust  vane  position- 
swirl  effect.  The 
position  for  max 
thrust  will  be  called 
VPMT 

4 

100%  ON  LY 

2.3.4:  VPMT;  1: 

VPMT  to  VPMT+  max 
available  throw  in 
Positive  direction,  in 
approximately  5 
degree  increments  for 
each  TS.  Repeat  in 
NEGATIVE  direction. 

As  above  + 
distance  from 
center  of  vane 
to  x-axis  of 
AROD. 

To  quantify  Single  vane 
control  power 

NOTE:  discontinue 
when  vane  stalls  (i.e. 
thrust  falls  ofO 

159 


Appendix  A 


5 

100%  ONLY 

1/3:  VPMT;  2/4; 

VPMT  to  VPMT+  max 
available  throw  in 
Positive  direction,  in 
approximately  5 
degree  increments  for 
each  TS.  Repeat  in 
NEGATIVE  direction. 

M 

Double  vane  control 
power;  vanes  deflected 
so  as  to  both  cause 
rolling  moment  in  same 
direction. 

NOTE:  discontinue 
when  vane  stalls  (i.e. 
thrust  falls  off) 

6 

100% 

50% 

75%,  in  that 
order 

Repeat  event  5  using 
all  four  vanes.  AT 
THREE  DIFFERENT 
THRUST  SETTINGS. 

M 

Total  vane  control 
power. 

NOTE;  discontinue 
when  vane  stalls  (i.e. 
thrust  falls  off) 

7 

100%,  chop  to 
75% 

VPMT 

rpm  initial, 
rpm  final,  T 
initial  and 
final,  TS  initial 
and  final,  TIME 
between  rpm 
initial  and 
final. 

Attempt  to  characterize 
thrust  response  during 
engine  deceleration. 
Qualitatively 
determine  if  the  spool 
down  sounds  linear. 

8 

75%  max  rate 
power 

application  to 
100% 

VPMT 

u 

Attempt  to  characterize 
thrust  response  during 
engine  acceleration. 
Qualitatively 
determine  if  the  spool 
up  sounds  linear. 

NOHES;  (A)  Prior  to  start  up  check  all  AROD  and  thrust  stand  components  for 
security.  Ensure  fire  extinguisher  is  available.  Brief  test  team  on:  Noise,  FOD, 
and  prop  arc  hazards. 

(1)  Initially  computer  input  value  for  throttle  setting  (TS)  will  be 
used  to  increment  engine  power.  Following  initial  runs  a  decision  will  be 
made  concerning  linearity  of  TS  vs.  T  (Thrust).  If  feasible,  TS  will  be  used  for 
follow-on  events. 

(2)  Abbreviations: 

TS:  Computer  input  Throttle  Setting 

rpm:  Revolutions  per  minute 

T:  Thrust  reading  in  LB  read  directly  from  scale. 

TSG:  Thrust  Scale  geometry,  i.e.  Moment  arm,  slider  location 
sufficient  to  calculate  actual  thrust. 

M:  Rolling  moment  reading  in  lb  read  directly  from  Moment 
scale. 

MSG:  Moment  Scale  geometry,  i.e.  Moment  arm,  slider  location 
sufficient  to  calculate  actual  Moment. 

Pa:  Ambient  pressure 
Ta:  Ambient  temperature 

VS:  Vane  positions  for  each  vane  as  input  from  computer. 
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TABLE  A.2 


AROD  ENGINE  TESTS:  4/16/93 


f  5 


;  5  5  S  5  Z  C  S  5  C  tf  5  2  T  T  ? 


isccsegeeistsaS 

^::r«ss2«aSeS 
e  a  «>  <*  • 


8SSS88888S88288SS=tSSS3 


?  it8!8!BStHtlllllinilttllltttltlttttltn 

}!}!l}}ll}}|}ll}}}li}tflll}|}l}lillllltl} 


i 


MoauiooooeaHntoaaaaaaaaaaaaaaaooaaaKVte^Naope 

*|SS8§i|i88H5388888SS8S38S83i|eSiEsi:SgS§i 

<<«a<ooeeooooopopppepooopoppooooapaaoeeee 

zxszxvxsc)4)tf)«»itf«etgMVMKVKM«wv8«wtca9a9aia9aaaaa 

I 

i  aa8S888S888888888888888888888888888S88888 

iiiiiiiiiiiiiiiiiiiiiiiiiiHiiiiiiliiiHii 

3  ijiir . 

I 

111 

2*  ^pp^r<-f<i«^apKpp 

! 

I 


Kpppppp— e*p'^pp»*pPP*-waa22PPap»-»i 


^  ^1*1 

ill  l!!l 


liliililili  nl 


ill  jf]! 


o^wp^pPKPPP-'ejpaapKPPP^gpappjjppp-gpjpPKppp-wp 
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TABLE  A.3 

AROD  ENGINE  TESTS:  4/27/93 


n  OOOOOOOOOOOOeOMOiAOMOiAOI>^OiAMeiAO«AOl0KOM>«AOlAeo«OOMOiAeoOMOMP 

• 


C««^r^«-tAOiA04««>iAOo«A*<»0«AOlAOIA^<AOpp^PO«»AO»-««»eo»MPO*AOnO«etf>e 

^ 


z  000000004 


ow>m<siOYM*’4 


a 

a 


ooow>ootAo«Aooo«Ao<AOMoeoooo«^oooooeoooeoooeoo«AOMoooooooeo 

s  ss;;g?s:s;;;i;;s?;;ss7ss:!&%s;:s;:8ss8ssssss«i8s::s&&ss:::isss8 


sIsIHiHslalsssIsssssSSsSiSSsssHssliSli 


eooooooooeo 


S.  OOOOOOOOOOOOIOOIOOOOOOPOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

Q.  ooooeooooooooeoooooo  o.o  ooooooooeooeooooeooooooeoooooe 

I  8  8  8  8  8  8  8  8  g  S  8  ;  g  §  8  8  8  X  ;  1 1 1 1 1 1 1 1 1  g  ^  §  |  f  S  §  §  g  §  S  S  §  S  ^  ^  C  C  t 


I 


«o^‘•PO'-«««^«Ao^>o«kO^^•r^^tf>ors.o«»0^c«n•««>•r>.•g»o^<«r»^w>0^-ooo^v«r»Ytf)0 

nnr^nn^^^v^vo^^vwiMMiniAWYtAwtAiAooooooooooro^vroKf^Kr'^rvfor^'Ooooooo 
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TABLE  A.4 

COMMERCIALLY  AVAILABLE  ALUMINUM  TUBE  CONHGURATIONS 
AND  STRUCTURAL  ANALYSIS  SPREADSHEET 


Figure  B.2:  Archytas  construction  detail:  plate  on  inside  ot  du 
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'  rPOC.nftH  "CAMBER",  nODiriEO  FROM  PROGRAM  OF  SAME  NAME  AS  USED  IN  AE3501 
C  rilTS  VERSION  INCLUDES  THE  DERIVATIVE  OF  THE  AQUILA  UAV ' S  MEAN  CAMBER 
r  I  TNE.  NOTE:  THIS  IS  FOB  THE  UNTWISTED  AIRFOIL.  THE  ACTUAL  AIRFOIL 
'•  HAS  THREE  DEGREES  OF  GEOMETRIC  TWIST  AS  SPECIFIED  IN  AD-A068  (REPORT 
C  ON  AQUILA ) . 


(  '  Tin  R.B.  STOHEV  NPS  WRITTEN;  DEC  1992 


‘  IMIS  PROGRAM  ALLOWS  YOU  TO  INPUT  THE  CONTROL  POINTS  OF  THE  VORTEX  PANELS 

•  ANt)  FIND  THE  RESULTING  AOA'S  FOB  A  GIVEN  CAMBER  DISTRIBUTION.  THIS  FILE 

•  MAY  BE  FAStl.Y  MODIFIED  FOR  DIFFERENT  DISTRIBUTION  EQUATIONS. 


PROGRAM  CAMBEF, 

BEAL  Y(120),XCtl20),ZC( 1 20 ) . DZC (1 20 ) , ALP ( 120) ,X0( 120) ,XTO( 120) 
OPEN  (UNIT-2C,FILE='CAMBER.IN' , STATUS- ' UNKNOWN ’ ) 

OPEN  ( unit-2  I , FI LE- ' CAMBER .OUT'  .STATUS- ’ UNKNOWN'  ) 

PRINT  *. 'INPUT  THE  NUMBER  OF  CONTROL  POINTS  ON  PLANFORM:' 

PRINT  * , ' ( CHORDWISE  TIMES  SPANHISE  STATIONS)' 

READ  *,N 

PRINT  INPUT  THE  Y-COORDINATE  OF  THE  WING  ROOT' 

PRINT  *,'0.0  IF  PLANFORM  IS  WING  ONLY:' 

READ  * , YCORD 
no  I-l,N 

REAIM20,*)  XQd  )  ,  XTO(  I)  ,  Y(  1  ) 

IP(Y(I)  .LT.  YCORD)  THEN 
tlDAIA  -  I 
END  1  F 

PRINT  *,  Y( I ) , YCORO,NDATA,N 
END  DO 
rriPD  -  n.O 
NPAMEL  -  0 
HSTART  1 
DO  ) n  ? - 1 , NDATA 

CORD  =  CORD  +  2*ABS(XTQ( 1 )-XO(I ) ) 

PRINT*,  1  ,CORD, Y(  I  +  l ) , Y( 1  ) 

NPANEL  =  NPANEL  ♦  1 
IF( Y( 1+ 1 ) .EQ.  Y( I ) )  GOTO  10 
XLE-XO(NSTART)-( XTO(NSTART)-XO(NSTABT) )/2 . 0 
DO  J=NSTART,NSTART+NPANEL-1 

PRINT*,  J .NSTART. NPANEL, CORD 
XC(J)  -  ABS( (XTQ(J)-XLE))/CORD 
C  AQUILA  AIRFOIL*»**»*»*************** 

IF(XC(J)  .LE.  .3)  THEN 

DZC(J)-  1.1 372-32. 154* XC(J)+332. 91 *XC(J)**2 

*  1 44S . 0*XC( J )**3*2190.35*XC(J)**4 
ALP(J)-  ATAN(DZC(J)) 

ELSE 

r)ZC(J)-  0.25532-1. 6 307 8*XC(J)+2. 75505* XC(J)**2 

*  -1. 41616* XC(J)**3 
Al  P(  J  )-  ATAN(  DZC(  vl  )  ) 

END  I  F 


END  DO 
NSTART-. 7 
NPANEL  -■  0 
'"ORD  -  0.0 
10  CONTINUE 

DO  I  -  NDATA* 1 , N 

ALP( r )  -  0.0 
END  DO 

r  REMEMBER  THE  SLOPE  -  -AOA 

WPITF(21,20)  ( -ALP( I ) . I-l ,N) 
20  FORMAT! RFIO . 4  ) 

END 
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■Program  "xplot”' 


% 

% 

Bob  Stoney  3/24/93  Naval  Postgraduate  School 
^  A  program  to  inplement  a  survey  of  static  longitudinal  stability  using 
"s  formulas  contained  in  Anderson's  text,  specifically  eq'n  7.30  on  pg  386 
^  IMs  program  analyzes  a  CANAPD  design,  the  Archytas  UAV 
%  being  developed  at  NPGS.  llTe  basic  idea  is  to  plot  both  the  center 
%  of  gravity  and  aerodynamic  center  versus  canard  CHORD  (which  is  sinply 
%  the  canard  AREA  divided  by  it's  span.. .in  the  archytas  configuration 

the  span  is  considered  a  constraint  (to  fit  the  canard  between  the  "goalpost" 

%  canard  booms)  .  Where  the  lines  cross  is  zero  static  longitudinal  stability 
i  and  to  the  left  of  the  intersection  is  varying  degrees  of  static  margin. 

%  (items  mar)ced  with  "****"  are  those  that  will  be  frequently  changed  by  user) 

^  (items  mar)ced  with  are  those  that  are  a  vector) 

i  "FRL"  is  fuselage  reference  line,  the  leading  edge  of  the  mean  aero  chord 
^  (negative  is  forward) ,  which  is  9  inches  aft  of  tlie  top  of  the  ABOD  duct . 

%  For  the  center  of  gravity  portion,  the  following  items'  weights  are  accounted 

%  for:  The  ARCD,  the  wings  and  hardware,  the  electronics,  the  canard  bocm,  and 

%  the  canard  itself.  See  the  end  of  the  program  listing  for  a  definition  of  the 

^  weights. 

clear 

r-lq 

lold  off 


- Solicit  input  frcm  user  for  configuration: - 


.AAA AAAAAAAAAAAAAAAA AA AAAAAAAAAAAAAAAAAAAAAAAAA AAAAA AA 


dispC Input  config  (l=^aft  pod; 2=forward  pod)') 
CCNFIG=input ( '  '); 


aaaaaaaaaaaaaa 


aaaaaaaaaaaaaa 


VAAAAAAAAAAAAAAAAAA 


^ A AA  A  A  A A A A  A A A  AAA A A A A A  A A AAA  A A AAAAAAAAAAAAAA AAA A A A A A AAA  A A A A A A A  A A  A A  A  AAAAAAAAAAAAAA 

% - Set  vqD  vector  of  Canard  chords  to  calculate  Xcg  and  Xac: - 

^aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaAaaaaaaaaaaaaaaaaa 

n=100;  %  number  of  increments  of  canard  chord. 

cc=linspace(.01,1.5,n);  %  Canard  diord  in  feet  @@@@ 


A  A  A  A  A  A  A  A  A  A  A  A.  AAA  A  A  A  AAA  A  A  A  A  A  A  A  AAA  A  A  AAA  A  A  AAA  AAA  A  A  A  A  A  A  A 


SAAAAAAAAAA 


S  A  A  A  A  A  A 


—AERODYNAMIC  CENTER  PCKTION— 


^  A  A A A  A  A A A A A A A AAA A A A A A A A A A A A A A A A A A A A AAA A A AAA A A A AAA A A  A A A A  A A AAA A A AAA  A  AAAAAAAAAAAA^ 

%  Uses  equation  7.30  on  pg  386  of  Anderson's  "Intro  to  Flight" 
dispC Input  t)ie  "Canardlocation" — the  distance  in  inches  frcm  the  duct  to’) 
disp ( 'canard  c/4.  BEME^©ER:  forward  is  negative. .. .enter  negative  nunlDer') 
ansl=input ( '  ' )  ; 
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Cancirdlcx3tion=ansl/12;  %  Distance  from  duct  lip  to  canard  quarter  chord 
Xacc=Canardlocation-9/12;  %  x  position  of  canard's  aero  center  frcm  FRL  ****** 
Xacwt)=8/12;  %  x  position  of  wing/body  aero  center  (in  feet) 

bc=31.5/12;  %  Canard  span  in  feet  (limited  by  boons)  ****** 

ac2cl=5.3;  %  Canard's  2  dimensional  lift  curve  sl(^5e  (per  radian) 

e=0.75;  %  efficiency  factor 

ac3d=ac2d. / (l.+ac2d./ (pi.*bc./cc. *e) ) ;  %Canard's  3-D  lift  curve  slope  (/rad)  @@@ 

deda=0.0;  %  d (epsilon-canard) /d (alpha) - the  change  in  cfcwnwash  of  the 

%  canard  due  to  angle  of  attaclc  change.  Pef:  Roslcam  part  6  fig 
%  8.67  page  274. 

awb=3. 77075;  %  Wing/boc^  lift  curve  slope  (per  radian)  (fror.  Vortex  Lattice) 

Sc=bc*cc;  %  Canard  area  in  ^"^2  @(?@(?@0 

5=29.2;  %  l^n  wing  area  in  ^^^2  (Used  in  vortex  lattice  program) 

Cr=3.327;  %  Ftoot  chord  in  ft  (Used  in  vortex  lattice  program) 

Ct=1.917;  %  Tip  chord  in  ft  (Used  in  vortex  lattice  program) 

cbar=2/3*  (Cr+Ct-Cr*Ct/ (Cr+cLH ;  %  Wing  mean  aerodynamic  chord 

VH=Xacc.*Sc./ (dDar.*S) ;  %  000000 

Xacabai=Xacvto/dDar  +  VH. *ac3d./av4>;  %  eq’n  7.30  (here  hn  is  same  as  Xacabar) 


A  A  A  A  A  A  AAA  A  AAA  A  A  A  AAA  AAA  AAAAA^AA  A  AAA  AAA  A  A  AAA  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  AAAAAAAA 

S - CEMTEK  OF  GRAVLTy  PO^nON - 

^AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

===WEIGHTS:  (in  lb) - - - 

Wballast=input ('Enter  ballast  wei^t  in  lb  (located  at  bocm  tip) ' ) ; 

%  ITus  is  a  ballast  to  try  and  iiiprove  the  static  margin 
Warod=52;  %  APOD  weight  (tpdated:  2/7/93)  See  "equipment  list"  at  prog,  end 
wfuel=1.5*6.8;  %  Weic^it  of  fuel  (eq'n  shows  conversion  from  gals  to  Ibf) 

%  Electronics  pod  wei^t  is  target  weio^  given  to  electronics  design  group 
%  See  equipment  list  at  program  end  for  desicription  of  actual  APOD 
%  electronics  pod  weii^ts. 
if  0CNFIG=1, 

Wpoc:t=9.0;  %  Aft  mount  oonfiguration 

else 

Wpod=9.5;  %  forward  mount  oonfiguration  (slightly  heavier  due  to 

%  supports  and  aerodynamic  fairings) 

end  %  if 

Wwings=23.2;  %  weight  of  IVID  WINGS  see  equipment  list  (updated: 2/7/93  ) 
Xbocrmount=22/12;  %  Permament  length  of  bocm  embedded  in  wing 
%  root,  up  to  lip  of  duct. 

Bocmlength=XboctimoLint+abs(Canardlocation)  %  Total  canard  txxxn  in  ft 
NboaTi=2;  %  number  of  bocms 

Bocmden=0. 312;  %  Weight  of  iDOcm  p)er  linear  foot  frcm  Aluminum  tube  design 
Wbocms=Bocmlength*BocirKten*Nboam;  %  lotal  weight  of  all  bocms 

Canardden=0 . 4 ;  %  Canard  weight  pjer  square  foot  (Bluebird  wing  w/2  servos:  .68) 
%  (Don's  R/C  plane  w/no  servos:  0.58)  Too)c  75%  for  foarn/glass  construction. 
WcanarcJ=tx:*cc*Canardden;  %  Total  weight  of  the  canard  0@@@@@ 

Wservo=0.25;  %  Weight  of  servo  (typ)e:IC  servo  MS-747WB) 

>^1^7^11=  (Warod+Wfuel+Wpod+Wwings+Wbocms+Wservo+Wballast)  *ones  (cc)  +Wcanard;  %@@@0 
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%  ■fO'lENr  APMS;  (ft) .  Reference  is  leading  edge  of  MftC,  fwd  is  negative — ^ 

XcgARCD=2. 94/12;  %  For  arod. 

FuelconfigF=l;  %  l=extemal  fuel,  placed  between  struts  **** 

%  2=AR0D  fuel  storage  location 

if  Fuelconfig=l, 

XccfUEI;=>-l/12;  %  for  "external"  fuel — no  fuel  in  AROD  body 

else 

XogFUED=-6. 75/12;  %  for  use  of  AROD's  fuel  tanks 

end  %if 

%  Electronics  pod  c.g.  d^jends  on  configuration.  For  config  1  (aft  pod)  the  c.g. 

is  as  presently  configured  on  ver..cle  (^ril  1993) .  For  config  2  it  is 
%  asstmed  that  the  pod  c.g.  will  be  at  the  same  location  as  that  for  the  AROD. 
if  0CNFIG=1, 

XcgPOC>=21/12;  %  Aft  mount  configuration  in  ft 

elseif  CCNFIG=2; 

XogPOD=-26. 75/12;  %  forward  mount  configuration  in  ft 

end  %  if 

XcgwiNGS=5/12;  %  for  wings  (determined  by  laalancing  wing  w/ccrponents  2/7/93) 
XcgB0CMS=Xbocrrinount-Bocmlength/2;  %  for  bocms 
xc<^lANARD=Xacc*ones(cc)+cc/4;  %  for  canard  0@@@0 

XcgSEF!VO=XcgCANARD;  %  for  servo,  assumes  servo  and  canard  cg's  the  same  @@@(? 
Xballast=Xacc-cc/4;  %  Assumes  ballast  is  positioned  on  front  of  boon  @00 
^  ^  ^^CMENrS  (in  ft-lb)==  ■: 

Mballast=44ballast.*Xballast;  %  @@@@@ 

Nfar  od=War  od*XcgABOD ; 

Mfuel=Wf  uel  *XccfUEL; 

^^Xld=♦^bod*XcgPCD  ; 

M«.ngs=Vfc(ings*XcgS4INGS; 

^*xcIns=♦^bocII^s*XogB0CM5; 

Mcanard=Wcanard.  *XogCAl«VRD;  % 

Mservo=Wservo*XcgSERVD;  %  @@@@@00 

Mro^Al;=(^brod+Mfuel+^^pod+^Wings+M^octns) ‘ones (cc)+Mcanard+Mservo+MDal last;  %  @000 

%  -  —Total  CENTER  CF  GRAVITy — 

XcgTOT\L=MIUI7\L,/WI01A^  %  000000 

XcgK7rAIbar=XcgTOrAL./cbar;  %  A/C's  eg  location  normalized  by  wing  @@@@ 

% - FIND  liXATION  CF  15%,  10%  AH)  5%  STATIC  fftRGINS - 

^ A A A A A A A A A  A A A AAA A A A AAA A A AAA A A AAA A A AAAAAAAAA A A A A A A A A A A A  A A  A AAA A A A A A A A A A  A A  A  A  A A A  A A  A 

plotflag5=0;plotflagl0=0;plotflagl5=0;  %  a  plotting  flag,  see  plotting  section 
for  i=l:n,  %  Note:  second  nutber  must  match  linspace  #  points 
3i(i)=Xacabar(i)-XogTDIALbar(i);  %  Static  margin 
if  abs 04 (i) -0.15) <=0.001,  %  to  identify  15%  static  margin  '/alue 

SM15=SM(i);  %  The  15%  static  margin  values 

XacabeurlS^Xacabard) ;  %  " 

XogTCri7UJDarl5=XcgTOrAIbar(i);  %  " 
ccl5=cc(i)*12;  %  " 

plotflagl5=l;  %  to  enable  polyline  plot  later 
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acl5=ac3d(i) ;  %  Cl  alpha  at  15%  SM 

Scl5=Sc(l);  %  Canard  area  at  15%  SM 

Cancdl5=cc(i) ;  %  Canard  chord  at  15%  SM 

Wcanarcil6=Wcanard(i) ;  %the  WT  of  the  canard  when  SM=15%  (put  on  plot) 
WIorAL15=WCTAL(i) ;  %  the  total  weight  when  SM=15%  (put  on  plot) 
elseif  abs(SM(i)-0.1)<=0.001,  %  to  identify  10%  static  margin  value 

SM10=SM(i);  %  The  10%  static  margin  values 

XacabarlCNXacabar (i) ;  %  " 

XcgTCa7\IJbarl0=XcgTOTAlJ3^  %  '• 
ccl0=cc(i)*12;  %  " 

plotflaglO=l;  %  to  enable  polyline  plot  later 
acl0=ac3d{i) ;  %  Cl  alpha  at  10%  SM 

SclO=Sc(i);  %  Canard  area  at  10% 

WcanardlOWcanardd) ;  %the  WT  of  tfie  canard  when  SM=10%  (put  on  plot) 
wrCTALlOWOTAL(i) ;  %  the  total  weight  when  3^10%  (put  on  plot) 
elseif  abs(SM(i)-0.05)<=0.00l,  %  to  identify  5%  static  margin 

SM5=SM(i);  %  The  5%  static  itargin  values 


Xacabar5=Xacabar(i) ;  %  " 

XcgTOmrj3ar5=XcglX7I7\Lbar(^  %  ” 

cc5=cc(i)  *12;  %  "  (converted  to  inches,  li)ce  plot) 

plocflag5=l;  %  to  enable  polyline  plot  later 
ac5=ac3d(i);  %  Cl  alf^  at  5%  SM 

Sc5=Sc(i);  %  Canard  area  at  5%  SM 


end  %if 


erxi  '»for 


••AAAAAAA*" 


- COMPARE  CG  AND  AERO  CENTER  CN  A  PLOT - 


V  A  A  A  A  A  A  A  AAA  A  A  AAA  A  A  A  A  AAA  A  AAA^' 


saaaaaaa** 


^AAA  A  A/* 


SAAAAAAAAAAAAAAAAAA 


piot(cc*12,xacabar,  '-\cc*12,XcgKTrALfcar, ) 
grid; xlabel  ( '(Sanard  chord,  cc,  in'); 

ylabel ('Non-dimen  CG  and  Neutral  Point  locations,  h,  hn,  aft  of  f^C  l.e.') 
title ( 'Canard  Sizing  Based  on  Static  Longitudinal  Stability') 
if  plotflag5=l, 
liold  on 


X5=[cc5  cc5);  %  for  polyline  function 

Y5=[Xacabar5  XogTOTAIiaarS] ;  %  " 

polyline (X5,Y5, ' — ')  %  draw  5%  static  margin  line 

texty=  (XacabarS-XogTOTALbarS)  /2+XcgTarAIbar5; 
text (cc5, texty, '5%  34')  %  put  label  on  5%  SM  line 

end  %if 

if  plotflaglC)==l, 
iTold  cn 

X10=[ccl0  cclO];  %  for  polyline  function 
Y10=(Xacabarl0  XcgTOTAUarlOJ ;  %  " 

polyline (XIO, YlO, ' — ')  %  plot  10%  static  margin  line 

tey.tY=  (XacabarlO-XcgTOI'ALbarlO)  /2+XcgTOrAIiDarlO; 
text  (ccl0,texty, '10%  34')  %  put  label  cn  10%  SM  line 

eixi  %if 
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if  plotflaqi  S=l, 
hold  on 

X15=(ccl5  ccl5J;  %  for  polyline  function 
Y15=(Xacabarl5  XcgTOrALbarlS];  %  " 

polyline (XI 5,  yi5,  '  — ')  %  plot  15%  static  nargin  line 

texty=  (XacabarlS-XcgTOIMJDarlS)  /2+XcgrOTAIiDarl5; 
text (ccl5,texty, '15%  34')  %  put  label  on  15%  34  line 

end  %if 

i  Plotting  "niceties": 

text (cc(10)*12,Xacabar (10) , 'A/C  Neutral  Point') 
text(cc(10)*12,XagrarALbar(10), ‘A/C  Center  of  Gravity') 

Ltop=sprintf ('Ballast  =  %g' ,  Wballast) ;text (.45, . 9, Ltop, 'sc' ) 

Ll=sprintf ( 'Wing  area  =  %g' ,S) ;text ( .45,  .87,l1, 'sc' ) 

L2=sprintf ( 'Wing/body  a  =  %g',awb) ;text(.45, .84,L2, 'sc') 

L3=sprintf ('Can. ("from  duct)  =  %g',Canarcilocation*12) ;text(.45,  .81,L3,  'sc') 
L4=sprintf  CXacv*)  =  %g',Xacv*i*12);text(.45,  .78,L4, 'sc') 

L5=sprintf ('Canard  span  =  %g',bc*12);text(.45, .75,L5, 'sc') 

L6=sprintf ( ' XcgABOD  =  %g' ,  XcgABCD*12) ; text ( . 45, . 72, L6, ' sc ' ) 

L7=sprintf  CXccpCD  =  %g' ,Xo^>0D*12)  ;text(.45,  .69,L7, 'sc' ) 
if  plotflagl5=l, 

L8=sprintf  ('Xacabar  (.15SM)  =%g'  ,Xacabarl5)  ;text  ( .  45,  .  66, L8,  'sc' ) 
L9=sprintf  CCG  (in)  (.15SM)  =  %g',XcgTarALbarl5*cbar*12)  ; 

text  ( .  45,  .  63,  L9, '  sc ' ) 

elseif  plotflaglO=l, 

L8=sprintf  ('Xacabar  (.lOSM)  =  %g',XacabarlO)  ;text  ( . 45,  .  66,L8, 'sc' ) 
end 

E«:cp=sprintf  CConfig:  %g ', CONFIG) ; text (.7,  .90,Rtqp,  'sc') 

Rl=sprintf  ( 'Warod  =  %g',Warod)  ;text (.7,  .87,Rl, 'sc' ) 

B2=sprintf ( 'Wfuel  =  %g' ,Wfuel) ;text (.7, .84,R2, 'sc' ) 

R3=sprintf ( 'Wpod  =  %g',Wpod) ;text(.7, .81,R3, 'sc') 

R4=sprintf  ( 'V^ings  =  %g' , Icings) , -text  (.7,  .78,R4, 'sc' ) 

R5=sprintf  ( 'Wbocms  =  %g',Wbocms) ; text  (. 7,  .75, R5, 'sc' ) 
if  plotflagl5=l, 

R6=sprintf  ( 'Wcanard( .  15SM)  =  %g'  ,Wcancu:dl5)  ;text  ( . 7,  . 72, R6,  ' sc' ) 
R7=sprintf  ( 'Wtotal  ( .15SM)  =  %g'  ,wrorAL15)  ;text  ( .  7, .  69,  R7, '  sc' ) 

R8=sprintf ('cc(.1534)  =  %g',ccl5);text(.7, .66,R8, 'sc') 
elseif  plotflaglO=l, 

R6=sprintf  ('Wcanarcl(.134)  =  %g',WcanardlO);text(.7,  .72,R6, 'sc') 
R7=sprintf(' wtotal  (.ISM)  =  %g',WIUEAL10);text(.7.  .69,R7, 'sc') 
end 

if  plotflaglO=l, 

R9=sprintf ('cc(.lOSM)  =  %g',ccl0) ;text ( .7, .63,R9, 'sc' ) 
end 
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% - Print  oiJt  for  Prof  H:  for  5,10,15%  SM  gives  ac3ci,  Xacc,  Sc 

^  A  yy  ^  A  A  At  A  Ay>N  A/y /yA  A.yy  A  AA  A  ^yv  A  A  .iNyV/N/^  /y  A  A>\  A>/V^>\  AAA 

disp { ' Conf iguratioi ' ) 
if  CCNFIG=1, 

disp  ('Pod  Aft') 
el  seif  C0NFIG=2, 

disp ( ' Pod  Fwd,  gyros  aft ' ) 
elseif  0CNFIG=3, 

disp ('Pod  Fwd, gyros  fwd') 
end 

if  plotflagl5=l, 
dispCFor  15%  SM: ' ) 

dispC  ac3r  Canard  ("duct)  c  chord  (in)  Sc  (ft^2)') 

outl5=[acl5,C:anardlocation*12,ccl5,Scl5) ;disp(outl5) 
end 

if  plotflaglO=l, 
dispCFor  10%  ' ) 

dispC  ac3D  CanardC'duct)  c  chord  (in)  Sc  (ft''2)') 

out 10= [aclO, Canardlocation*12, cclO, SclO] ; disp (outlO) 
exl 

if  plotflag5=l, 
dispCFor  5%  SM: ' ) 

dispC  ac3D  CanardC'duct)  c  chord  (in)  Sc  {ft''!)') 

out5= [ac5, Canardlocation*12, cc5, ScS] ;disp (out5) 

end 

^A  A  A  A  A  AAAAA  AA  AAAAA  AA  A  A  AA  AA  A  AAAAAAA  AAAA^-'y^A.AAAAAAAAAAA  A  AA  AAA  A  A  AA  A  A  A  A  A  A  A  A  A  AAA  A  A 

% - EQUIPMENT  LIST  (DEFINES  WHAT  IS  INCLUDED  IN  WEIGHTS) 

^A AAA AAAAA A AA AAA AAAAAAAAAAA A AAAAA AAAAAAAAAAAAAAAAAAAAAAAAAA AA AAA A AA AAAA AAA AAA 

'?=  ABOD  Weight  includes  (2/7/93) : 

%  Airframe  "Andy",  in  NPGS  configuration  (i.e.  with  no  avionics  )oetween  struts) . 
%  •  Carplete  engine  and  prop  assembly  (recently  run)  with  intaJce 
%  filters  installed  on  carbs  and  inlet  guide  vanes  installed. 

%  •  4  vanes,  8  centerline  supports  (2nd  set  of  4  installed  by  LXIDR  Moran  to 
%  steady  aft  mounted  electronics  pod) 

%  •  One  servo  between  struts  (fuel  management  servo?) .  71/2"  aft  of  FPL 

%  •  4  vane  servos. 

%  •  4  landing  gear  legs  (in  AROD  configuration,  not  yet  lengthened  for  Archytas) 
%  and  ccnposite  landing  gear  "hoop". 

%  •  NO:  pod  (or  associated  electronics),  fuel,  primer,  roll  caqe,  or  forward 
'•  "bullet"  included. 

%  Weight:  50.5  lb  +  630g  (f'or  "lioop") 

%  CG  position  (from  FRL) ;  with  hoop:  11  15/16"  w/o  hoop:  1)  1/8" 
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%  POD  WEIGHT  includes  (4/4/93) : 

%  •  Frame  (that  holds  circuit  cards);  s)d.n  (4"  ring,  8"  ring,  8"  cap) 

^  •  Electronics:  3  ~7"  cards  (fiber  cxjtic  card  removed);  1  MHz  DC-DC  converter; 
%  6  blade  power  supplies  in  tray;  ccnplete  gyro/acclercmeter  pacicage 
%  •  Board  that  smooths  generator  output  (the  one  with  t)Te  big 
%  orange  capacitor) . 

%  •  For  forward  mounting  includes  bullet  support.s  (4),  front  cap  (with  fan  and 
%  other  electronics. 

%  WEIGOT:  Aft  mount:  10  lb  Forward  mount:  11.72  lb 
%  Local  (X:  is  5  inches  aft  of  forward  lip  cf  gyro  sectiion  s)cin. 

'i;  See  also:  Figure  entitled  "Electronics  pod  dimensions — Forward  pod  config" 

%  (A  MacDraw  document) 

%  WINGS  weight  includes  (2/7/93) : 

%  •  RH  Aquila  wing  w/o  harci^are,  RH  shoulder;  RH  tip 

i  •  Forward  test  spar,  with  4  plates  for  holding  end  and  2  plates  for 
*  attaching  to  shoulder's  spar,  plus  8  bolts  to  put  it  all  together. 

%  •  Aft  test  spar,  with  aft  grips  (grip  weight  estimated) . 

‘  TCTAL  weight  of  two  wings  with  associated  corponents:  23.2  lb 

. Program  "xplot"-  - 
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%================— ==—=Pn3gram  ■  incidence"  == 


%  Works  cnly  for  15%  Static  Margin . 

%  Must  first  run  "3<plot" 

%  #####  indicates  quantity  to  be  changed  ty  user 


\  yv  A.  yv  y\ /V /N /v /s  ^  ^  ^  ^ /s/v  /V A  yv  yy /V  y\  yvyv /s /s  >< 


A  /\/\^yN/vyvyv^  /y/\/\/y/\/vyyy 


%  First ,  identify'  if  the  15%  static  irargin  ccnditicn  vas  possible  between  chord 
%  lengths  of  0  and  18 " : 


/Vyy  ^  yy  >\  yy  yv  ^>S  yy /\  y 


V  ^  yv^  ^  ^  ^  ^  ^ y\ ^  ^  yy  yy /y  yy yy  yy  yy  yy  y* 


clg 

hold  off 

if  plotf]agl5==l, 

q^l5;  %  nuiber  of  incranrents  for  incidace  angle  ####### 

r=2;  %  incranent  step  size  ####### 

for  i=l:r;q, 

initia’  incidenoe=-2;  %  in  degrees  ####### 
ic=  ( initialincidaice+-i-l)*pi/180; 

incidencedeg(i)=ic*  180/pi;  %  Ihe  incidence  for  this  iteration,  in  deg 

%  Find  slcpe  of  On  vs.  alpha  curve: 

lc=Xacc- (XcgTaEMJ3arl5)*c±Br;  %  Nfcment  arm  frcm  total  og  to  canard 

%  aero  center  (ft) 

Scl5=ccl5/12*bc;  %  ccl5  is  the  canard  chord  length  @  15%  SM  (in  inches) 
%  be  is  the  canard  span  in  ft  (see  "arderscnlcngxplGt") 
%  Sc  is  canard  area  in  ft^2 
VH=lc*Sci5/  (c±ar*S)  ; 
ac3d=ac:2d/  (l+ac2d/  (pi*bc/  (cx:15/12)*e) )  ; 

hacv4?=0.24;  %  Ihe  aero  carter  of  wing/boeV  nen-dim  ty  char 

OTBlpha  ( i )  =avio*  ( (XcgTOBOtorlS-hacvi))  -VH*ac3d/awb)  ; 

QraevheO. 00761;  %  Fran  vortex  latticre 

QtO  ( i )  =Qnacvh>-VH*ac3d*  ic; 

nfc50;  %  ninber  of  increments  for  next  value: 
alphaABS=lin^ce(0, 15,m)*pi/180; 

Qrogd,  :  )=QTO(i)+QTalpha(i)*alphaABS; 

x(i)=alphaABS (20)*  180/pi;  %  for  labelling  the  incidence  angle 

y(i)=QTc:g(i,20) ;  %  for  labelling  the  incicience  angle 

end  %  for  blexrk 

ch^(XcgTOTAIhBrl5-hacvb) 
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PLOT  QtDg  VS.  ACA 


iV/VAA  AAA  AAAAAA  AAAA  AAAA  A  AAAAAA  A/\/NyVA  AA 


plot  (alphaAES*  180/pi, Qnog) 

grid; xlabeK 'Absolute  angle  of  attack,  at^iia  a,  deg'); 
ylabel  (' Pitching  nnrent  coefficient  abcut  c.g.,  Qitg'); 
titleCAQA  vs.  Qtog  for  varying  incidence  angles') 

%  Label  the  Qncg  vs.  ADA  lines  for  each  incidence: 
for  i=l:r:q, 

Label=^rintf  { ' ic  =  %g' ,incidenced0g(i) ), -text  (x(i) ,y(i) , Label) 
end 


%  Put  urportant  info  cn  the  plot; 

Ll=Sprintf  ( 'Config  =  %g' .OONFTG)  ;teBCt(  .7,  .9,L1, ‘sc’ ) 

[j2=sprintf( 'Ballast  =  %g‘ .Vtellast)  ;te9d:(  .7,  .87,L2, 'sc* ) 

L3=^rintf  { 'h-hacvto  =  %g',XogTDrAIiDarl5-hac\4));te9et(  .7,  .84,L3,  'sc' ) 
L4=^rintf  { 'Wtotal  =  %g'  .WTOI'ALIS)  ;text(  .7,  .81,L4,  'sc* ) 


L5=£printf  ( 'oc(  .15SM)  =  %g' ,ccl5)  ;text(  .7,  .78,L5, 'sc' ) 

L.6=^rintf{ ‘Can.  ('fran  duct)  =  %g'  ,Canardlocaticn*12)  ;text  ( .7,  .75,16,  'sc' ) 
L7=qprintf  COG  (in)  (.15SM)  =  %g ' ,  XogTOrAIbarl5*(±ar*12 ) ; 

text(.7, .72,L7, 'sc') 

%  Put  lines  on  the  plot  for  Cruise  (CR)  and  Power  /tproach  (PA)  angle  of  attac)c 
xlineCR=  [2 . 28  2.28];  %  2.28  degrees  is  the  Cruise  angle  of  attac)c 

xlinePA=[10.61  10.61];  %  10.61  degrees  is  the  PA  angle  of  attack 

scales=axis;  %  ■axis"  will  return  the  axes  presently  in  use 
yline=  [scales (3)  scales(4)];  %  Sare  for  both  CR  and  PA 
polyline (xlineCR,yline)  ;polyline(xlinePA,yline) ;  %  Draws  vertical  lines  at 

%  the  desired  CR  and  PA  ADA's 
text(xlinfiCR(l),yline(l),  'Desired  CR  ADA') 
text  (xlinePA(l)  ,yline(l) ,  'Desired  PA  ADA' ) 


else 

di^(  'Program  incidence  not  run  — 15%  S.M.  point  did  not  exist  for  this' ) 
disp(  'ccnfig' ) 

end  %  if  to  see  if  15%  static  margin  point  existed 


A  A A A A A AAA A A AAAA AAA AAAAAAA AA A AA 


AAAAAAAAAAAAAAA AAAAAA AAA AAA A A A 

END  FRODFAM  •INCIEEICE'' 


N  AAA  A  AAAA  A  A  AAA  A  A  A  AAAA  AAAA  AAA  AAAAAAAAAAAAAAA  A^i 


\AA AAA A AAAA A AAA 
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APPENDIX  E 

ASSUMPTIONS  USED  IN  SIMUIATION  PROGRAM 

A.  GENERAL 

This  appendix  contains  many  of  the  assumptions  made  during 
formulation  of  the  AWA-series  3  DOF  simulation  programs.  The 
format  of  this  appendix  is  informal,  stating  the  assumption,  it's 
location  in  the  program,  justifying  information  and,  in  some  cases, 
suggestions  for  further  refining  the  assumption.  The  order  of  the 
assumptions  roughly  follows  the  order  reached  in  the  program 
AWAl. 

B.  ASSUMPTIONS 

(1)  ASSUMPTION:  Wing  area=29.2  ft'^2 
LOCATION:  archytasS  function 

JUSTIFICATION:  From  vortex  lattice,  Aquila  documentation  and 
actual  measurement  of  the  wing. 

SUGGESTIONS  FOR  REFINEMENT:  The  duct  will  certainly  add  to 
the  effective  wing  area-ducts  are  well  known  for  their  capability  to 
provide  lift.  A  study  of  duct  aero  might  provide  a  better  estimate  for 

S. 
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(2)  ASSUMPTION:  Cdo=.06 
LOCATION:  archytas  function 
JUSTIFICATION:  Prof  Howard's  estimate. 

SUGGESTIONS  FOR  REFINEMENT:  Should  account  for  fan,  shroud, 
etc.  Wind  tunnel  test 

(3)  ASSUMPTION:  Iyy=7.4  slugs-ftA2 
LOCATION:  archytas  function 

JUSTIFICATION:  Aquila's  was  10.4  according  to  ref  8.  Reduced  it 
25%  for  lighter  archytas 

SUGGESTIONS  FOR  REFINEMENT:  calculate  using  swing/time 
period  method 

(4)  ASSUMPTION:  Cmq»-1.47 

LOCATION:archytas  function,  used  in  control  function  for  pitch 
damping 

JUSTIFICATION:  This  is  Aquilas  Cmq,  as  per  Aquila  documentation. 
SUGGESTIONS  FOR  REFINEMENT:  Vortex  lattice  program,  using 
canard  and  shroud. 

(5)  ASSUMPTION:  Canard  and  vane  deflections  are  step  inputs,  no 
oscillations,  no  time  constant. 

LOCATION:  user  input 

JUSTIFICATION:  probably  not  true.  Depends  on  servos. 
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SUGGESTIONS  FOR  REFINEMENT:  once  flight  quality  servos  are 
obtained  they  could  be  tested  under  load  to  determine  fidelity  and 
time  constant.  May  need  to  be  modeled  as  a  noise  vector  in  the  state 
equations.  Note:  LCDR  Moran  and  CAPT  Kuechenmeister  have 
modeled  the  servo  response  for  their  thesis  work. 

(6)  ASSUMPTION:  Sea  level,  standard  day  density 
LOCATION:  main  program 

JUSTIFICATION:  simple,  probably  good  for  Monterey  operations. 
SUGGESTIONS  FOR  REFINEMENT:  could  modify  program  for  temp 
affects 

(7)  ASSUMPTION:  Aircraft  coefficient  of  lift  beyond  stall  as  per 
profile  discussed  in  Section  VI. 

LOCATION:  main  program,  after  "Ydot  vector  initialization"  block 
and  again  in  iteration. 

JUSTIFICATION: 

SUGGESTIONS  FOR  REFINEMENT:  I  could  find  no  high  alpha 
information  during  my  literature  search  for  similar  low  sweep  angle 
configurations.  Wind  tunnel  tests  would  further  refine  this  or  a  more 
thorough  lit  search. 

(8)  ASSUMPTION:  Drag  model  at  high  alpha  is  an  assumed  profile. 
LOCATION:  main  program 
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JUSTIFICATION:  Drag  varying  as  the  square  of  alpha  matches  well 
with  the  pre-stall  region  and  seemed  reasonable  for  the  post  stall 
region  (i.e.  drag  kept  going  up,  quickly,  past  stall). 

SUGGESTIONS  FOR  REFINEMENT:  Wind  tunnel/lit  search/CFD. 
Account  for  entire  vehicle, 

(9)  ASSUMPTION:  Canard  pitch  authority  is  simply  based  on  a  Cl 
vs  alpha  (using  overall  vehicle  alpha  and  an  extrapolated  NACA  0009 
airfoil  times  dynamic  pressure  (same  as  vehicle's)  times  distance  to 
c.g. 

LOCATION:control  function 

JUSTIFICATION:  Good  first  approximation.  No  downwash  to  wing 
from  canard  is  probably  a  good  approximation  in  that  the  duct  may 
ingest  the  wake  of  the  canard, 

SUGGESTIONS  FOR  REFINEMENT:  Modify  airfoil  characteristics 
once  configuration  is  known.  Wind  tunnel  fore  i  measurement  for  Cl 
vs,  alpha  and  flow  visualization  for  canard-wing  downwash  (and 
flow  energization)  effects.  Lit  search? 

( 10)  ASSUMPTION:  No  contribution  to  pitching  moment  from  alpha 
dot  or  q  dot 

LOCATION:  control  function 

JUSTIFICATION:  neglected.  Is  often  neglected  in  conventionally 
configured  aircraft  where  the  effect  would  be  more  expected  (i.e. 
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time  lag  of  downwash  transport)  so  I  neglected  it  here,  feeling  that 
the  canard's  contribution  to  downwash  over  the  wing  would  be  small. 
SUGGESTIONS  FOR  REFINEMENT: 


(11)  ASSUMPTION:  Gyroscopic  effects  not  considered 
LOCATION: 

JUSTIFICATION: 

SUGGESTIONS  FOR  REFINEMENT:  Model  them. 
Kuechenmeister  has  done  this  in  his  thesis  work. 


CAPT 
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PROGRAM  AWAl 
LCDR  R.B.  Stoney 
Naval  Postgraduate  School 


Call  function  ArchytasS" 
to  retrieve  constants 


ArchytasS 


S.K.Cdo.lyy.Lc.Lv.Lw, 

Sc,Sv,cbar,Cmq 


Solicit  user  for: 
#  iteration 
DELTA! 

Weight 
Control  Mode 
Suspend? 


Default 

available 


SUSPEND? 

(for  user  to  make 
control  changes) 


'How  many  secs 
between  suspensions? 
Read:  "UPDATE”  ! 


UPDATE*0 


Control  Mode  ^ 
(options:  1  =fixed  q 
2=vane  &  canard  . 


Solicit: 

Canard  deflection 
Vane  deflection 
Throttle  setting 


Next  page 


Next  page 
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"Look  up" 

Drag  coefficient 
(if/then/else)  | 


WRITE  TO  OUTPUT  FILES: 

^OUTFILE(1  ,:)=[x,h,V,Vdot,THETA, ALPHA, counter  #,T,xdot,hdot 
gamma, Energy  height,  energy  height  dot,q,qdot,canard  defi, 
vane  defl,time] 

EXTRA(  1  ,:)=[counter#,gammadot,Q,row,CI,Cd,x/h/\//gamma 
at  half  step,Mcgw,Mcgc,Mcgv,Mq,L,terma,termb,Q,CI] 

**This  is  an  output  file  used  for  troubleshooting 


ST0PFLAG=1 
(this  is  a  counter  to  end  the  program  when 
the  use.r-inpute(i  #  iterations  is  reached) 


While  looD 

B 

"v  STOPFLAG»999 

A 

STOPFLAG~=999 


PT=PT+1  (counter) 


UPDATECOUNT=UPDATECOUNT=1  (suspend  option) 


T 


Next  page 
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CALL  FUNCTION  "YPRIMES' 


Current  Y(=[x,h,V, gamma] 
THETA, T.S.M,Q.CI.Cd.ALPHA 


Passed  back; 

Ydot(*[xdot,ydot. 

Vdot.gammadot] 


"yprimes 

equations 

VII.A.B,C,G 


Half  step  forward 
differencing  (equation 
VII.XX) 


As  before,  write  to  output  files: 
OUTFILE(PT,:) 
and 

EXTRA(PT,:) 


h<=0  ? 

(i.e.  the  vehicle 
hit  the  ground) 


PRINT:  "Ground  impact" 


Next  page ' 


While  loop 


STOPFLAG=999 
(will  kick  us  out 
of  while  loop) 
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i 


Pass  in:  q.DELTAT, THETA, T, ALPHA, {Y}, 
lyy,Q,S,Sc,Sv,cbar,Lc,Lv,Lw, row, Dele, De 


Call  function  "control  i" 

Iv.Cmo 

■I 

Function  "controls" 

Pass  back: 

(see  next  page) 

qdot,q,THETA 
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- PROGRAM  "AWAl" 


?  LiCDR  Bci)  Stoney,  Naval  Postgraduate  Schcol,  Monterey 
^  Written  for  APPLE  NRCINIOSH  ® 

i  Last  i.pdate:  12  MAY  93  (Throttle  nodel  and  vane  model  ipdated) 

29  Nfeiy  93  (PRCFILE  question  added  at  beginning  of  program 
’?  to  optimize  continuously  vpdated  plot  for  different 

flic^t  profiles  (hover  versus  horizontal  flight) 

Three  degree  of  freedon  (pitch,  x,  z)  simulation  for  ARCHYTAS  Vertical 
»  Attitude  Takeoff  and  landing  (VATOL)  UAV.  Includes  inprovesnents  over  PMS3, 

^  in  that  the  equations  of  motion  (contained  in  '*yprime4’’  and  "control 3"  do  not 
%  have  any  small  angle  assjiptions.  Also  incorporates  autoscaling  plot  feature 
%  for  contincusly  updated  dcwnrange  distance  vs.  height  plot. 

^  Thanks  to  Colin  Cooper  (aka  "Mr  MMIAB")  of  the  EC  department  for  help  with 
%  making  the  graphics  look  nice. 

Note:  •'**»****••  means  an  area  open  to  further  progrartming  (may  be 
%  accotpanied  by  artpllf’ying  informaticxi  concerning  work  to  be  done) 

clear  %  Clear  all  variables  from  system-wipe  the  slate  clean 

clc  %  Clear  the  camand  window 

clg  %  Clear  the  graph  window 

hold  off  %  To  release  any  previous  graph  settings 


%  PROGRAM  SCHEME 

%  The  algorithm  is  performed  as  follows:  (line  30) 

%  •  Aircraft  constants — gathers  needed  constants  for  particular  aircraft 

fron  a  function. 

%  •  Stcurting  values — initializes  model's  position,  orientation  and  velocity 

•  Iteration — Utilizes  predictor  and  corrector  to  increment  desired 

quantities . 

%  Predictor — first  guess  at  the  solution  of  the  governing  Dif  EQ's 

%  Corrector — Corrected  solutions  to  the  governing  Dif  BQ's 

•  Output  file  management — sets  ’jp  a  matrix  of  the  most  inportant 

%  quantities,  which  are  recorded  at  each  iteration. 

%  •  Aircraft  control — Uses  a  function  to  control  system  irputs. 

%  •  Final  output — To  let  the  user  extract  desired  information. 
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•  Fi-inctions  This  program  utilizes  the  following  created  functions; 
"archytas3" — aircraft  data  for  the  archytas  UAV 
"vprtmes4"— calculates  y'  sclutior."  to  co'jpled  system  equations 
’’ccntrol3" — provides  adrcraft  control  inputs 
"labeler" — piovi  Jes  labels  for  the  plots 

"THROTILE" — converts  user  input  of  throttle  position  (in  %)  to 

Thrust  (in  Ibf)  based  on  engine  test  data  (see  "IT  vs.  T") 

•  Other  programs  to  help  plot  inputs  and  response; 

"plotter" — ^plots  control  inputs  and  vehicle's  response 
"Nplot"  used  for  troubleshooting.  Shews  mcments  into  conservation 
of  an^^lar  mcmentum  equation  contained  in  "control 4" 

"Movie" — replays  the  vehicle's  position  and  pitef 

attitude  during  tiie  most  recent  run  or  use  "load  workspace" 


L1-1;L2=1;L3=1,-L4=1;L5=1;  %  Flags  for  upeteted  plotting  routine 

- - Aircraft  constants - 

( ? ,  K,  CDo,  lyy ,  Lc,  Lv,  Im,  Sc,  Sv,  Qaar ,  Qnq)  =arch^'tas3; 

.tee  Archytas  function  for  description  of  call  variables.  (lirte  60) 

For  plot  optimization,  as)c: - 

(.iispC  ') 

di?[M'What  tipe  of  maneuver?  l=takeoff-transit’on  or  tKver;  2=level  flight') 
d' ST  Cor  horizontal  flight  to  landing  transit : ) 

PR.)FlLE=input('  ' )  ; 

- Starting  values - 

dispC  ') 

dispCInput  as  a  rew  vector:  [iiteraticans, delta  t.  Weight  (Ibf ),')  ; 
disp( 'Control  mode(l=fixed  pitch  rate  or  2=cont  canard/vane) , and' ; ; 
disp( 'suspend  capability  (l=on,0=cff) ] .  Should  be  five  inputs.'); 
dispCOR  TYPE  1  FCB.  DEFAULT  VALUES=(500  0.05  100  2  1)’) 

INITTALeinput  ( '  ' ) ; 

if  INITIAL=1,  %  DEFAULT  VALUES: 

ITEBS=500;  %No.  of  iterations 

DELTAT=0.05;  %Time  increment  between  calculations 

wt=100;  %Aircraft  weight  in  Ibf 

CCNIMDDE>=2;  %Control  Mode  (l=const  q;  2=user  inputs  Del  c&v) 

SUSPEND=1;  %Suspend  flag  (to  view/control  du’-ing  run) 

else,  %  (User  ^jecified  values: ) 

I'ITRS=INITIAL(1) ;  %No.  of  iterations 

DELTAT=INITIAL(2) ;  %Time  increment  between  calculations 

Wt=INI’nAL  (3) ;  %Aircraf t  weight  in  Ibf 

CCN1M3DE=INITIAL ( 4 ) ;  ^Control  Mode  (1-con.stant  q;  Pmiser  inputs  I'^l  cSv) 
SUSPEND=INITIAL(5) ;  %Suspend  flag  (to  view/control  during  run) 

end 
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- Block  to  set  flag  ("UPDATE”)  to  alert  program  to  periodically 

' - halt  program  for  view/control. ..  .or  will  allow  uninternpted  run 

if  SUSPEND=1,  %  (line  90) 

dispC Input  the  #secs  (oetween  suspensions'); 

UPDATE=irput  ( '  ' ) ;  %  Time  period  Isetween  '/iew/crntrol 

elseif  SUSPEND==0, 

UPDATB=0; 

else 

dispCYou  did  not  enter  a  0  or  1  for  suspend  option — default  is  O') 
UPDATEH); 

end 

- End  UPDATE  set  loop - 


- Control  input  block - (line  101) 

if  C0NIM3DE=1,  %  Constant  pitch  rate  mode  (mostly  for  troubleshooting) 
ciispC  [Input  the  desired  pitch  rate  in  deg/ sec  (dcwn  is  neg)  and  TP]  ') 
qdeg=input ( '  ' ) ; 

q=qdeg(l)*pi/180;  %  Pitch  rate  converted  to  rad/ sec 

qdot=0;  %  Because  coristant  pitch  rate 

Delc“0;  %  Canard  angle  set  to  0  for  output  file 

Delv=0;  %  Vane  angle  set  to  0  for  output  file 

TP=qdeg(2);  %  Tbrottle  position  in  percent 
else, 

if  CC>mCDE~=2, 

dispCYour  input  for  control  mode  was  not  1  or  2.. default  is  2'); 
else 

dispC Input  as  row  vector  initial  canard  deflection  (LE  dcwn  is  +), ') 
dispCVane  deflection  (TE  ddwn  is  +)  and  Throttle  position  ') 
disp (' (percent — 0  is  idle  100  is  full  throttle):  [Dele,  Delv,TP]') 
dispCOR  TYPE  1  PC®  DEFAULT:  (-2  0  100]  ')  %  see  default  loop  below 

CONnNIT=input ( '  '); 
if  a:iNITNIT=l, 

Delc=^2*pi/180;  %  Default  canard  deflection  in  rad (line  120) 

Delv=0;  %  Default  vane  deflection  in  rad 

TP=100;  %  Default  throttle  position  in  percent 

else 

Delc=CCNriNIT(l)*pi/180;  %  Canard  deflection  converted  to  rad 
Delv=O0NnNIT(2)*pi/180;  %  Vane  deflection  converted  to  rad 
TPKXNITNITO) ;  %  Throttle  position  in  percent 

end  %OCNnNIT=l? 

cf=0;  %  Pitch  rate  set  to  0  for  first  entry  in  output  file 

qdot=0;  %  Pitch  accel  set  to  0  for  first  entry  in  output  file 

end  %C0NIMIS:-=2? 
end  %OCNIMDra:  if,  else 


[T|=THPaiTLE(TP) ;  %  Ccill  function  that  converts  Throttle  position  to  thrust 

^ - End  Control  input  block - 
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CF32.174; 

NW«)t/g;  IWass  in  slugs 

Pl-1;  %  Counter 

UPDATEC0UNT=1 ;  %  Counter  for  program’s  option  to  suspend  run  for  view/irput 

FiCW=0. 002 37688;  %  Density  in  slugs/ft''3 

- The  Y  vector  (corprised  of  4  elements  that  constitute  the  depenctent 

^variables  in  the  differential  equations)  is  as  follows: 

%  Y(l)=x  position  in  ft,  Y(2)=height  in  ft,  Y (3)=Velocity  ft/sec, 

»  Y(4)=qartfna  (flight  path  angle  in  radians).  C>r:  [Y]  =  [x,h,V,garrina] Transpose 
%  This  vector  will  be  overwritten  throughout  the  program,  with  required  values 
i  (line  145) 

^  being  retained  in  the  OOTFILE  matrix  (to  be  introduced  below) . 
dispC  ') 

dispCInput  as  row  vector  (deg)  (hi,Vi,garma(i) — deg,and  theta(i) — deg)  ') 
dispCNote:  Vi  must  not  be  zero  or  a  DIVIDE  BY  ZERO  error  will  occur.'; 
dispCCB  TYPE  1  FX»  DEFAULT:  [15.5  89.5  90]') 

INITIAL2=input ( '  ' ) ; 

Y(1)=0;  %lnitial  x  positicxi  in  ft 

if  INITIAL2=1,  %  i.e.  default  selected  by  user. 

INITIAL2(1)=15;  Y (2)=INITIAL2 (1) ;  %  These  are  all  spelled  out 

INITIAL2(2)=0.5;  Y (3)=INITIAL2 (2) ;  %  because  they  are  used  in  a 

INITIAL2 (3)=89. 5;  Y(4)=INI'nAL2(3)*pi/180;  %  configuration  definition  later 

INITIAL2(4)=90;  THETA=INITIAL2 (4) *pi/180; 

else 

Y(2)=INITIAL2(1);  %  Initial  hei^t  in  ft 

Y(3)=INITIAL2{2);  %  Initial  velocity  in  ft/sec 

if  INITIAL2(3)>89.5  &  INITIAL2{3)<90.5, 

INITIAL2(3)=89.5;  %  to  prevent  instability  in  nunerical  scher s  due  to 

%  singularity  in  gamna  dot  equation  contained  in  function  ";,primes" 
end 

Y(4)=INITIAL2 (3) *pi/180;  %  Initial  flight  path  angle  in  rad 

THETA=INITIAL2(4)*pi/180;  %  Initial  Pitch  attitude  in  rad 

end  %  Default  if  locp 

% - End  Y  vector  initialization - 

% - YDCTT  vector  initialization - 

% - Tlie  YDOT  vector  is  as  follows: 

%  YDOT(l)=Groundspeed  in  ft/s,  YD0T(2)=Rate  of  cliirti  in  ft/s, 

%  YDOT(3)=Acceleration  along  velocity  vector  in  ft/s^2 
%  YDar(4)=time  rate  of  change  of  fli^t  path  angle  in  radians/ sec 
%  Or;  [YDCT]  =  [Xdot,Hdot,VdDt,ganiiia  dot ] transpose 
This  vector  will  be  overwritten  throughout  the  program,  with  required  values 
toeing  retained  in  the  OUTFILE  matrix  (to  be  introduced  below) . 

YDOT  (1)=Y  (3)  *cos  (Y(4) ) ;  %Initial  groundspieed  in  ft/s 

YTOr(2)=Y (3) *sin(Y(4) ) ;  %  Initial  rate  of  climb  in  ft/s 

YDOT (3) =0.0;  %  Initial  acceleration  along  vel  vector  in  ft/s^2 
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'ifDOr(4)=O.0;  %  Initial  fit  path  angle  rate  in  rad/sec 

% - End  ydot  vector  initialization - 

%  Record  configuration: 

CONFIG=[Wt,  MXTAT,  INITIAL2  (1) ,  INITIAL2  (2) ,  INITIAL2  (3) ,  INITIAL2  (4) )  ; 

(>=0.5*B0W*y(3)^2;  %  Dynamic  pressure 

ALPHA=raETA-y(4);  %ACA  in  radians 

%  Fran  "Fincil  Archytas  Cl  vs.  ADA”  plot: 
if  (ALPHA<=*-1.57)  1  (AIPHA>=1 . 57) , 

CL=0; 

di^( 'Warning:  WING  aoa  greater  than  +/-  90  deg') 
el seif  ALPHA>-1.57  &  ALPHA<-0.72, 

CI/*-l .  083-0 . 69*ALPHA; 
el  seif  ALJ?WC>=-0.72  &  AIi>HA<^0.285, 

CD=2 . 81+33. 21  *AIPHA+105. 38*AU>HA^2+138 .  l*ALPHA^3+64 . 9*ALPHA''4; 
elseif  ALPWC>=-0.285  &  AIPfft<=0.29, 

CD=9 . 52e-2+3 . 38+ALPHA; 
elseif  ALPHA>0.29  &  ALPHA<0.9, 

CL=-2. 037+25. 54+ALPHA-72.15*ALPHA^2+82. 59*ALPHA'' 3-33. 4*ALPHA''4; 
elseif  AIPHA>=0.9  &  AIPHA<1.57, 
a>l .  875-1 , 194+ALPHA; 
end 

CD=0. 0608+0. 0543*ALPHA+0.910*ALPHA"2;  %  Coefficient  of  drag  as  fen  of  aoa 

%  for  post-stall  drag  increases  which  CD=f(CL^2)  doesn't  account  for. 

Es=Y(3)"2/(2*g)+Y(2);  %  Energy  height 

Esdot=Y(3)*YDOr(3)/g+YDar(2) ;  %  Time  rate  of  change  of  specific  energy — for 


% - BEGIN  FILE  ^®)^IAGE^EN^ - 

%  The  outfile  consists  of  the  following: 

%  OOTFILECrow#,  :)  =  [X,H,V, VcJot,lheta,Alpha,Counter#, Thrust, Groundspeed,Rate 

%  of  climb,  flight  path  angle,  energy  height,  time  rate  of  change 

%  of  Eng  Ht, pitch  rate,pitch  accel, canard  defl,vane  defl,time,TP] 

OUTFILE (1,  :)  =  (Y(l),Y(2),Y(3),YDar(3),THErA*180/pi,AIPHA*180/pi,PT,T,  .  .  . 

YDOTd)  ,YDCT(2),Y(4)*180/pi,Es,Esci3t,q*180/pi,q^t*180/pi,Delc*180/pi,  . . . 
Delv*180/pi,DEI,TAT*PT,TP] ;  %  Note  continuation  mar)cs  (...) 

EXTRA(PT,  :)  =  [PT,YDar(4),Q,RCW,CL,CD,0,0,0,0,0,0,0,0,0,0,0); 

%  Troubleshooting  info  (O's  for  sane  items  'cause  undefined  so  far) 
XHPlJ0rrCOUNrER=0;  %  Used  in  continuously  update  x-h  plot 
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ST0FFLAG=1 ; 

while  STCPFLAG~=999,  %  i.e.  terminates  icxp  when  STCPFLAG=999  (assigned) 

PT=PT+1;  %  Increment  counter 

UPDAIECXXJin'KJPDArECXJUMr+l;  %  Increment  "su^nd"  option  counter 
XHPIi7ItX)UNrai!=XHPI/?ItX^^  %  Used  in  continuously  update  x-h  plot 

ALPHA=THETA-Y(4);  %AQA  in  radians 

s  Predictor  and  intermediate  values 

[YDOr]=yprimes4 (Y, THETA, T,S,M,Q,CL,CD,ALPHA) ;  %  Call  to  the  D.Eqn  FXinction 
for  i=l:4, 

Y1 (i)=Y(i)+DELTAT*YDar(i) ;  %Euler  first  step  forward  predictor 

YHALF(i)=Y(i)+DELTAT/2»YDar{i);  %  Euler  half  Step  forward  predictor 

end  ^send  for 

%  Corrector  and  new  value  formulaticn 
[  YDOT]  =ypr  imes4  (YHALF ,  THETA,  T,  S,  M,  Q,  CL,  CD,  ALPHA)  ; 
for  i=l:4, 

Y2 (i ) =YHALF (i ) +DELIAT/2*YDOr (i) ;  %Richardson  extrapolation 

Y(i)=2*Y2(i)-Yl(i);  %Correct ion  made 

end  %end  for 

Es=Y(3)''2/(2*g)+Y(2);  %Energy  height 

Escbt=Y (3) *YDOr (3) /g+YDOT (2) ;  %New  TRC  of  energy  height 
!;>=0.5*RCW*Y(3)''2;  %  Dynamic  pressure 

ALPHA=THETA-Y(4);  %AQA  in  radians 

%  From  "Final  Archytas  Cl  vs.  AOA"  plot: 
if  (ALPHA<=«-1.57)  |  (AIPHA>=1 , 57) , 

CL=0;- 

disp( 'Warning:  WING  aoa  exceeds  +/-  90  deg') 
elseif  ALPHAO-1.57  &  AIPHA<-0.72, 

CI^l .  083-0 . 69*ALPHA; 
elseif  AIPHA>^0.72  &  AlJPHA<=-0 . 285, 

CL=2. 81+33. 21*AIi’HA+105. 38*ALPHA^2+138.1*ALPHA''3+64 . 9*ALPHA''4; 
elseif  ALPHA>=-0.285  &  ALPHA<=0.29, 

CL=9 . 52e-2+3 . 38*ALPHA; 
elseif  ALPHA>0.29  &  ALPHA<0.9, 

CL=-2. 037+25. 54*ALPHA-72.15*AIi>HA''2+82. 59*ALPHA''3-33. 4*ALPHA^4; 
elseif  ALPHA>=0.9  &  AIi’HA<1.57, 

CL=1 . 875-1 . 194+ALPHA; 
end 

CD=0. 0608+0. 0543*AIPHA+0.910*ALPHA''2;  %  Coefficient  of  drag  as  fen  of  aoa 

%  for  post-stall  dreig  increases  which  CD=f(CL"2)  doesn't  account  for. 
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% - Ojtput  file  update  (add  latest  iteration) - 

OUTFILE  (PT,  : )  =  (y  (1) ,  y  (2) ,  y  (3) ,  yOOTO) .  raETA*180/pi,  ALPHA*180/pi,  PT,  T _ 

yDar(l).yDOr(2).y(4)*180/pi,Es,Esdot,qri80/pi,qciot*180/pi,Delc*l80/pi. . . . 
Delv*180/pi  DELTAT*PT.TP] ; 

% - Aircraft  control - 

***Constant  q  or  Delc/Delv  control,  depending  on  CONTMODE  selection*** 
if  y(2)<=-0  %Stops  0  ground  impact 

disp( 'Ground  Inpact — and  that  is  bad') 

ST::FF1AG=999;  %Stops  while  locp 

end  %end  if 

[ THETA,  q,  q±3t,  Mcgw,  Mcgc,  Mcgv, Mq,  L,  Q,  CL)  =control4  (q,  DELTAT,  .  .  . 

THETA,  T,  ALPHA,  y,  lyy ,  Q,  S,  Sc,  Sv,  Cbar,  Lc,  Lv,  Lw,  RCW,  Dele,  Delv,  Qrq,  CCMIMCDE )  ; 
%Call  control  function 

if  PT=ITCRS,  %  i.e.  if  #  of  iterations  is  through 

STC®>FLAG=999;  %  Kicics  out  of  while  loop 

dispCAll  iterations  corrpleted') 
end  %end  if 

EX1KA(PT,  :)  =  [PT,yDar(4),Q,RCW,CL,CD,yHALF(l),yHALF(2),yHALF<3), . 

yHALF  ( 4 ) ,  Mcgw,  Mcgc,  Mcgv,  Mq,  L,  Q,  CL)  ; 

%  Extra  info  for  troubleshooting 


%x>ocooooooooooooooooooo<xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx)ooooooooooooooooo<xxxxx 

% - — Contin\x>usly  vpdated  plot  for  user  reference - 

%xxxxxxxxxxxxxxxxxxxxxxjooooooooooooo<xxxxxxx>OQQOOQQooooooQO<xxx>oooooooooooooo<xx 

TINE=.05;  %  Put  in  here  so  that  the  plot  wouldn’t  update  every  time  step... slow 

% - 

% - PKFILE  1  (takeoff  trans)  SECTION  CCNTINUOUS  PLOT - 


if  PPOn:LE=l,  %  i.e.  takeoff  transition 


%  Ibis  plot  places  a  circle  ("polymark")  at  the  vehicle’s  x/h  position  ^d 
%  draws  a  line  ("polyline")  at  an  angle  corresponding  to  a/c’s  pitch  attitu^. 
if  (XHPl£rrCOUNreR*DELTAT)>=TIPE,  %  i.e.  causes  this  plot  to  update 

if  ((y(l)<l00  &  y(2)<100)  &  L1=1),  %  TMs  if  locp  sets  plot  axes  for  ipdated 

hold  off 

V=f-25  100,0,100) ; axis  (V) ;  %  Axis  limits  (xmin,  xmax,  ymin,  ymax) 


plot (OOTFILE ( : , 1) ,COTFILE ( : , 2) , 'o’ ) 

xlabel (’Downrange  distance,  x,  ft’ ) ;ylabel (’Height,  y,  ft') 

L1=0;  %  set  the  plot  flag  to  zero  to  prevent  redrawing  of  plot 


elseif  (((y(l)>=100  s  y(l)<200)  t  (y(2)>=100  &  y(2)<200))  &  L2— 1), 

hold  off  . 

V=  [-25, 200,0, 200) ; axis (V) ;  %  Axis  limits  (xmin,  xmax,  ymin,  yrrwaX) 

plot (OUTFII£  ( : , 1) ,CUTFILE ( ; , 2) , ’o’ ) 

xlabel (’ Downrange  distance,  x,  ft ’);ylabel (’Height,  y,  ft') 


L2=0  ; 
hold  on 
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elseif  (((Y(l)>=200  S  Y(l)<400)  I  (y{2)>=200  &  Y(2)<400))  S  (L3=l)), 
hold  off 

V=  [-25, 400,0, 400] ;  axis(V);  %  Axis  limits  (xmin,  xiibx,  ymin,  ytnax) 
plot  (0(jrFII£(:,l),0UTFILE(:,2).  'oM 

xlabel (' Downrange  distance,  x,  ft' );y label ('Height,  y,  ft') 

L3=0; 
hold  on 

elseif  ( ( (y(l)>=400)  I  (y(2)>=400))  &  (L4=l)), 
hold  off 

v= [-25, 1000,0,1000] ;  axis(V);  %  Axis  limits  (xmin,  xmax,  ymin,  ymax) 
plot  (OUTFILE  ( : ,  1) , OUTFUE  ( : ,  2) .  'o' ) 

xlabel  (' Downrange  distance,  x,  ft' )  ;ylabel  ( 'Heicflit,  y,  ft') 

L4=0; 
hold  on 

elseif  (((Yd)  >=1000)  I  (Y  (2)>=1000) )  &  (L5=l)), 
hold  off 

V= [0,2000,0, 2000] ;  axis(V);  %  Axis  limits  (xmin,  xmax,  ymin,  ymax) 
plot  (OLTTFIIE  ( : ,  1 ) ,  OOTFILE  ( : ,  2 ) , ' O ' ) 

xlabel ('Downrange  distance,  x,  ft ');ylabel( 'Height,  y,  ft') 

L5=0; 
hold  on 

end  %  end  if 

%plot  (OmriLE  { : ,  1 ) , (XTIFILE  ( : ,  2) ) 
polymark  (OUTFILE  (PT,  1) , OUTFILE  (PT,  2) o' ) 

dv=(  (V(2)-V(l) )  /50)  *oos  (THETA)  ;dy=(  (V(4)-V(3) ) /50)  *sin(THETA)  ; 
polyline  ( [OITFILE  (PT,  1)  -dx  OITTFILE  (PT,  1)  +dx] ,  [CUTFIIE  (PT,  2)  -cfy  -  •  • 

OUTFILE (PT, 2) +<:^],  '-w')  %  This  shows  tlie  vehicle's  pitch  attitude 

XHPLanj00NrER=0;  %  To  reset 

end  %  end  if  for  '' (XHPIiOTaxjNTEK*DELTAT)  >=0 . 1'' 

% - PKFILE  2  (horiz.  or  landing  transition)  SECTICN  OF  CX:*TriNUOUS  PLOT 

elseif  PPOTI1£=2,  %  i.e.  horizontal  flight  optimized  display 
if  (XHPLCTCOUNrER*DELTAT)>=<rD«:,  %  i.e.  causes  this  plot  to  update  every  TIME 
if  Y(3)>=40,  %  Velocity  >  40?  If  so,  leaves  "big  picture  plot",  otherwise 

%  focuses  on  vehicle  for  control  to  hover .... 
if  ((Y(l)<1000  &  Y  (2) <1000)  &  Ll=l) 
hold  off 

V=[-25, 1000,0, 1000] ;axis(V);  %  Axis  limits  (xmin,  xmax,  ymin,  ymax) 
plot (CUTFILE ( : , 1 ) ,  OJTFILE ( : ,  2 ) , ' o ' ) 

xlabel ('Downrange  distance,  x,  ft' );ylabel(' Height,  y,  ft') 

L1=0;  %  set  the  plot  flag  to  zero  to  prevent  redrawing  of  plot 


)x)ld  on 
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elseif  (((Y(l)>=1000  &  Y(l)<2000)  I  (Y(2)>=1000  &  Y(2)<2000))  &  L2=l), 
hold  off 

V=[-25, 2000, 0,2000] , 'axis (V) ;  %  Axis  limits  (xmin,  xmax,  ymin,  ymax) 

plot (arrFIIZ(:,l),CHnFILE(:,2),  'o') 

xlabeK'Dcwnrange  distance,  x,  ft' ) ;ylabel  ('Height,  y,  ft') 

L2=0  ; 
hold  on 

elseif  (((Y(l)>=2000  &  Y(l)<4000)  I  (Y(2)>=2000  &  Y(2)<4000))  &  (L3=l)), 
hold  off 

(-25,  4000, 0, 4000] ;  axis(V);  %  Axis  limits  (Mtiin,  xmax,  ymin,  ymax) 
plot  (ajrFILE(:,l),CIUTFlLE(:,2),  'o') 

xlabeK'Dcwnrange  distance,  x,  ft' ) ;ylabel ('Height,  y,  ft') 

L>--0; 
hold  on 

end  %for  this  if;  "if  ((Y (1X1000  fi  Y(2)<1000)  &  U=l)" 
elseif  (Y(3)<40)  S  (L4— 1), 
hold  off 

V=[Y(1)-100,Y(1)+100,0,1.5*Y(2)];  axis(V);  %  x  +/-  100  and  0  to  l.Sheight 
plot (arTFIIZ  ( : , 1) ,OUTFILE ( : , 2) , 'o' ) 

xlaloel ( 'Dcwnrange  distance,  x,  ft ') ;ylabel ( 'Height,  y,  ft') 

L4=0; 
hold  on 

end  %  end  if  for  Airspeed  >  40?  (has  elseif  too) 
polymar)t  (OOTFIIE  (PT,  1)  ,OUTFILE  (PT,2) ,  'o' ) 
dx= ( (V (2 ) -V (1 ) ) /50) *oos (THETA) ;cfy= ( (V (4) -V ( 3) ) /50) *sin (THETA) ; 
polyline  ([afrFILE(PT,l)-dx  OUTFILE(PT,l)+dx] ,  [0UTFILE(PT,2)-dy  . . . 

0UmLE(Pr,2)+cfy], '-w')  %  This  shews  the  vehicle’s  pitch  attitude 

XHPIi3rixxUNreiM);  %  To  reset  (will  allow  plot  to  wor)c  in  another  .1  sec 
end  %  end  if  for  ''(XHPIOimJNlER*DEl,TAT)>=<).l'' 
end  %  if  for  which  PRCffTIJ: 


%XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXmXXXXXXXXX)OOOQOOOQOOOOOOOOaXXXXX 

% - END  Continuously  updated  plot  for  user  reference - 

%XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX5O0OO0O00QOOOOOQO00OQQ0<XXXXXXXXXXXXXXXXXXXXXXX 


% - Suspend  loop  (to  view  progress /update  control) - 

if  UPDATEOOUNr*DELTAT>=UPDATE,  %  i.e.  "Time  to  suspend?" 
UPDATEOOUNTK);  %reset  it  for  next  suspend 

ans=irput (' Program  suspended,  want  a  plot?  l=yes  0=no'); 
while  ans~=K), 

hold  off; 

axis([l  2  3  4]);axis; 

U=1 ;  L2=l ;  L3=l ;  L4=l  ; 
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di^( 'These  are  the  quantities  stored:') 


dispCl: 
disp('2: 
dispC  3; 
disp('4:  Accel 
dispCS:  Theta 


X  posit 

height 

^peed 


11: 

12: 


Fit  Path  Angle 
Energy  Ht  (Es) 
13:  Esdot 
14:  q 

qdot  ') 


16: 

17: 

18: 

19: 


Dele' ) 

Delv' ) 
time' ) 

Thrtl  Posit' ) 


15: 


6:  Alpha 
7:  Point# 

8:  Thrust 
9:  Groiindspd 
10: Rate  of  Cl 
%  Aslc  user  for  desired  x  and  y  plot  values: 
xanswer=irput { ' Type  #for  what  you  want  on  x-axis'); 
yanswer=input (' Type  #for  what  you  want  on  y-axis'); 
xanswer<l  |  xanswer>19  I  yanswer<l  I  yanswer>19,  %  (i.e.  mi sentry) 
xanswer=18;yaTswer=18;  %  just  plots  time  vs.  time  as  error  indication 


end 


plot  (OinriLEf :  .xanswer)  .O/TFILEf:  .vanswert ) 

[XIABEL,  YIAEEL]  =labeler  (xanswer,  yanswer )  ; 
xlabel  (XIABEL)  ;yla)3el  (YIABEL)  ; 
title  (['Vfeight:  ' ,  nuin2str  (Wt) , 'Ibf '  ] ) 
pause 

%  As)c  user  if  they  want  another  plot: 
ans=irput( 'Another  plot?  0=40' ); 
end  %  end  suspend  while  loop 
if  OCNIM3DE=2, 

MAT2=(Delc*180/pi,Delv*180/pi,TP] ; 

disp( 'Present  values  for  Delc,Delv  (in  deg)  &  TP:'); 

disp(^«^T2) 

ans2=input('WaiTt  to  change  these  settings?  04IO'); 
if  ans2~^, 

NlSJiCCNI^ir^xjt  ( 'Enter  as  row  matrix,  new  [Dele  (deg) ,  Delv  (deg) ,  TP  (%)]')  ; 
NEWOONTsi2e=size  (NEWOONT) ; 
chec)c=4EWOONrsi2e  (2) ; 
if  c±iec)c--=3, 

di^CYou  misentered  the  new  (Control — default  is  previoiis  input') 
else, 

Delc=f)EWXNT  (1)  *pi/180; 

Delv=NEWOCNr  (2)  *pi/180; 

TP=NEWC0Nr(3); 


end  %  this  if  block  prevents  program  terminaticxi  by  user  misentry. 
end  %ans  2  if 
elseif  CCNIMXE=1, 

dispC Present  value  for  pitch  rate  (deg/ sec)  &  TP:') 
disp(q*180/pi)  ;di^(T)  ; 

ans3=input( 'Would  you  like  to  change?  0=tio'); 
if  ans3~=0, 

NEWQT=input ( ' Irput  new  q,  in  deg/sec  (down  is  neg)  &  TP  as  [  ] ') 

cf4EWOr(l)*pi/180; 

TP=4®«3T(2); 
end  %  ans3  if 
end  %  0CNIMX>E=2  or  1? 

[T]=THRC)rTLE(TP);  %  call  to  function  "THPOTTLE" 

end  %  end  new  control  irput  if  loop 
end  %  end  suspend  if  loop 

- End  SUSPEM5  loop  block - 
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end  %end  the  iteration  hhlle  loop 
%  Di^lay  configuraticxi  to  user: 
disp ( ' Configuration : ' ) 
foimat  ccmpact 

di^C  wei^t  Delta  T  hi  Vi  Gannai  Ihetai  ') 

disp(OCl«’IG) 

format  short 

% - Final  Output - - - 

ansl=input  ( 'Want  a  plot?  l=yes  O=no');  %Kic)t  out  flag 

while  ansl~*0i 
hold  off; 

axis([l  2  3  4]);axis; 

L1=1;L2=1;L>1;I/4=1; 

disp ('These  are  the  quantities  stored:') 

dispel:  X  posit  6:  Alpha  11:  Fit  Path  Angle  16:  Dele') 

disp{'2:  heiqjit  7:  Point#  12:  Energy  Ht  (Es)  17:  DelV) 

di^('3:  ^3eed  8:  Thrust  13:  Esdot  18:  time') 

di^('4;  Aooel  9:  Groundspd  14:  q  19:  Thrtl  Posit') 

dispCS:  Theta  10:Rate  of  Cl  15:  qdbt  ') 

%  As)c  user  for  desired  x  and  y  plot  values 
xanswer=irputeiVpe  if  or  what  you  want  on  x-axis'); 
yanswer=irput ( 'Type  ifor  what  you  want  on  y-axis'); 
disp ('Note:  will  have  to  add  own  axis  labels,  etc’)  %******at  this  point 
plot  (0(TIFILE(:,xanswer),anFn£(:,yanswer) ) 

(XLABEL.YLABEL] -labeler  (xanswer,  yanswer); 
xlabel  (XLABEL)  ;ylabel  (VIABEL) ; 
title  (('Weight:  ',niin2str(Wt), 'Ibf  ]) 
pause 

%  As)c  user  if  they  want  another  plot: 
ansl=irput  ( 'Another  plot?  0=#O' ) ; 
end  %  this  vMle  loop 

% - ^Final  instructions  (help  for  storing  and  retrieving  data) - 

di^CRecomended  oenmands  to  save  data:') 
disp  ( '  »f ilename'OUTFILE; ' ) 

di^('»filenameC=OCt®TG;  (to  keep  a  record  of  the  configuration)') 

disp('»save  filename  filename;') 

disp('»save  filenameC  filenameC;') 

disp ('Then,  when  you  want  to  xise  it  again:') 

disp('»load  filename  or  filenameC  (filename  is  a  18  by  iiters  matrix  for  a 
particulcir  set  of  conditions — specified  on  filenameC') 

%  NOTE:  also  useful  are  programs  "filesaver"  and  "resurrect" 

% - NJTE:  See  adso  "M4A  assuiptions"  (Microsoft  word  docunent) 
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'Function  archytas 


f unction [S,  K, CDo,  lyy ,  Lc, Lv,  iM,  Sc,  Sv, Ooar, Qmq, xcg]  =archytas3  ( )  ; 
yvircraft  data  function  ****iir**Note  these  are  all  WAGS  at  this  point*** 

^  Used  by  program  PMS3 
%  LCDR  Bob  Stoney 
%  Last  ipdate:  8  May  93 
S=29.2;  %  Wing  Area  in  ft^2 

e=0.8;  %  Efficiency  factor 

AFM.25;  %  Aspect  ratio  (vortex  lattice  output) 

K=l/ (pi*AR*e) ;  %  Constant  for  CD  calculation  (CI>=CDo+K(CL''2) ) 

CDo=0.06;  %  Rick's  guess 

Iyy=7.4;  %  Moment  of  inertia  (slugs-ft''2)  Aquila=10.4  rept  AD-A068  345 

Lc=4.66;  Is  Absolute  Distance  in  ft  frcm  canard's  aero  center  to  a/c  c.g. 
L\*=1.5;  Distance  in  ft  fron  Vane's  aero  center  to  a/c  c.g.  ***Refine 
Lw=0.758;  %  Distance  in  ft  frcm  wing's  aero  ctr  to  a/c  c.g.  (Frcm  program 

%  "5^1ot") 

i^******************QjgQ^  MM  •  ****31. 5"x12"s=2 . 6  ft''2************************* 
Sc=2.6;  %  Area  of  canard  in  ft^2  (From  program  "Xplot") 

Sv^l.lll;  %  Area  of  two  (2!!)  vanes  in  ft^2  (Each  is  8"  by  10") 

CR=3.33;  %  Root  Chord  in  ft 

Ct=1.92;  %  Tip  chord  in  ft 

lambda=Ct/CR;  %  Taper  ratio 

Cbar=2/3*CR* (1  +  lanibda^2/ (1+lambda) ) ;  %  MAC 

Qik^-1.47;  %  Pitch  darrping  (1/rad)  Ref:  Aquila  rept  AD-A068  345***** 
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Ion  "IHCnUE' 


%  A  ccnpanlon  progrem  to  "ANAL"/  the  3  DCF  simulator. 

%  This  function  receives  the  user  Input  throttle  position  (TP)  (in  % — CW.dle, 

%  100=full  throttle)  and  converts  It  into  thrust  (T)  In  Ibf,  based  on  data  from 
%  engine  tests  on  4/16/93  and  shewn  curve  fit  on  the  crldcet  grs^^  file 
%  "Ihrust  vs.  IP" 

%  Bob  Stoney,  ^lav^Ll  Postgraduate  School,  created:  5/12/93 
%  Last  i;pdate:  5/12/93 

function (Tl=THBani£  (IP) ; 

flacpl; 

while  flag==l, 

if  1P>100  I  IP<0/ 

dispCTHE  raBOnUE  POSmON  YOU  INPUT  WAS  <0  CR  >100,  PE-INPUT  0-100') 

TP»“irput  ( '  ' ) ; 

flag=»l; 

else 

T=21 . 215- .  40269*IP+ .  10932*lP^2-2. 1907e-3*TP^3+l .  7515e-5*TP"4-5 . 1479e-8*IP''5; 

flag=^;  %  i.e.  will  )d.c)c  us  out  of  while  loop  and  bac3t  to  main  prog, 
end  %  IP  out  of  bounds  if  loop 
end 

t  _  III  program  "IHK/ITU" 
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i; 

i: 


'function  ''yprimes4"' 


I/:r»R  Bob  Stoney,  M>GS. 

Contains  the  equations  defining  conservation  of  linear  mcmentun, 

^  used  in  program  AWA-series 
'-i  Written;  Aug  92 
Last  ipdate:  4  May  93 

Tmprovernent  over  previous  editions: 

crarrma  dot  equation  (YDOr(4))  doesn’t  have  the  mall  angle  assunption  of 
%  previous  programs. 

^  See  AWAl  for  description  of  this  function! 

function  [yDOTl=yprimes4  (Y,THE17\,T,S,M,Q,CL,  CD, ALPHA)  %  [passed  out]  (passed  in) 
The  ydot  vector  is  as  follows :  (xdot,  hdDt,  vdot,  ganmadot) 

% 

%  This  function  calculates  the  y  solutions  to  the  coupled  system  equations: 
YDCrr(l)*Y(3)*cos(Y(4)); 

YDar(2)=Y(3)*sin(Y(4)); 

D=Q*S*CD;  %  Drag 

g=32.174;  %  gravitational  constant 

YDOT  ( 3) =T*oos (THETA-Y (4) ) /M-D/M-g*sin (Y  (4) ) ; 

L=g*S*CL;  %  Lift 

TEPM1=  (L-D*tan  (Y  (4) ) )  /  (M*Y  (3) )  ; 

TEPM2=T*  (tan(Y(4))*cos(ALPHA)+sin(ALPHA))/{M*Y(3)); 

'IT:RM3=-g/(Y(3)*cx)S(Y(4) )  ); 

TEBM4=-YD0T (3)  *tan  (Y  (4) )  /Y  (3)  ; 

YDOT  { 4 )  ==TEPm+TEPM2+TEPt43+’rePM4  ; 


=end  function  ”yprixnes4"' 
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%  I03R  Bob  Stoney.  Last  vjpdate;  12  May  93.  (vane  cxxitrol  power  re-rodeled  based 
%  on  engine  test  data) . 

%  This  function  is  for  aircraft  control.  It  is  essentiedly  a  model  of  Newton's 
%  second  law  ai^lied  to  oonservatlon  of  aingular  mcmentum  around  the  pitch  axis. 

%  lids  function  is  a  oonpanion  to  AWU.,  a  point  mass  simulator  for  the  Archytas 

%  Change  over  control  3:  Iigproved  Vane  control  power  model 
%  and  Cl  vs.  model. 

%  CAimcjN:  "Q"  is  d^tamic  pressure,  "q"  is  pitch  rate. 

function[THElADLrr,qout,q±jtout,Mcgw,M2gc,Mcgv,Mg,L,Q,CLl=ccncrol4. . . 

(q.DELlAT,THE17V,T,AIf’HA,y,Iyy,Q,S,Sc,Sv, . 

Ooar.Lc.Lv.Lw.BCW,  Dele,  Delv,  Cinq,  OCNIMXE)  ; 

% - Constant  pitch  rate  section: 

if  (XINIMX)E=1, 

THETACXnMlETA+q*DELTAT; %  Note:  for  0CNIMXE5=1  q,qdot  have  alreac^  been 

%  specified 

qout=q; 

cutout =0;  %  Constant  pitch  rate. 

Mcqw=0;  %  These  effects  not  included.  Passed  for  troubleshooting, 

Mcgc=d;  %  These  effects  not  included.  Passed  for  troubleshooting. 

Mcgv=0;  %  These  effects  not  included.  Passed  for  troubleshooting. 

Mq=0;  %  These  effects  not  included.  Passed  for  troubleshooting. 


% - Using  canard  and  vane  control: - 

else,  %i.e.  if  Dele  and  Delv  ^Decified 

%  — Wing  ocntributicxi; 

%  LIFT: 

%  Ffon  "Final  Archytas  Cl  vs.  ADA"  plot: 
if  (ALPHA<=-1.57)  |  (AIPHA>=1.57) , 

dispCWami.ig:  WING  aoa  greater  than  +/-  90  deg') 
elseif  AIfHA>-1.57  &  ALPHA<-0.72, 

CL=-1 . 083-0 . 69*ALPHA; 
elseif  AIi'HA>=>-0.72  &  ALPHA<=-0 . 285, 

CI;=2 . 81+33. 21  *AIPHA+105 . 38*ALPHA''2+1 38 . 1*ALPHA''  3+64 . 9*AIf>HA^4  ; 
elseif  ALf>»O»-0.285  &  AU>HA<=0.29, 

CLf=9 . 52e-2+3 . 38*AU>HA; 
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elseif  ALPHA>0.29  &  ALPHA<0.9, 

CL=^2 . 037+25 . 54*ALPHA-72 . 15*AlPHA''2+82 . 59*ALPHR''3-33 . 4»ALPHA^4; 
elseif  ALPHA>0.9  &  ALPHA<1.57, 

CL=1 . R75-1 . 194*ALPHA; 
end 


“’i  Qtiac  contriJbution:  (Ref:  NACA2412  in  Anderson.) 
QiBc=*-0 . 05; 

%  Add  Lift  and  Qmac  contributions: 


L=Q*S*CL; 

Mcgw=  -L*oos  (AIPHA)  *Lw  +  Onac*Q*S*Cbar; 

% - Canard  contribution;  %Presently  using  fit  frcm  main  wing  to  show  )vLcSi  aoa. 

%  From  "Final  Archytas  Cl  vs.  AOA"  plot: 

if  {{ALPHA-Delc)<=-1.57)  |  ( (ALPHA-Delc)  >=1 . 57) , 

Clc=0; 

ciisp( 'Warning:  Canard  aoa  greater  than  +/-  90  deg') 
elseif  (ALPHA-Delc) 5-1. 57  &  (AIPHA-Delc)<-0. 72, 

Clc=-1 . 083-0 . 69* (ALPHA-Delc) ; 
elseif  (ALPHA-Delc) >*-0.72  S  (ALPHA-Delc) <=>-0.285, 

Clc=2.81+33.21*(AIPKA-Delc)+105.38* (ALPHA-Delc) ^2. . . 

+138 . 1* (AIP>HA-Delc) ^3+64 . 9* (ALPHA-Delc) ^4; 
elseif  (ALPHA-Delc) >=-0.285  S  (AIi>HA-Delc)<=0.29, 

Clc=9. 52e-2+3. 38* (ALPHA-Delc) ; 
elseif  (AIPHA-Delc)>0.29  &  (ALPHA-Delc) <0.9, 

Clo=»-2. 037+25. 54*  (AlPHA-DelC)-72. 15* (ALPHA-Delc) "2 . . . 

+82 . 59*  (AIPHA-DelC)  ''3-33.4*  (ALPHA-Delc)  ''4; 
elseif  (ALPHA-Delc)  >=0.9  &  (ALPHA-^>elc)<1.57, 

Clc=l . 875-1 . 194*  (ALPHA-Delc) ; 
end  %end  if  for  Clc  loop 

I>fcgc=Q*Sc*Clc*Lc*cos  (ALPHA)  ; 

% - ^Vane  oontribution: 

%  Note;  "Delv",  the  vane  deflection,  is  assuned  to  be  0  at  0  control 
%  power.  Any  deviation,  plus  or  minus,  produces  a  pitching  mcment. 

Vanenvodel=2;  %  Used  to  ccnfiare  the  two  different  types  of 

%  vane  control  power  models.  l="old" — NACA0009  data  with 
%  velocity  measurements  from  engine  tests.  2="new" — mcment 
%  data  from  engine  tests  (only  applies  for  90-120  lb  but  left 
%  as  is  for  now)  converted  to  pitch  direction. 

%  NOTE:  Vanemodel=2  was  used  for  eill  runs  used  for  analysis  in  the  thesis 
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if  Vanenodel-=l, 

if  (Delv)<— .56, 

Clv=0; 

dispCWaming:  Vane  local  aoa  less  than  -.56  (low  end)') 
elseif  {(Delv)>-.56)  &  < (Delv)<=^.23) , 

Clv=l.  0419+22. 5*  (Delv)+67.754’*  (Delv)''2+56. 417*  (Delv)  "3; 
elseif  ((Delv))>-.23  &  ((Delv)<.2), 

Clv=2.2772e-20+5.8809* (Delv) ; 
elseif  (Delv) >=.2  &  (Delv) <.65, 

Clv«=-3 . 2096e-2+12 . 433*  (Delv) . . . 

-37.257* (Delv) "2+27 . 979* (Delv) "3; 
elseif  (Delv) >.65, 

Clv=0;  %***Beyond  specified  alpha  t)»re's  no  lift*** 
disp('Waming:Vane  local  aoa  gtr  .65  Oiigh  end)') 
end  %«yl  if  for  Civ  loop 

%  From  Thrust  tests,  the  following  relationships  apply  to  convert 
%  from  thrust  to  rpm  and  then  to  velocity  at  the  vane  (Vv)  (ft/sec) : 
BEM=1771 . 2+78 . 657*T-0. 26205*T"2; 

Vv^36.47-2.016e-3*PPM+1.867e-6*RPM"2;  %  f  (AQA,V?)  ***** 
Mcgv=-Lv*0 . 5*B0M*Vv"2*Sv*Clv;  %  Pitdiing  Mcment  due  to  vane 
elseif  Vanemodel=2, 

Newvanelift=2*  (0.267*Delv*180/pi) ;  %  must  convert  to  degrees 

%  'cause  that's  how  curve  was  fit. 

Mcgv=»-Newvanelift*Lv; 
end  %  if  for  vanerodel  type 

% - Pitch  dartping: 

^trQTC[*Q*S*abar"2/  (2*Iyy*Y  (3) )  ; 

% - Total :  Pitch  acceleration  (by  conservation  of  angular  mcmentm) : 

q±3tout=>  (Mcgw+Mcgc+Mcgv)  /Iyy+Mq[*q;  %q  dot  to  send  to  main  program 
qDut=q+qdDtout*DEL'IAT;  %q  to  send  to  main  program 


THEn7VDUT==TOETA+qput*DELTAT;  %pitch  angle  to  serxl  to  main  program 

end  %erid  if 

%  See  also:  AWA  assurptions  (Microsoft  word  document) 

^  ■  - 

I . .  UP  ocxrtrol  4  ~ 

-  - 
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'Function  "labeler"' 


This  function  provides  axis  labelling  for  AWA  series  of  programs, 
i  LCDR  Boh  Stoney,  Naval  Postgraduate  School 
Last  update:  12  NFVY  93  (addition  of  throttle  position) 

fui:cticn (XLABEL, YIABEL] =labeler (xanswer,yanswer) 

^  The  pass  in  list  contains  xanswer  and  yanswer,  which  is  what 
the  user  as)<ed  to  be  plotted  on  an  x-y  plot.  The  pass  out 
%  variables  XLABEL  and  YLABEL  are  text  strings  to  be  placed  on 
%  the  desired  plot, 
if  .xanswer==l, 

XlABEL=’Dcwnrange  distance,  x,  ft’; 
elseif  xanswer=2, 

XLABEL=' Height,  y,  ft'; 
elseif  xanswer=3, 

XLABEL^ ' Speed,  V,  f t / sec ' ; 
elseif  xanswer=4, 

XLABEL= 'Acceleration,  Vdot,  ft/sec^2'; 
elseif  .xanswer^S, 

>'1ABEL=' Pitch  attitude,  theta,  deg*; 
elseif  xa;iswer=6, 

XLABEL=' Angle  of  attack,  alp)ha,  deg’; 
elseif  .xanswer=7, 

XLABEL= ' Point  number  of  iteration’; 
elseif  y.answer=8, 

XLABEn>'Thurst,  T,  Ibf; 
elseif  xanswer=9, 

XLABEL= ' Grcundspeed,  ft /sec ' ; 
elseif  .xanswer=10, 

X!T'£EI^'Rate  of  Climb,  ft/sec'; 
els'll  f  ;-:answer=ll, 

XLA3EL='Flt  j3ath  angle,  gairma,  deg'; 
elseif  xan.swer=12, 

.XIABEL=' Energy  hei^t,  Es,  ft'; 
elseif  xanswer=13, 

XLABEL= ' Es  dot ,  f t/sec ' ; 
elseif  .xanswer=14, 

XLABEL=' Pitch  rate,  q,  deg/sec'; 
elseif  xanswer=15, 

XlABEL-'Pitch  acceleration,  qdot,  deg/sec^2’; 
elseif  .xanswer=16, 

XLAPEIe= ' Canard  deflection  (TE  down  +),Delc,  deg'; 
xanswer— 17, 

XIABEl-'Vane  deflection  (IE  up  +),Delv,  deg'; 
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e  1  i  f  xanswer=l  8 , 

XLABEL=’Tijre,  t,  sec'  ; 
elseif  xanswer=19, 

X1ABEL=’ Throttle  position,  TP,  percent’; 
end  %  end  if  locp  for  XLABEL 

. - Y  axis  label: - 

if  yanswer^l, 

YLAPEL='Dcwnrange  distance,  x,  ft'; 
elseif  yanswer=2, 

iXABEL=’ Height,  y,  ft'; 
elseif  yanswer=3, 

yiABEl^' Speed,  V,  ft/sec'; 

elseif  yanswer=4, 

yiABEL= 'Acceleration,  Vdot,  ft/sec''2'; 
elseif  yanswer=5, 

yiABEl/=' Pitch  attitude,  theta,  deg'; 
elseif  yanswer==6, 

yLABEL=' Angle  of  attack,  alpha,  deg'; 

; elseif  yanswer=7, 

YTABEL=' Point  nunber  of  iteration’; 
elseif  yanswer=8, 

YIABEL^'Thurst,  T,  Ibf; 
elseif  yanswer=9, 

ilABEL^'Groundspeed,  ft/sec'; 
elseif  yanswer=10, 

yiABEL='Rate  of  Climb,  ft/sec'; 
elseif  yanswer^ll, 

YLABEL='Flt  path  angle,  gamma,  deg'; 
elseif  yanswer=12, 

Y1ABEL=' Energy  heic^it,  Es,  ft'; 
elseif  yansv««r=13, 

YIABEL='Es  dot,  ft/sec'; 
elseif  yanswer=14, 

YLABEL='Pitch  rate,  q,  deg/sec'; 
elseif  yanswer=15, 

YlABEL='Pitch  acceleration,  qdot,  deg/sec''2'; 
elseif  vanswer==16, 

yiABEL= ' Canard  deflection  (TE  down  +),Delc,  deg'; 
elseif  yanswer==17, 

YlABEL='Vane  deflection  CIE  +),Delv,  deg'; 
elseif  yanswer=18, 

yiABEL='Tijne,  t,  sec'; 
elseif  yanswer=19, 

YLABEL=' Throttle  position,  TP,  percent'; 
end  %  end  if  locp  for  XLABEI, 


^  .  r^nd  function  "labeler"- 
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■ 

■-  - - =inisc.  program  "plotter"' 


^  Program  used  in  conjunction  with  the  AKA  series  programs. 

^  Uses  subplot  to  plot  most  Important  parameters  on  a  two 
*  sheets  of  paper. 

%  To  use,  siJTply  typie  plotter  at  the  MATIAB  ccnmand  carrot  » 

!=  LCDR  Bob  Stoney,  NPGS 
Created  12  SEP  92 
clg 

hold  off 

^ - Aircraft  parameters;  Alpha,  speed,  theta,  q - 

subplot  (221) , plot  (CXJTFIUE  ( : ,  18)  ,CXnriLE( : ,  6) ) 

xlatjel  (’Time,  t,  sec' )  ;ylabel  ( 'Angle  of  attac)c,  alpha,  deg'), -grid 

title (( 'Weiqht=' , num2str  (Wt) , '  Ibf  ')) 

snip]  Ml  (22?) ,  plot  ((XnrilJ!:  ( : ,  18)  ,OinTlLE  ( ; ,  3) ) 

xlabel  ('Time,  t,  sec ');yla]oel(' Speed,  V,  ft/sec' ) ;grid 

subplot  (223) ,  plot  (OOTFIIZ  ( : ,  1 8)  .OUTFUE  ( : .  5) ) 

xlatjel ( 'Time,  t,  sec' ) ;ylabel ( 'Pitch  attitude,  theta,  deg’) /grid 

subplot  (224) ,  plot  (Oirn:TI£  ( : ,  18)  ,<3JTFILE  ( : ,  14) ) 

xlabeK'Time,  t,  sec' ) /ylabel ('Pitch  rate,  q,  deg/sec' ) /grid 

pause 

cla 

nold  off 

^ - Ctontrol  parameters:  Canard,  vane,  thrust,  and  x  vs.  h - 

subplot (221 ) , plot (OOTFILE ( ; , 18) , OUTFILE ( : ,  16) ) 

xlabeK'Time,  t,  sec' ) /ylabel ('Canard  defl..  Dele,  deg’)/grid 

title (( 'Weight=', num2str (Wt), '  Ibf'J) 

subplot  (222) ,  plot  (OnriLE  ( : ,  18) ,  OUTFILE  ( : ,  17) ) 

xlabeK'Time,  t,  sec' ) /ylabel ('Vane  deflection,  Delv,  deg’)/grid 

subplot (223) ,  plot (OUTFILE ( : , 18) , OUTFILE ( : , 19) ) 

xlatjeK'Time,  t,  sec' ) /ylabel ( 'Throttle  posit,  TP,  %')/grid 

subplot  (224),  plot  (OUTFILE  ( :,1),  OUTFILE  (:,  2)) 

xlabel  ( ' Downrange  distance,  x,  ft ' )  / ylabel  ( '  Hei^t,  h,  f t ' )  /  grid 

pause 

hold  off 


%■  -  ■  -L  ■■  ■  -  g— — 

■  ■  — ■  -  -  end  program  "plotter"' 

^ -  ■  —  -  _ _ 
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Lsc.  program  "Movie' 


'Movie'  LCDR  Bob  Stoney,  NFS.  1  Nferch  93 
last  update;  30  May  93 

'  Ccmpanion  program  to  AWAl  that  shows  the  vehicle's  position  and  pitch 
'  at’ itude  during  the  most  reoent  run. 

clg 

otf 

ansl=size  (KEEPmRK)  ; 
ansla=ansl (1) ; 

dispCWhat  axes  would  you  like  on  the  movie?  [xmin,  xmax,  ymin,  ymax]  ') 
dispC Enter  1  for  default  to  last  axis  used  on  AWAl  run.') 
mo'yiea.xis=input  ( '  ' ) ; 
if  movieaxis=l, 

axis(V);  %  i.e.  DEFAULT 

else, 

a.xis  ( [movieaxis(l)  ,mcfvieaxis(2)  .movieaxis  (3)  ,movieaxis(4)  ] ) ; 

%  Axis  limits  ([xmin,  xnax,  ymin,  ymax]) 
ei«,i 

plot  (0LnFILE(l,l),0inFILE(l,2),  'o') 

xlabel CDownrange  distance,  x,  ft' ) ;ylabel(’ Altitude,  h,  ft') 

Ll=sprintf ('Delta  T  =  %g',DELTAT) ;text(.2, ,85,L1, 'sc' ) 

L2=sprintf ('Weight  =  %g',Wt);text(.2, .8,L2, 'sc') 

STEPSIZE=5;  %  Interval  of  plotting  the  vehicle's  position  and  pitch  attitude 
for  PT=1 : STEPS! ZE:ansla; 

^  Scale  the  pitch  attitude  to  the  last  scale  used  by  main  program 
xScaler= (Vkeep (ansla, 2) -Vkeep (ansla, 1) ) / (Vkeep (PT, 2) -Vkeep (PT,  1 ) ) ; 

%to  scale  pitch  attitude  line  in  x  direction  to  last  plot  coordinates 
dx=KEEPdxdy (PT, 1 ) *XScaler; 

YScaler= (Vkeep (ansla, 4) -Vkeep  (ansla, 3) ) / (Vkeep (PT, 4) -Vkeep (PT,  3) ) ; 

%to  scale  pitch  attitude  line  in  y  direction  to  last  plot  coordinates 
d.-=KEEPdxdy  (PT,  2)  *YScaler; 
end  scaling 

polymark  (CXJrFII£  (PT,  1)  ,0UTFILE{PT,2) ,  'o' ) 
polyline  ([CWTFTLE(PT,l)-dx  CXm’ILE(PT,l)+dx]  . . . 

.  (an r'lLE  (PT,  2)  OUTFILE (PT,  2)  +dy] , '  -w' ) 

^  This  -Viows  the  vehicle's  pitch  attitude 
end 


=end  program  "Movie"’ 
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% - Program  ^^lot - 

%  Program  to  aid  in  troubleshooting  the  control  function. 

Simply  generates  a  plot  of  the  Manents  about  the  eg 
%  due  to  Wing,  canard,  vane  and  pitch  datiping  vs.  time. 

%  Also  shows  {on  another  plot)  the  theta,  q  and  qdot  results  of 
the  control  function. 

To  use:  type  "Nplot"  at  the  cormand  carrot. 

%  LCT)R  Bob  Stoney 
%  Created  13  SEP  92 
clg 

hold  off 
sutplot(211) 

plot(OUTFILE(;,18),EXTRA(:,ll),  '-',OinFILE(:,18),EXTRA(:,12),  '  — - 

OOTFILE ( : , 18) , EXTRA ( : , 13) , ' ,OUTFILE ( : , 18) , EXTRA ( ; ,  14) ,  ’ o’ ) 
xlabel  ('Time, t, sec' )  ;ylabel  ('NWing-  C — V*  qoo’);grid 
title (( 'Weight®', num2str (Wt), '  Ibf  ')) 
subplot  (212) 

plot  (OUTFILE  ( : ,  18) ,  OUTFILE  ( : ,  5) , '  -  ’ ,  dTITTLE  ( : ,  18) ,  CUTFILE  (:,14),’  — ' . 

OUTFILE  ( : ,  18) ,  (XTTFILE  ( : ,  15) , ' : ' ) 
xlabel ('Time, t, sec' ) ;ylabel ('theta-  q —  cjdot. . ’) ;grid 
title (( 'Weight® ', nir\2str (Wt ), '  Ibf']) 
pause 

%  — ■  -  — -  —  ■■  -■  .  -  - 

% -  -  ■•■  end  program  "^plot"=  -- 
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^  "filesaver". .  .this  program  is  used  to  save  information  frcm 
program  AWAl. 

LCDR  R.B.  Stoney,  Naval  Postgraduate  School 
>  May  1993 

R27=arrFILE; 

R27E=E<TRA; 

R27C=OONFIG; 
save  B27  b27; 
save  B27E  R27E; 
save  R27C  R27C; 
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i  "resurrect". . . .a  program  to  load  the  results  of  a  run  of  program 
PWAl  that  were  saved.  See  also  "filesaver” 

«  l£DR  R.B.  Stoney,  Naval  Postgraduate  School 
i  May  1993 

^  user:  change  the  norioer  on  each  line  for  the  appropriate  file 
Note;  can  use  "save  workspace"  if  desired  but  that  takes  up  a  LOT 
%  of  memory 


load  R24  %  "OUTFILE"  matrix 

load  R24C  %  "COUfTiG"  matrix 
load  R24E  %  the  "EXTRA"  matrix 
OUTFILE=R24; 

Wt=R24C(l); 

dispCFile  ntmber  R24  has  been  recovered') 


%  •  end  program  "resurrect"^ 
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Throttle  posit,  TP.Qinard  deflection.  Dele,  deg  Pitch  attitude,  theta,  cWgtgle  of  attack,  alpha,  deg 
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Figure  H.l 

Archytas  simulation  "Low  vane  deflection  pushover  from  hover" 


Appendix 


Throttle  posit,  TP,®inard  deflection.  Dele,  deg  Pitch  attitude,  theta,  ddgigle  of  attack,  alpha,  deg 


Throttle  posit,  TP.^nard  deflection.  Dele,  deg  Pitch  attitude,  theta,  d^gle  of  attack,  alpha,  de 
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Figure  H.3 

Archytas  simulation;  "Low  speed,  low  altitude  climbing  gentle 

pushover" 
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Throttle  posit,  TP,®inard  deflection.  Dele,  deg  Pitch  attitude,  theta,  cWgigle  of  attack,  alpha,  deg 
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Figure  H.4 

Archytas  simulation:  "Vertical  acceleration  to  gentle  pushover" 
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Throttle  posit,  TP,  %  Canard  defl..  Dele,  deg  Pitch  attitude,  theta,  degAngle  of  attack,  alpha,  deg 
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Figure  H.5 

Archytas  simulation:  maximum  effort  slow  flight,  unsuccessful 

attempt 

at  level  transition  from  horizontal  flight  to  hove 
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Throttle  posit,  TP,  %  Canard  defl..  Dele,  deg  Pitch  attitude,  theta,  deg  Angle  of  attack,  alpha,  deg 
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Figure  H.6 

Archytas  simulation:  Transition  to  horizontal  flight 
followed  by  "zooming"  transition  back  to  hover 
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